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PREFACE 


The  Mississippi  Sound  and  Adjacent  Areas  Study  was  authorized  by  congres¬ 
sional  resolutions  of  1  February  1977  and  10  May  1977.  The  main  purpose  of 
the  study  as  stated  in  the  resolutions  is  to  determine  whether  the  present  and 
proposed  dredged  material  disposal  methods  for  maintenance  and  construction  of 
the  various  projects  in  Mississippi  Sound  should  be  modified  in  any  way  in  the 
interest  of  economic  efficiency  and  environmental  quality.  The  resolutions 
request  an  investigation  by  the  US  Arm'  Corps  of  Engineers  of  various  dredging 
techniques  and  the  possibility  of  developing  a  coordinated  program  for  the  re¬ 
gion,  with  appropriate  consideration  of  ecological  factors.  The  region  under 
study  is  defined  to  include  the  body  of  water  and  adjacent  land  and  estuarine 
areas  extending  from  Chandeleur  Sound  and  Lake  Borgne,  on  the  west,  along  the 
Mississippi  and  Alabama  coasts  to  the  eastern  shore  of  Mobile  Bay,  on  the  east. 
It  is  bounded  on  the  north  by  Interstate  Highway  10  and  on  the  south  by  the 
120-ft  bottom  contour  of  the  Gulf  of  Mexico. 

The  numerical  model  investigation  described  herein  was  authorized  by  the 
US  Army  Engineer  District,  Mobile  (SAM).  This  study  was  conducted  at  the 
US  Army  Engineer  Waterways  Experiment  Station  (WES)  in  the  Wave  Dynamics 
Division  (WDD) ,  Hydraulics  Laboratory,  under  the  direction  of  Mr.  H.  B. 

Simmons,  Chief  of  the  Hydraulics  Laboratory,  Dr.  R.  W.  Whalin,  former  Chief 
of  Che  Wave  Dynamics  Division,  and  Mr.  C.  E.  Chatham,  Jr.,  present  Chief  of 
the  Wave  Dynamics  Division. 

The  investigation  was  performed  and  this  report  prepared  by  Dr.  R.  A. 
Schmalz,  WDD.  Ms.  Mary  Ann  Leggett  retuned  the  Gulf  Tide  Model  in  order  to 
provide  revised  tidal  constituent  information.  Ms.  J.  I.  Jones  prepared  the 
geophysical  data  and  developed  Appendix  A. 

The  cooperation  of  and  coordination  with  SAM  personnel ,  including 
Messrs.  Dennis  McCann,  Maurice  James,  Dru  Barrineau,  and  Timothy  Phillips, 
enabled  the  work  to  remain  focused  on  district  needs  and  provided  for  an 
effective  transfer  of  the  technology  developed.  The  numerical  computations 
associated  with  this  work  were  performed  on  CYBER  175  and  CRAY  1  computers 
located  at  the  Air  Force  Weapons  Laboratory,  Kirtland  Air  Force  Base,  New 
Mexico.  The  numerical  model  and  all  datasets  were  transferred  to  SAM. 

WES  provided  assistance  in  installing  the  model  on  the  Boeing  Computer  System 
for  use  by  SAM. 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can 
(SI)  units  as  follows: 


Multiply 

By 

cubic  feet  per  second 

0.02831685 

feet 

0.3048 

feet  per  second 

0.3048 

miles  (US  statute) 

1.852 

square  feet 

0.09290304 

square  feet  per  second 

0.09290304 

square  miles  (US  statute) 

0.4470400 

be  converted  to  metric 

_ To  Obtain _ 

cubic  metres  per  second 
metres 

metres  per  second 

kilometres 

square  metres 

square  metres  per  second 

square  metres 


NUMERICAL  MODEL  INVESTIGATION  OF  MISSISSIPPI  SOUND 


AND  ADJACENT  AREAS 


PART  I:  INTRODUCTION 


1.  Mississippi  Sound  and  its  adjacent  waters  as  shown  in  Figure  I-I 
comprise  an  extremely  productive  aquatic  ecosystem.  Major  navigation  channels 
maintained  by  dredge  operations  within  Mississippi  Sound  are  located  at  Gulf¬ 
port,  Biloxi,  and  Pascagoula.  In  order  to  protect  the  ecosystem  as  well  as 
provide  for  efficient  navigation  a  study  was  authorized  by  congressional  reso¬ 
lutions  of  the  Senate  Committee  on  Environment  and  Public  Works  and  the  House 
Committee  on  Public  Works  and  Transportation.  The  US  Army  Engineer  District, 
Mobile,  has  specific  authority  and  jurisdiction  to  conduct  the  study. 

2.  A  numerical  modeling  approach  is  required  to  develop  in  an  efficient 
manner  a  quantitative  knowledge  of  tidal  circulation  and  salinity  distribution. 
The  US  Army  Engineer  Waterways  Experiment  Station  (WES)  was  requested  to  apply 
WIFM,  a  two-dimensional  vertically  integrated  model  (Butler  1980).  In  addi¬ 
tion,  model  capability  was  to  be  extended  to  include  the  prediction  of  salin¬ 
ity.  In  order  to  calibrate  and  verify  the  model,  Raytheon  Ocean  Systems  was 
contracted  by  the  Mobile  District  to  collect  prototype  velocity,  temperature, 
conductivity,  and  meteorological  data.  The  National  Oceanic  and  Atmospheric 
Administration  (NOAA)  was  responsible  for  obtaining  tidal  elevation  data. 

Texas  A&M  was  contracted  to  develop  a  numerical  Gulf  Tide  Model  (GTM)  to  pro¬ 
vide  accurate  tidal  elevation  information  at  the  boundary  of  the  WIFM  model 
global  grid. 

3.  This  report  is  structured  in  the  following  manner. 

a.  In  Part  II,  both  the  Raytheon  Ocean  System  and  NOAA  Data  Collec¬ 

tion  Programs  features  are  presented.  Harmonic  analysis  results 
are  presented  in  Part  III  for  both  water  surface  elevations  and 
currents.  The  GTM  is  outlined  and  the  results  for  the  0^  ,  , 

Pj  ,  M^  ,  and  S^  tidal  constituents  are  presented  in  Part  IV. 

Major  features  of  the  salinity  algorithm  incorporated  within 
WIFM  are  discussed  in  Part  V'  in  terms  of  the  equations  con¬ 
sidered,  the  numerical  approximations,  and  the  effective  disper¬ 
sion  coefficient  formulation. 

b.  The  remaining  parts  of  this  report  comprise  the  numerical  model¬ 
ing  effort  and  are  presented  in  the  order  in  which  this  work  was 
performed . 
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ssippi  Sound 


Initially,  a  global  grid  representing  Mississippi  Sound  and  ad¬ 
jacent  areas  was  constructed  to  interface  with  the  GTM  grid  en¬ 
compassing  the  entire  Gulf  of  Mexico.  The  development  of  the 
global  grid  depth  field,  barrier  island  conf iguration ,  fresh¬ 
water  inputs,  and  location  of  the  gaging  stations  is  presented 
in  Part  VI . 

In  order  to  determine  periods  for  study  on  the  global  grid,  data 
over  the  entire  180-day  survey  period  were  reviewed.  Part  VII 
outlines  the  data  review  in  terms  of  both  hydrodynamic  and  salin¬ 
ity  simulation  requirements.  Water  surface  elevation  and  cur¬ 
rent  tidal  constituent  constants  developed  in  the  harmonic 
analysis  represent  astronomic  influences  only;  i.e.,  all  meteoro¬ 
logical  effects  are  assumed  to  be  removed  in  the  filtering  pro¬ 
cess  employed  in  the  analysis.  In  contrast,  no  harmonic  analysis 
was  conducted  for  measured  salinity  values  and,  as  a  result,  the 
meteorological  effects  have  not  been  removed.  The  meteorological 
effects  are  assumed  to  consist  of  wind  only;  pressure  anomalies 
were  not  considered.  Therefore,  in  the  simulation  of  salinity, 
the  wind  distribution  over  the  grid  must  also  be  specified. 

Since  this  represents  an  additional  input,  salinity  simulations 
were  considered  after  all  hydrodynamic  simulation  work  had  been 
completed.  The  hydrodynamics  were  considered  first  on  the  global 
and  subsequently  on  a  refined  grid  encompassing  the  Pascagoula 
Channel . 

In  Part  VIII  the  20-25  September  1980  period  is  studied  in  detail 
in  order  to  identify  a  suitable  time  for  calibration.  The  first 
five  days  (20-24  September)  were  selected  as  the  calibration  pe¬ 
riod.  In  Part  IX  detailed  data  for  the  period  12-16  June  1980 
are  presented  as  a  separate  verification  period. 

Global  grid  hydrodynamics  are  presented  in  Part  X.  Bottom  fric¬ 
tion  mechanics  are  calibrated  for  the  20-24  September  period  and 
subsequently  verified  for  the  12-16  June  1980  period.  A  hypo¬ 
thetical  regional  dredge  plan  involving  increasing  the  size  of 
Sand  Island  was  considered  by  modifying  the  depth  field  in  the 
vicinity  of  this  feature  in  the  global  grid.  Hydrodynamics  for 
the  20-24  September  1980  period  were  simulated.  By  comparing 
the  flow  structure  with  the  Sand  Island  complex  in  place  with 
the  original  flow  structure  it  was  possible  to  dirpctly  deter¬ 
mine  the  extent  of  influence  of  this  type  of  change  within  the 
global  grid. 

The  methodology  to  be  employed  in  constructing  refined  grids 
around  alternative  navigation  channel  projects  is  developed  in 
Part  XI  along  with  the  development  of  a  refined  grid  around  the 
Pascagoula  Ship  Channel. 

Refined  grid  hydrodynamics  are  presented  in  Part  XII.  Bottom 
friction  mechanics  calibrated  and  verified  on  the  global  grid 
were  employed  in  a  simulation  of  the  20-24  September  1980  period. 
Water  surface  elevations  developed  in  the  global  grid  were  used 
to  drive  the  refined  grid  boundary. 
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i.  All  salinity  simulation  work  is  presented  in  Part  XIII.  The  20- 
24  September  1980  period  was  considered.  Wind  speed  and  direc¬ 
tion  were  specified  as  meteorological  inputs.  A  constant  drag 
coefficient  of  0.001  was  employed  to  develop  the  surface  wind 
stress,  and  the  friction  mechanism  previously  developed  was  used 
to  implement  bottom  stress.  Simulation  results  for  both  the 
Flux-Corrected  Transport  and  the  Three  Time  Level  Explicit 
schemes  are  presented  for  the  global  grid.  Effective  dispersion 
coefficients  in  each  coordinate  direction  were  calibrated.  The 
Flux-Corrected  Transport  global  grid  salinity  levels  and  water 
surface  elevation  saved  at  the  boundary  of  the  refined  grid  were 
used  to  supply  the  boundary  conditions  for  a  Flux-Corrected 
Transport  simulation  on  the  refined  grid.  Previously  calibrated 
global  grid  effective  dispersion  c  ue  f  f  i  *.  i  eul  s  were  used  in  the 
refined  grid  simulation. 

j.  Study  conclusions  and  recommendations  for  additional  simulation 
work  are  presented  in  Part  XIV.  The  procedures  for  constructing 
the  input  dataset  for  the  Pascagoula  Channel  refined  grid  are 
outlined  in  Appendix  A.  This  appendix  is  intended  to  provide  a 
guide  for  District  use  in  developing  refined  grids  encompassing 
the  Biloxi  and  Gulfport  Ship  Channels.  A  general  tide  genera¬ 
tion  program  is  documented  in  Appendix  B.  This  program  was  used 
to  generate  the  predicted  tide  based  on  the  five  major  con¬ 
stituents  and  compare  the  predicted  tide  levels  with  the  unfil¬ 
tered  and  filtered  water  levels  and  currents  at  each  gaging  sta¬ 
tion.  In  addition,  the  tide  generation  portion  of  the  program 
was  incorporated  within  WIFM  in  order  to  generate  water  surface 
elevations  at  the  boundary  of  the  global  grid.  A  cubic  poly- 
nominal  feathering  technique  is  presented  in  Appendix  C  in  order 
to  smoothly  transition  from  the  zero  elevation  and  velocity 
state  to  the  predicted  tide  levels  on  the  boundary.  In  order  to 
efficiently  specify  linearized  advection  conditions  around  bar¬ 
rier  islands,  subroutine  ADVBAK  was  developed  as  documented  in 
Appendix  D  and  incorporated  within  WIFM. 


4.  In  the  numerical  work  presented  in  this  report  a  nested  grid  model¬ 
ling  philosophy  is  utilized.  Simulation  results  over  the  outer  nest  are  used 
to  provide  the  boundary  conditions  for  the  next  inner  nest.  The  nesting  is 
taken  to  level  two  ( GTM  drives  WIFM  global  grid,  which  drives  WIFM  refined 
grid).  Inherent  in  this  approach  is  the  premise  that  changes  made  on  the 
inner  grid  nest  do  not  cause  hydrodynamic  and/or  salinity  fluctuations  which 
propagate  to  the  boundary  of  the  inner  grid  nest,  thereby  destroying  the  in¬ 
tegrity  of  boundary  values.  As  part  of  t  ti  l  s  study,  in  order  to  demonstrate 
the  validity  ol  this  nested  grid  approach,  the  region  of  extent  of  fluctua¬ 
tions  in  hydrodynamics  and  salinity  generated  by  changes  ill  the  global  grid 
and  refined  grids  was  investigated.  Prior  to  t h i s  study,  the  region  of  extent 
of  navigation  channel  induced  changes  on  flow  si  rue  line  was  not  known.  Thus 
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it  was  not  known  whether  contemplated  deepening  and  widening  alterations  for 
the  Biloxi  Ship  Channel  might  influence  Pascagoula  Channel  conditions  and  vice 
versa.  A  major  concern  was  that  to  study  any  one  single  channel  alteration  an 
extremely  large  refined  grid  encompassing  all  navigation  channels  would  need 
to  be  employed.  The  numerical  work  presented  was  focused  to  address  these 


issues . 


PART  II:  DATA  COLLECTION  PROGRAM  ESTABLISHMENT 


5.  WES  Wave  Dynamics  Division  personnel  participated  in  a  series  of 
meetings  held  at  WES  with  the  Mobile  District  and  their  consultants  to  develop 
prototype  data  collection  requirements  for  calibration  and  verification  of  the 
numerical  model.  Preliminary  meetings  were  conducted  in  March  and  April  1979 
prior  to  formal  study  authorization.  In  March  1980,  station  locations  and 
data  collection  procedures  were  finalized. 

Raytheon  Ocean  Systems  Program 

6.  Raytheon  Ocean  Systems  collected  oceanographic  and  meteorological 
data  over  the  period  23  April  through  20  October  1980.  Approximately  40  cur¬ 
rent  meters  and  18  conductivity-temperature  instruments  were  deployed  among 
the  21  stations  shown  in  Figure  1 1 —  1 .  Wind  speed  and  direction,  air  tempera¬ 
ture,  and  air  pressure  were  measured  at  the  five  meteorological  stations  shown 
in  Figure  II-2.  Bottom  pressure  measurements  were  recorded  at  the  three  deep- 
sea  sites  shown  in  Figure  11-3.  Salinity  and  temperature  depth  profiles  were 
acquired  at  three-week  intervals  for  the  stations  shown  in  Figure  II-4.  Hydro- 
graphic  surveys  of  the  barrier  island  passes,  major  navigation  channels,  and 
bay  passes  as  shown  in  Figure  I I —5  were  conducted  during  the  period  23  June 
through  11  July  1980. 

7.  Deployment  procedures,  quality  control,  data  reduction,  instrument 
specifications,  and  data  return  are  as  reported  by  Raytheon  Ocean  Systems 
(1981) . 


NOAA  Data  Collection  Program 

8.  NOAA  and  their  contractors  obtained  tidal  elevation  data  over  the 
Apri 1-November  1980  period  at  the  stations  shown  in  Figure  1 1  —  3 .  Six-minute 
and  hourly  data  were  edited  and  furnished  to  WES  on  magnetic  tape.  Tide  gage 
reference  levels  were  connected  to  NGVD  (1929)  to  provide  for  a  consistent 
geodetic  reference.  The  relations  between  mean  lower  low  water  (MLLW)  and 
NGVD  (1929)  as  well  as  the  exact  station  locations  are  shown  in  Table  1 1  —  1 . 
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Table  1 1-  1 

NOAA  Tide  Station  Locations 


Station 

No . 

Mean  Lower 
Low  Water 
ft* 

873-1269 

-0.29 

873-1952 

-0.21 

873-5184 

(-0.29)** 

873-5523 

-0.24 

873-5587 

-0.40 

873-7048 

(-0.32)** 

874-0199 

-0.24 

874-0405 

-0.26 

874-1196 

-0.30 

874-2221 

-0.29 

874-3081 

-0.44 

874-3495 

(-0.40)** 

874-3735 

-0.36 

874-4586 

-0.42 

874-467  1 

(-0.40 )** 

874-4756 

-0.28 

874-6819 

-0.39 

874-7437 

-0.27 

874-9704 

(-0.27)** 

8  7  6  -  0  4  1  2 

(-0.28)** 

876-0596 

(-0.28)** 

876-0742 

(-0.28)** 

Latitude 

(N) 

Longitude 

(W) 

30°  14 

.9 

87°40 , 

.  1 

30° 18 . 

.2 

87°44 . 

1 

30°  15 . 

.0 

88°04 . 

5 

30°26 . 

.6 

88°06 . 

8 

30°  1 5 

.5 

88°06 . 

8 

30°42 

.3 

88°02 

4 

30°20 . 

.5 

88°24 . 

4 

30° 1 2 , 

.2 

88°26 . 

5 

30°20 . 

.4 

88°32 . 

0 

30° 14. 

.  1 

88° 39. 

7 

30°23 . 

.2 

88°46 . 

4 

30°25 . 

.2 

88°49 . 

7 

30°23 . 

.4 

88°  5 1 

4 

30°23 

.3 

88°57 . 

.8 

30°24 

.8 

88°58 . 

.5 

30° 12. 

.  7 

88°58 . 

.3 

30° 18 . 

.6 

88° 14 . 

,  7 

30° 19 , 

.5 

89° 1 9 . 

.5 

30°  14 

.  7 

89° 36 . 

.8 

29°  1  2 

.  3 

89°02 . 

.  2 

29°29 , 

,  6 

89° 10 . 

4 

29°49 . 

,4 

89° 16. 

,2 

_  Location  _ 

Gulf  Shores,  AL 
Bon  Secour,  AL 
Dauphin  Island,  AL 
Fowl  River,  AL 
North  Poi nt ,  AL 
Mobile,  AL 

Grand  Batture  Island,  MS 

Petit  Bois  Island,  MS 

Pascagoula,  MS 

Horn  Island,  MS 

Davis  Bayou,  MS 

Old  Fort  Bayou,  MS 

Cadet  Point,  MS 

Broadwater  Marina,  MS 

Popps  Ferry  Bridge,  MS 

Ship  Island,  MS 

Pass  Christian,  MS 

Bay  Waveland,  MS 

Pearl  River  at 
Burl i ngton ,  MS 

Nort  h  Pass ,  LA  + 

Breton  Island,  LA t 

Comfort  Island,  LA 


MLLW  referenced  to  NGVD. 
Approximated  from  nearby  stations. 
Not  shown  in  Figure  I  I -3. 


PART  III:  PROTOTYPE  DATA  ANALYSIS 


9.  A  harmonic  analysis  was  performed  on  both  water  surface  elevation 
and  current  data.  The  analysis  techniques  are  presented  in  detail  by  Outlaw 
(1983).  We  outline  briefly  the  techniques  and  results  of  the  analysis  for 
water  surface  elevations  and  currents  in  turn  below. 

Water  Surface  Elevations 

10.  The  harmonic  analysis  for  surface  elevation  tidal  constituents  was 
conducted  for  22  NOAA  tidal  elevation  stations  and  three  Raytheon  pressure 
cells  employing  hourly  data.  Station  numbers  and  locations  are  as  shown  in 
Table  1 1  —  1  .  The  NOAA  stations  are  representative  of  tidal  elevations  in  the 
nearshore  region  and  along  the  Mississippi  Sound  barrier  islands.  The  pres¬ 
sure  cell  data  are  representative  of  the  tide  in  the  Gulf  approximately 

27  miles*  south  of  the  barrier  islands. 

11.  The  harmonic  analysis  included: 

a.  Editing  to  remove  data  spikes,  insertion  of  missing  data,  and 
subtraction  of  the  mean  from  the  data  record. 

b.  Filtering  to  remove  high  and  low  frequency  trends  from  the  data. 

c.  Harmonic  analysis  for  tidal  constituents. 

The  NOAA  surface  elevation  data  had  been  edited  by  NOAA  and  was  continuous 
over  the  record  length.  During  the  analysis,  pressure  cell  data  were  con¬ 
verted  to  surface  elevation  and  corrected  for  barometric  pressure  changes 
prior  to  removing  the  mean. 

12.  A  digital  band-pass  filter  was  applied  to  eliminate  frequencies  in 
the  data  outside  the  semidiurnal  to  diurnal  tidal  frequency  range. 

13.  In  the  harmonic  analysis,  the  tidal  elevation  at  a  given  station  is 
represented  as  follows: 

h(t)  =  H  +  ^  f.H.  cos  +  (V  +  u).  -  K.) 

°  ^  1  1  \  T  i  °  1  7 

i=l 

A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  12. 
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where 


The 


h  ( t  )  = 
t  = 

H  = 

o 

f  = 

1 

l  V  +  u )  = 

o  1 

I  - 
i 

n  - 
H  = 

i 

K  - 
1 

const i tuent s 


tidal  elevation  as  a  function  of  time 
time 

mean  of  the  filtered  data  record 
node  factor  for  constituent  i 
ecju  i  1  1 1>  r  l  uni  argument  for  constituent  i 
period  tor  constituent  i 
number  of  constituents  considered 
mean  amplitude  for  constituent  i 
local  epoch  for  constituent  1 
included  in  the  analysis  and  their  periods 


are : 


Harmonic  Constituent 


Symbo 1 


Period,  Ji^r 


Principal  lunar  diurnal  01 

Lunisolar  diurnal  K1 

•  Principal  solar  diurnal  PI 

Smaller  lunar  elliptic  Ml 

Small  lunar  elliptic  J1 

Larger  lunar  elliptic  Q1 

Principal  lunar  M2 

,  Principal  solar  S2 

|  Larger  lunar  elliptic  N2 


25.82 
23.94 
24.07 
24.84 
23.  10 
26.87 
12.42 
12.00 
12.66 


14.  The  mean  amplitude,  H  ,  and  local  epoch,  K.  are  determined  by 
minimizing  the  variance  between  the  filtered  record  and  the  predicted  record 


given  in  the  equation  above. 

15.  Mean  elevations  of  the  NOAA  Data  Record  relative  to  mean  lower  low 


water  are  given  in  Table  1 1 1-1.  The  constituent  amplitude  and  phases  are  pre¬ 
sented  in  Tables  I  1 1  - 2  and  II1-3,  respectively.  The  length  of  the  data  record 
and  the  root  mean  square  (RMS)  error  are  given  in  Table  111-4.  The  principal 
tidal  constituents  are  the  diurnal  constituents  01  ,  K1  ,  PI  ,  and  the  semi¬ 

diurnal  constituents  M2  and  S2 


Cu r  rent s 

16.  Till*  harmonic  analysis  foi  currents  was  .  otidui  t  ed  f"t  all  observed 
data  separately  for  the  east-west  and  u<  •  r  t  h  -  sou  t  fi  t omponenl s  In  t  h  i  s  manner 
analyzed  current  components  correspond  in  orientation  to  the  uumerit.il  model 
current  components  thereby  f.uilitating  i  ompa  i  i  .  n  s 


Table  1 1 1  -  1 


Mean  E levation  of  Observed  NOAA  Data  Record 
Relative  to  Mean  Lower  Low  Water 


Stat  ion'" 

Eleva  t i on 

873-1269 

0.71 

873-1952 

0.69 

873-5184 

-- 

873-5523 

0.73 

873-5587 

0.96 

873-7048 

0.85 

874-0199 

0.85 

874-0405 

0.88 

874-1196 

0.88 

874-2221 

0.90 

874-3081 

1.07 

874-3495 

-- 

874-3735 

1 .01 

874-4586 

1 .04 

874-4671 

-- 

874-4756 

0.90 

874-6819 

0.99 

874-7437 

0.87 

874-9704 

-- 

876-0412 

0.56 

876-0595 

0.71 

876-0742 

0.92 

Refer  to  Table  I  I  - 1  for  T  designations 
used  in  1 i gu re  I  I  -  3  . 
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tituent 


designations  used  in  Figure  II- 


Elevation  Local 


1  ■  CO  '£'  f  "J  ''-I  _■  7.  rf  _•  'ij  £1  r  _  ,j-|  .Vi  f  .V.  -|J  ,v,  T, 

j-i  —•  ,t,  r-  •:>:>  .7.  y.  e  t  l  >j-.  .7.  — .  £1  it.  t  0*1  '7-  — •  :o 

CM  ^  |V(  *tj  1  .'.*1  Lj'l  t  f  _  y“  rf  !j-.  jj*|  f'A  |V,  ,"lj  .—I  y  J*.  Tj.  £1  ',J  .‘ij  ,V. 

*  i-.-i  .V,  ,v,  C.J  AA  ;V,  v,  AA  —  'O 


'0  .t.  iT'  1"  v  •£•  •-*  — 1  i'iJ  -•  t  !T'  ■  j'1  j'.'  -•  <t  .;o  -:o  y  .'|J  7. 

cm  ^  J,  y  J-  .7.  y  cij  u"  -  [  .  7'  7.  t  rn  fij  .7.  7.  1-.-,  n,  t  c.j 
c/:  ,-,j  ,-1  j  ,‘ij  CO  — •  y  aa  I  -  t  ."'J  iV|  T  "1  "Vi  /.j  £i  .'ij  O'1  J.  —1 

r.-i  iv,  ,  OJ  'v'  '•'<  O'i  '"ij 


,  J-,  .  J~i  iT,  IT1  f  _  I  Vi  *£»  i'»J  -Vi  £|  “■  f«.  r  -  »*0  OJ  «— 1  f-  — ‘  T— I  .V,  U”1  ‘*0 

•=•  Aj  ^  co  •£•  <- ’t  Aj  y  aa  =•  t  t  y  j-  £•  t  y  y  >7.  u*>  ~  oj 

21  ,-|J  f  -  u'i  aj  -|J  AA  rt  OJ  AJ  AJ  ,  A  IJ-I  rf  .Aj  f  „  ,T.  ,-|J 

.V,  ,v,  AJ  C  A  AA  Aj  .'A  .  A  AA  AA  AA  AA  AA  ^  .AJ 


CO  r  -  ~  -  '£' 


O  CO  -•  CO  t  - 


Li  AJ  ijA  r-._  _■  f..  -■ 


'  £1  AA  _i  'L'  ‘L'i  O'l  O  '  I)"'  '  0  CO  I  -  AA  CO  L'  U'i  f  -  <7  •  U~  ~T  U"'  ~ 

CO  CO  CO  :  '  0  '?•  OJ  '7  ■  CO  CO  CO  CO  O'l  .7  CO  CO  '7 1  '7  ■  '7  ■  — •  T  cr. 

,-»j j  (-,j  i  oj  co  oj  ivi  oj  co  oj  *’'j  ru  oj  r,j  oj  oj  oj  oj  oj  cu  co  ro  oj  o j  ro 


f...  f  ._  yj  CO  CO  CO  '  L‘  U'.'  AA  — ■  AA  CO  -H  f...  S  If  AA  f  ..  f  -  It  —I  .Ti  CO 

ro  iT  .  o  CO  AJ  ^  O  IT.  f .  *T  f...  rj-  CO  O  f  -  £1  U'I  CO  •■£'  U'I  U'I  CO  O  AJ  U'i 

AA  AJ  L'  CO  CO  U'I  |7-.  AJ  >7'  CO  CO  CO  *7'.  L'  0  If'  f  -  •£>  rf  <3>  -<  i-*  CO  rT 

AA  Aj  AA:  CO  CO  CO  CO  C'J  CO  CO  CO  CO  CO  CO  0 j  CO  AJ  CO  AJ  AJ  ^  CO  CO  CO  Co 


Ij-f  f...  ,  A.  T  CO  OJ  OJ  f:  rf  y  CO  ‘  lit  AA  rf  y  rf  |f)  UO  CO  CO  |T.  AJ  £i 

.J*  CO  CO  Tf  AJ  AJ  *.£l  f  -  .7:  ..£■  iT,  :.£|  AA  -<  i-  O  f--  ’t  U*'  CO  ''0  £l  H  £i 

y  *f  aa!  7  y  co  oj  uo  aj  y  u  5  y  y  y  y  t  y  0  y  co 

AA  AA  aa  CO  CO  CO  CO  CO  CO  CO  CO  ‘O  CO  AA  ro 


co  co  li"'!  ■£' 


C'J  CO  »r  CO  f-  Aj  rf  y  if.  aa  .t.  o  CO  U'  J.  aa  f-._  rf 


yr  rf  .“'|  OJ  i“i  'L'  f-  CC-  '.‘A  rf  ,T.  AA  T*.  rf  ij'.  .£i  *7 .  ’t  t  i.*A  |7.  £i  j“t 
O  O  '7'.  AA  AJ  rf  C  Ij'l  -  -■  -r  AJ  ,AJ  ,AJ  AA  O  AA  y  j"'  CO  AA 

AA  A-,  rri,  C'J  C'A  CO  CO  CO  aa  CO  CO  CO  CO  CO  aa  aa  aa  ro  CO  CO  CO  CO  AJ  CO  AA 


i,  c* .« ro  ». '  i  oj  o  >  *. 1  • 


i',j  ^  ,J-,  .£|  -Ti  r-_  iT.  iT.  u*  «7*  OJ  L*  — ’  ’t'  Vi  CO  T.  .Vi  •  rf  «T.  £1  rf  »V» 

•  •■I  •■"••**. *"*•"**  ■■■••••• 

^  iT.  co  'J*1  •ro  '■£»  u-'  rf  r_  ..o  r  -  u  •  l-  j  «.u  oj  r.  rf  it,  ^  r  _  u  -o  r  -  'ij  *  *.» 

^  o  Oj  -7-  oj  oj  .  o  rf  -r-  -r-  •  v*»j  *ij  'ij  ^  oj  ro  *-»  J*  ^  oj 

,vi  i'i'i  o"r  o j  ro  r*’*  ro  ro  ro  ro  ro  ro  ro  ro  ro  ro  ro  ro  ro  ro  ro  ru  ro  ,V| 


U’i  ** 


0  CC*  *-*  |7‘*  r--  rf  7  • 


1.1  ..IJ  U  •-- 


v  ’.i.'  r*  Li  * 


I  r— <  I  •  ■  .  __ 

lol  <7-i  r--  £'  CA  f  -  C'A  AA  ->  -•  (■-  U*I  u'i  7 •  r  -  L'  5‘  U'I  .7.  -•  O  CC  •£'  AJ  S' 

)  |  iT.  iT*  7.  Oj  w-t  f—i  i'iJ  ‘H1  ‘'*'i  iT.  7.  '7'  7  ’  — '  — '  •—  —  >7'  ^  i“ij  'jj  .■  i-h 

oj  oj  ru  oj  *0  ro  ro  ro  ro  ru  ru  ru  oj  ro  *.*,  ,v,  (*o  oj  ro  o-»  ro  ru  ro  -o 


*  .  . 

»7*  OJ  rf  ro  ?■•-  CO  -7*  u“  0  »-*  ^  U“'  U~  0  *-*  '7  -  f  -  *T  -.U  u  *  n- 

£i  y“i  ijj  ru  co  7-  — ■  *7  ■  ru  >7*  co  co  r-  u_i  »-*  *’o  n  ^  <7* 

ru  7- »-»  jo  go  c-HfH  i-ij  r-  rf  f  _  »r « 'i-  r  -  >0  "7  f-  Lf'  n 

•-4  j'i  u“i  ijO  r  .  “  •  7-4  OJ  CO  ro  ».*0  rf  rf  rf  £i  r*-  «7-  — '  — '  — ' 

I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I.  • 

CA  <r  aa  aa  aa  aa  aa  aa  i  t  7  7  7  7  7  t  7  t  7  t  7  'f  T  l' 

<n  cm  .  -  r-.  r  -  f  -  r  -  r  -  r  -  r-.  r.  r  -  t  -  f  -  r-  r  _-  t-  r-r  r  -  *;•-  r:  r: 

O-t  D«  iji  ,Ti  ,v»  'l’i  .*i*i  Yi  .  .Vi  i'i*'  O';  .‘l’i  i*i‘i  I'M  CO  CO  CO  CC1  •  0  '-0  •-,.t 


Refer  to  Table  1 1  —  1  for  T  designations  used  in  Figure  II- 


Station 

Analysis  Record  Length,  days 

RMS  Error 

873-1269* 

182 

0.06 

873-1952 

182 

0.13 

873-5184 

182 

0.09 

873-5523 

141 

0.12 

873-5587 

182 

0.10 

873-7048 

182 

0.15 

874-0199 

182 

0.10 

874-0405 

182 

0.08 

874-1196 

182 

0.10 

874-2221 

182 

0.09 

874-3081 

182 

0.16 

874-3495 

182 

0.17 

874-3735 

182 

0.16 

874-4586 

182 

0.16 

874-4671 

182 

0.18 

874-4756 

122 

0.10 

874-6819 

182 

0.16 

874-7437 

182 

0.17 

874-9704 

182 

0.17 

876-0412 

182 

0.07 

876-0595 

146 

0.06 

876-0742 

182 

0.11 

P22 

182 

0.05 

P23 

182 

0.05 

P24 

159 

0.09 

*  Refer  to  Table  1 1  —  1  for  T  designations  used  in  Figure  1 1 - 3 . 


17.  The  harmonic  analysis  included: 

a.  Editing  of  the  current  data  and  selecting  periods  which  con¬ 
tained  the  longest  continuous  data  return. 

b.  Filtering  to  remove  high  and  low  frequency  trends. 

c.  Determination  of  the  tidal  constituents  for  both  current 
components . 

18.  The  same  harmonic  analysis  procedures  as  outlined  previously  for 
water  surface  elevations  were  employed  in  the  current  component  analysis. 
Velocity  component  means  are  given  in  Table  III-5.  Amplitude  and  phases  for 
each  component  are  given  in  Tables  1 1 1 —6  through  1 1 1 - 9 .  In  these  tables,  S,  M, 
and  B  denote  surface,  middepth,  and  bottom  meter  locations,  respectively.  At 
stations  where  only  one  meter  is  located  it  is  denoted  by  M. 


Table  I I I - 5 

Station  Mean  Velocity  Components 


Station 

1980  Analysis  Period 

Mean , 

fps 

Start 

End 

North/South 

East/West 

Vl-M* 

9 

Apr 

17 

Jun 

0. 

01 

N 

0. 

16 

E 

V2-M* 

9 

Apr 

17 

Jun 

0. 

01 

S 

0. 

1 1 

E 

V3-S* 

20 

Apr 

17 

Jun 

0. 

00 

0. 

02 

E 

V3-B 

20 

Apr 

17 

Jun 

0. 

05 

N 

0. 

03 

E 

V4-S* 

9 

Apr 

17 

Jun 

0. 

02 

N 

0. 

02 

W 

V4-B 

9 

Apr 

17 

Jun 

0. 

08 

N 

0. 

00 

V5-M 

20 

Apr 

12 

Jun 

0. 

04 

N 

0. 

05 

E 

V6-S* ,** 

10 

Apr 

17 

Jun 

0. 

01 

N 

0. 

12 

E 

V6-B** 

10 

Apr 

14 

Jun 

0. 

04 

N 

0. 

05 

E 

V7-S 

10 

Apr 

17 

Jun 

0. 

02 

N 

0. 

12 

E 

V7-M* 

20 

Sep 

21 

Oct 

0. 

01 

S 

0. 

02 

E 

V7-B 

2 

May 

17 

Jun 

0. 

04 

N 

0. 

01 

E 

V8-S 

10 

Apr 

17 

Jun 

0. 

08 

N 

0. 

13 

W 

V8-M* 

10 

Apr 

17 

Jun 

0. 

08 

N 

0. 

18 

w 

V8-B 

2 

May 

30 

May 

0. 

05 

N 

0. 

03 

w 

V9-S* 

11 

Apr 

19 

Jun 

0. 

14 

S 

0. 

06 

E 

V9-B 

2 

May 

19 

Jun 

0. 

00 

0. 

00 

V10-S*,** 

10 

Jul 

21 

Oct 

0. 

02 

S 

0. 

01 

E 

V10-B** 

10 

Jul 

21 

Oct 

0. 

.  14 

N 

0. 

04 

E 

VI 1 -M* 

3 

May 

15 

Jun 

0. 

.03 

S 

0. 

.02 

E 

V12-S* ,** 

5 

May 

19 

Jun 

0, 

.01 

S 

0. 

.  10 

E 

V 1 2 - B** 

16 

Apr 

19 

Jun 

0. 

.01 

N 

0. 

.  10 

E 

V13-S* 

11 

Apr 

20 

Jun 

0 

.01 

S 

0. 

.05 

E 

V13-B 

11 

Apr 

20 

Jun 

0 

.05 

N 

0. 

.03 

E 

V14-M* 

11 

Apr 

20 

Jun 

0 

.07 

S 

0, 

.11 

E 

V 1 5  -  S'-'"’' 

3 

May 

20 

Jun 

0 

.05 

N 

0, 

.03 

W 

V15-M*,** 

5 

May 

20 

Jun 

0 

.  12 

N 

0 

.07 

W 

V 1 5 - B  ** 

16 

Apr 

20 

Jun 

0 

.  13 

N 

0 

.11 

W 

V16-S 

16 

Apr 

7 

May 

~  "" 

Vlb-B 

16 

Apr 

1 

May 

-- 

0 

V17-M* 

12 

Aug 

22 

Oct 

0 

.02 

N 

.  00 

V18-M* ,** 

12 

Aug 

22 

Oct 

0 

.02 

N 

0 

.07 

W 

V19-S* 

16 

Apr 

20 

Jun 

0 

.06 

S 

0 

.06 

E 

V19-B 

16 

Apr 

20 

Jun 

0 

.03 

S 

0 

.05 

E 

V20-M* ,** 

17 

Apr 

20 

Jun 

0 

O 

OO 

N 

0 

.05 

W 

V21-S 

22 

Apr 

20 

Jun 

0 

.51 

S 

0 

.33 

W 

V2 1 -M* 

22 

Apr 

20 

Jun 

0 

.34 

s 

0 

.  19 

W 

V21-B 

22 

Apr 

15 

Jun 

0 

.29 

s 

0 

.10 

W 

Stations  considered  in  the  two- 
shown  in  Figure  II-l. 

dimensional  numerical  modeling  effort 

as 

These  stations  are  designated  VS  in  Figure  II-l. 
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Table  Iir-6 

North/South  Velocity  Component  Mean  Amplitude  (f p s ) 


Station 


Vl-M 

0.01 

02-M 

0.04 

03-S 

0.20 

03-B 

0.25 

04-S 

0.33 

V4-B 

0.19 

V5-M 

0.08 

V6-S 

0.04 

06-B 

0.06 

07-S 

0.08 

07-M 

0.08 

07-B 

0.07 

08-S 

0.29 

08-M 

0.17 

08-B 

0.08 

V9-S 

0.55 

V9-B 

0.38 

010-S 

0.23 

010-B 

0.22 

Yll-h 

0.09 

012-S 

0.02 

V12-B 

0.07 

V13-S 

0.11 

V13-B 

0.13 

014-M 

0.51 

V15-S 

0.34 

V15-M 

0.26 

015-B 

0.22 

0 1 7 -M  < 

0.22 

018-M 

0.56 

019-S 

0.08 

V19-B 

0.06 

020-M 

0.07 

021-S 

0.55 

V21-M 

0.34 

021-B 

0.26 

_ Constituent 

Ml  J1  01 


0.02  0.04 

0.01  0.01 

0.02  0.03 

0.02  0.02 

0.03  0.05 

0.01  0.01 

0.01  0,03 

0.01  0.02 

0.02  0.02 

0.01  0.01 

0.01  0.03 

0.03  0.05 

0.02  0.05 

0.02  0.01 

0.03  0.03 

0.03  0.05 

•  03 


02 
03 
0.04 
05 
03 
02 
03 


01 

0.01 

01 

0.01 

06 

0.03 

04 

0.05 

05 

0.05 

06 

0.03 

02 

0.04 

01 

0.02 

03 

0.03 

02 

0.02 

02 

0.01 

02 

0.02 

12 

0.08 

08 

0.05 

02 

0.03 

21 

0.15 

1 5 

0.06 

08 

0.04 

07 

0.05 

04 

0.03 

02 

0.02 

02 

0.03 

02 

0.05 

03 

0.06 

16 

0.10 

12 

0.07 

07 

0.06 

06 

0,04 

09 

0.05 

31 

0.14 

04 

0.03 

02 

0.02 

03 

0.02 

19 

0.07 

9  r> 

0.06 

21 

0.08 

Table  III-7 

North/South  Velocity  Component  Local  Epoch  (deg 


Table  III-8 

East /West  Velocity  Component  Mean  Amplitude  (f_ps) 


v  1  -  m 

0 . 33 

0.4 

,  v  2  --  h 

0 . 36 

0.4 

v  3  -  S 

0. 48 

0.6 

B  03 -B 

0.41 

0.5 

i  o  4  -  s 

0 . 4  7 

0 , 6 

V  4  -  b 

0.33 

0 . 4 

■  ■  Z  -  H 

0 . 24 

0.3 

V  6i  ~  '~i 

0.27 

0 . 3 

'.'6  -  B 

0.19 

0 . 2 

v  .  <  *  3 

0 . 28 

0 . 3 

I  -'7-M 

0.21 

0 . 1 

\  '  -  Zi 

0,20 

0.2 

.3  S 

0.42 

0.4 

V8-h 

0 . 40 

0 . 4 

:  oa-B 

0 . 20 

0 . 2 

V  9  -  S 

0,07 

0 . 0 

L  09-B 

0.06 

0 . 0 

*  V 1  0  -  s 

0 . 23 

0 . 1 

V10-B 

0 . 28 

0 . 2 

•m-H 

0 . 08 

0.0 

V12-S 

0.12 

0. 1 

V12-B 

0.08 

0. 1 

VI3-S 

0 . 08 

0. 1 

|  y  i  3  -  s 

0.03 

0.0 

l  4  -  h 

0.11 

0.1 

vi  r»  -  s 

0.06 

0.0 

0  1  5  -  M 

0  .  OS 

0.0 

V  1  5  -  B 

0.0  4 

0.0 

; .  017-rt 

0 .04 

0.0 

r  018-M 

0.28 

0.2 

*  019-S 

0.16 

0 . 1 

■  V19-B 

0.13 

0.1 

r  02-0-M 

0.29 

0.3 

[  V21-S 

0.30 
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Table  III-9 


East /West  Velocity  Component  Local  Epoch  (deg) 


Vl~h 

136.1 

124.4 

263.8 

350 . 8 

101.3 

344 . 9 

240.4 

281.5 

02 -M 

128.7 

126.7 

230.6 

16.0 

108.8 

340.5 

'  238 . 6 

275 . 4 

V3-S 

129,4 

126.0 

189.0 

84.2 

105.8 

316.9 

n  n  c"  q 

i-  ♦  / 

289 . 5 

V3-B 

126.1 

126 . 0 

196.7 

54 . 3 

99 . 4 

322 . 2 

231.5 

279 . 7 

04-S 

100.1 

100.7 

212.5 

301 . 6 

65 . 4 

286 . 2 

186.3 

237 . 7 

04- B 

106.0 

91 . 7 

119.4 

220.5 

50.3 

284 . 8 

198.5 

178.8 

05-M 

97.3 

95.1 

350.7 

274.2 

56 . 4 

282.6 

178.2 

198.5 

06-S 

109.1 

104.1 

o  'l  y  t 

19.8 

63 . 0 

291.0 

216.5 

278 . 9 

06-B 

105.1 

103.4 

240.9 

302.7 

1 . 0 

289 . 6 

228.9 

286 . 2 

07- S 

96 . 1 

104.8 

218.5 

14.3 

89 . 1 

284 . 1 

199.3 

53 . 0 

07-M 

93.8 

107.2 

154 . 9 

249 . 7 

53  .  1 

271.7 

177.4 

141.6 

07- B 

95.6 

95 . 9 

155 . 7 

120.2 

93 . 9 

272 . 1 

199.1 

209 . 8 

08-S 

38 . 9 

50 . 2 

342 . 4 

251 . 3 

28 . 6 

217.9 

153.4 

161.7 

08-M 

43  .  1 

51.4 

181.8 

322.1 

50 . 5 

212.2 

152.0 

12  5.9 

08-B 

28 . 4 

41.4 

296.5 

1  . 1 

349 . 6 

195.7 

14  2.6 

79 . 2 

09-S 

56 . 8 

45 . 8 

15.0 

172.3 

344 . 7 

2  2  4.2 

14  1.5 

80.4 

09-B 

53.0 

68 . 1 

25 . 3 

113.0 

36 . 0 

235.6 

144.6 

155.6 

oio-s 

240 . 2 

256 . 6 

260 . 0 

348.0 

194.6 

33.6 

292 . 3 

331.1 

oio-b 

235.5 

263 . 4 

279 . 6 

115.8 

221.5 

42.5 

331.4 

316.3 

0 1  1  -  M 

4 . 6 

3  2  6.2 

91.1 

357. 7 

65 . 0 

295 . 9 

39.0 

255 . 8 

012-S 

330 . 3 

347.1 

58 . 1 

192.1 

347 . 7 

204 . 7 

94.0 

136.1 

0  1  2-B 

340 . 8 

332.1 

44.2 

118.7 

243.6 

145.4 

90 . 0 

30 . 8 

013-S 

331.3 

312.2 

49 . 6 

174.4 

253.9 

219,7 

65  .  V 

141.1 

013-B 

328.9 

331.7 

264.0 

171.4 

314.8 

170.7 

91 .8 

106.2 

014-rl 

246 . 8 

2  4  5.3 

27.5 

231 .4 

241 .3 

78. 2 

324  .  B 

318.7 

015-S 

319.8 

2  9  7. 5 

8 . 0 

81.1 

354 . 3 

275.2 

245 . 6 

263 . 0 

015-rt 

285 . 2 

292 .0 

29^,9 

53 . 3 

309 . 2 

113.4 

313.6 

279 . 1 

V 1  5  -  B 

3  3  1.3 

38 . 6 

115.2 

203.7 

3  2  0.2 

132.2 

23 . 5 

291.2 

V17-M 

4  .  4 

3  06 . 0 

6  7.5 

328 . 4 

76 . 4 

151.0 

95 . 4 

64 . 1 

V18-M 

28.7 

30 . 3 

120.4 

356 . 0 

n  9  .  n 

169.5 

81.1 

91.8 

019-S 

2  7  7. 7 

2  6  6.2 

297.0 

343 . 5 

262 . 1 

20.0 

340.7 

341.2 

V  1  9  -  B 

2  66.1 

25  b .  / 

253 . 6 

340.2 

243.9 

10.3 

351.3 

3 . 6 

020-M 

2  8  6 . 8 

285 . 7 

260 . 4 

307.0 

243.4 

94 . 0 

41.9 

23 . 0 

V  2 1  -  S 

2  4  1.1 

24  7. 6 

23 . 7 

162.0 

185.4 

36 . 0 

323 . 6 

305 . 2 

V21-M 

241.7 

228 . 5 

330 . 0 

6.5 

157.0 

45.7 

47.0 

317.3 

021-B 

240.5 

221.3 

269.5 

41.2 

182.0 

31 .5 

45 . 4 

8.7 

PART  IV:  GULF  TIDE  MODEL  DEVELOPMENT  AND  RESULTS 


19.  Consider  the  depth  integrated  linearized  equations  of  motion  and 
continuity  for  a  homogeneous  and  incompressible  fluid  in  spherical  coordinates. 
Reid  and  Whitaker  (1981)  have  employed  alternating  direction  implicit  finite 
difference  approximations  to  these  equations.  In  an  application  to  the  Gulf 

of  Mexico  a  15-  by  15-minute  latitude  and  longitude  grid  employing  2228  compu¬ 
tational  cells  has  been  employed.  Wind  forcing  and  stratification  effects 
have  not  been  considered.  The  tides  are  lot  red  by  direct  tide  potential  and 
by  volume  transport  (flux)  at  the  Florida  Strait  and  Yucatan  Channel. 

20.  Tidal  forcing  (flux  potential  at  the  two  ports  and  the  effective 

amplitude  and  phase  of  the  direct  tide  potential)  is  determined  separately 
for  each  of  the  five  major  tidal  constituents  in  the  Gulf  (0,  ,  K,  ,  P,  .  M.,  , 
and  S^)  by  minimizing  the  sum  of  the  squared  errors  between  the  measured  con¬ 
stituent  amplitude  and  phase  and  the  model  constituent  amplitude  and  phase  for 
all  gaging  stations.  In  the  original  application,  Reid  and  Whitaker  considered 
20  stations  located  around  the  entire  Gulf  of  Mexico.  In  the  least  squares 
analysis,  each  of  the  20  stations  was  given  equal  weight.  In  order  to  improve 
results  along  the  northern  Gulf  coast  stretching  from  Atchafalaya  Bay  to  Pensa¬ 
cola  Bay,  stations  along  this  coastal  segment  were  weighted  more  heavily  than 
those  occurring  outside  this  segment.  The  least  squares  analysis  was  repeated 
for  all  five  constituents  with  the  final  results  for  the  amplitudes  (centi¬ 
metres)  and  phases  (Greenwich  epoch)  presented  in  Tables  IV-1  through  IV-5  for 
the  0,  ,  K,  ,  P,  ,  M„  ,  and  S„  tidal  constituents,  respectively. 
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PART  V:  SALINITY  ALGORITHM  DEVELOPMENT 


21.  Salinity  is  considered  as  a  conservative  (passive)  constituent. 

The  three-dimensional  passive  constituent  transport  equation  for  laminar  flow 
is  presented  followed  by  the  modifications  necessary  for  turbulent  flow.  The 
turbulent  flow  equation  is  depth  integrated.  The  transport  equation  obtained 
is  then  transformed  using  an  exponential  stretch.  Numerical  approximations  to 
the  transformed  equation  are  presented  followed  by  the  development  of  rela¬ 
tions  for  the  effective  dispersion  coefficients. 


Cons  L  i  t lie n t  T r a nsport  Equation  i n  Ca rtesian  Coordinates 


22.  The  constituent  transport  equation  is  given  for  laminar  flow  as 


(V.l) 


where 

s  =  concentration  of  the  material  of  concern 

D  -=  molecular  diffusion  coefficient  in  the  x  direction 
x 

D  5  molecular  diffusion  coefficient  in  the  y  direction 

y  * 

D  ~  molecular  diffusion  coefficient  in  the  z  direction 
z 

x,y,z  =  Cartesian  coordinates 

u,v,w  =  velocity  components  in  the  x,  y,  and  z  directions,  respectively 
t  =  time 


For  a  turbulent  flow,  the  eddy  dispersion  is  significantly  larger  than  the 
molecular  diffusion.  The  following  analogous  formula  holds  where  time  averag¬ 
ing  over  the  time  scale  of  the  turbulence  has  been  performed. 


(V.  2) 


where  ,  Ky  ,  and  are  turbulent  eddy  dispersion  coefficients.  Equa¬ 

tion  V.2  may  be  written  in  conservation  form  by  adding  s  times  the  continuity 
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equation  (namely,  zero)  to  the  left-hand  side  to  obtain 


9s  +  9(us )  +  9 (vs )  9(ws  ~ 

9t  9x  9y  9z 


a_ 

U 

.  3_  { 

K  9s  \ 

+  2L  1 

fK  *±\ 

3x 

\  x  9x/ 

9y  ' 

V  y  dy) 

dz  1 

*z  9z / 

(V.3) 


This  form  of  the  equation  is  then  depth  integrated  as  described  in  Schmalz 
(1983a)  to  obtain: 


9  /■  i_  \  9  ,,  v  9  ,  .  9 

9t  (h8)  +  9x  (hus)  +  9y  (hV5)  =  9x 


(“I  © +  ry  H  if)  <v  4> 


where  h  is  the  water  depth  and  K*  and  K*  are  effective  dispersion 

x  y 

coef f icients . 


Constituent  Transport  Equation  in  Transformed  Coord inates 

23.  The  transport  equation  is  transformed  from  x-y  space  to  a\~a2 
space  by  means  of  the  following  coordinate  transformation  as  considered  by 
Butler  (1980). 


X  =  3!  +  Vl 


(V.5) 


y  =  a2  +  b2a2 


y  -  a. 


(V.6) 


Then  for  an  arbitrary  hydrodynamic  variable  p(x,y,t) 


9£)  =  9£_^1  9p  _  9p  d°2 

9x  9a  dx  9y  9a^  dy 


(V.  7) 


I f  we  introduce  p^  =  dx/da^  and  p^  =  dy/da^  then 


9£_l_9g_  9p  _  1  9p 

9x  “  p]  90j  9y  P2  3«2 


(V.  8) 


24.  Transforming  Equation  V.4  in  x-y  space  to  a^-a.. 
the  following  result. 


space  we  obtain 


+  - 
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( V .  9  I 


(ds)  + 


1  M, 


M2  Mj 
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°2  M2 


,  «2 


where  d  i.s  used  to  indicate  water  deptli  in  place  oi  h  and 
(  )t  =  3/9t 


(  ^  =  a/acrj 
(  \2  =  a/a°2 


Equation  V.9  is  the  relation  that  is  the  subject  of  numerical  approx  imal  i  on . 


Numerical  Approximations 


25.  Schmalz  (1983a,  1983b,  1983c)  considered  several  alternate  tech¬ 
niques  for  approximating  the  linear  form  of  Equation  V.9.  The  Flux-Corrected 
Transport  (FCT)  scheme  was  selected  as  the  most  accurate  scheme  and  has  been 
incorporated  in  Lhe  Waterways  Experiment  Station  Implicit  Flooding  Model  (WIFM.l 
In  addition  a  Three  Time  Level  Explicit  Transport  scheme  was  also  incorporated 
in  the  model.  A  space  staggered  grid  as  shown  in  Figure  V.l  was  employed  in 
all  of  the  formulations.  The  datum  convention  is  presented  in  Figure  V.2. 

26.  Let  us  introduce  the  following  notation  as  a  prelude  to  the  approxi- 

mations.  Define  for  an  arbitrary  variable  F  ,  where  t  =  kAt  ,  y  -  nAy  , 

n ,  m 

x  -  mAx  : 
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(V. 10a) 


(V. 10b) 


(V. 10c) 
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Flux- co r rected  transpor t 

27.  Two  schemes  are  used  in  implementing  this  approach:  a  lower  order 
in  space  nonosci 1 latory  scheme  and  a  higher  order  in  space  scheme  subject  to 
oscillation.  In  the  method  implemented,  two  time  level  implicit  multiopera- 
tional  ADI  schemes  were  employed.  The  forward  time  upwind  space  (FTUS)  and 
forward  time  centered  space  (FTCS)  schemes  were  used  as  the  lower  and  higher 
order  in  space  schemes,  respectively,  and  are  discussed  in  turn  below.  Fi¬ 
nally,  the  necessary  flux  correction  procedures  are  developed. 

Leendertse  FTCS  mul t iope rati ona 1  scheme 

28.  The  following  finite  difference  equation  is  considered  as  an  ap¬ 
proximation  to  the  nonlinear  transport  Equation  V.9. 
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29.  The  above  equation  is  assumed  to  be  contained  within  the  alternat¬ 
ing  direction  difference  equations  presented  in  paragraphs  .10  and  31  below. 

For  the  linear  case  obtained  for  (lk)  =  (li.)  =  I  ,  K  ,  K  constants 

2  n  1  m  a  a 

Z  1 

in  space  and  time,  and  u  ,  v  ,  d  constant  in  space  and  time,  the  consti¬ 
tuent  intermediate  time  level  (k+l/2'-v  in  paragraphs  30  and  31)  may  be  elimi¬ 
nated  in  the  alternating  direction  approach  and  the  total  difference  equation 
obtained  equals  the  above  difference  equation  plus  some  higher  order  in  time 
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factorization  terms.  The  total  difference  equation  is  consistent  with  the 
linear  transport  equation.  For  the  nonlinear  case  considered,  it  is  not  pos¬ 
sible  to  eliminate  the  constituent  intermediate'  time  level.  Thus  the  exai  t 
form  of  the  factorization  terms  can  not  be  determined.  However,  their  numeri¬ 
cal  effect  can  he  tested. 

30.  The  approximations  for  the  X- Sweep  m...y  now  he  written  as  follows: 


II  we  place  all  terms  at  time  I  •  ■  \ »  I  K+l/2'  (,n  it,,-  l.-fl-hand  side  if  the 
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(V.  13) 
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Collecting  all  terms  in  Equation  V.12  at  time  level  k  ,  denoting  the  result 

as  B  ,  we  obtain  with  K  =  K 
m  y  a 
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(V. 17) 


11.  Collecting  all  results  we  obtain  the  following  interior  equation 
for  the  X-Sweep 
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32.  The  approximations  for  the  Y-Sweep  may  now  he  written  as  follows: 
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Expanding  Equation  V.19  by  employing  V.10  and  collecting  terms  at  time  level 
k+1  on  the  left  side  arid  leaving  terms  at  time  level  K+l/2'-'  on  the  right 
side  the  following  interior  equation  for  the  Y-Sweep  is  obtained. 


k+1  k+1  k  + 1 

a  s  +  a  s  +  a  s  .  ,  -  r 

n-  ,:u  ii-!, m  n,m  n,m  n+l,m  n+l,m 


if  re  K  -  K  and  K  =  K 

x  Uj 


.  -  k+1 

d  n- 1  /  2 , in 
,)n 


vk+1  (R  lk+1 

n-1/2 ,m  y'n-l /2 ,ro 


2(u2)n-l/2 J 


( V  .  2  0  > 


( V .  2 1 ) 


an+l  ,n 


v 


k+1 

n+1/2 


2 


,  m 


<K  iktk7 

y  n-f 1 ,  2  ,ro 


n+ 1  /  2 


( V .  22  ) 


+ 


.it 


k+1 

n  +  1  /2 ,  m 
n+1  /  2 


k+1 

n-1/2 


( u ) 


n-1/2 


m 


(V.2  i| 


'n.ra  2(„ 


it  — )  k+1/2*  k+l/2*  )  K+l/2*  k+l/2*l 

IcA^iij  L'  d*  '  n.Brtl/2  Un,o+-l/2  '  *  /  n.B-1/2  un,B-l/2j 


.  \(\Vu,  <£*%jJ2D . (I b-  <•-•" •• £#| 

1  m  I 


( V . 24 ) 


Leendertse  FTUS  multioperational  scheme 

33.  The  following  finite  difference  equation  is  considered  as  an  approx¬ 
imation  to  the  nonlinear  transport  Equation  V.9: 
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=  0  at  (n,m) 


34.  The  following  upwind  difference  operators  used  in  the  above  equa¬ 
tion  are  defined  at  (n,m)  as  follows: 
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35.  For  the  linear  case  Mu)  =  (p„)  =  1.0  ,  K  ,  K  ,  u  . 

|_  1  m  2  n  °]  a2 

and  d  constant  ,  the  constituent  intermediate  time  level  in  the  alternating 
direction  approach  may  be  eliminated.  The  difference  equation  obtained  is 
consistent  with  the  linear  transport  equation  and  equals  the  above  difference 
equations  plus  some  higher  order  in  time  factorization  terms.  For  the  non¬ 
linear  case  considered  here,  it  is  not  possible  to  eliminate  the  constituent 
intermediate  time  level.  Therefore  the  exact  form  of  the  factorization  terms 
can  not  be  determined.  However,  their  numerical  effect  can  be  assessed. 

36.  This  scheme  is  similar  to  the  standard  ADI  technique  except  that  up¬ 
wind  differencing  is  employed  for  the  advective  terms.  The  necessary  modifica¬ 
tions  for  the  X-Sweep  are  shown  in  Table  V- I  and  those  employed  for  the  Y-Sweep 
are  given  in  Table  V-2. 
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Table  V-2 

Y-Sweep  Modifications  FTUS 
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Flux  correction  procedures 

37.  If  the  factorization  terms  are  ignored,  the  schemes  above  may  be 
written  in  the  following  flux  format. 


,k+l  I  k  ck 

d  s  =  d  S 
n,m  n,m  n,m  n,m 


[toi(M,)./l«2(M2)„] 
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m+1/2  n,m-l/2, 
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where  t  =  kAt  ,  x  =  £  ((Jj)  .Ao^  ,  y  =  £ 

i  i 

S  =  concentration  at  location  (n,m)  at  time  level  k 
n  ,m 

AOj(|J  )  =  x  space  step  at  m 
^02^2^m  E  y  space  step  at  n 


I  =  general  index  at  time  level  k+1  ,  which  we  set  to 
H  or  L  for  the  higher  or  lower  scheme, 
respectively 


Fn±l/2,m±l/2 


fluxes  through  the  appropriate  cell  faces  of  cell 
(n,m).  Form  is  dependent  upon  the  finite  difference 
formulation 


We  observe  from  Equation  V.27  that  the  difference  between  the  higher  and  lower 
order  scheme  at  (n,m)  may  be  written  as  follows: 


(C  -  S»,J-  -K'V.^'V.d  1  [(fn+l/:,n  -  Fii+ 1  /  2 ,  m  ) 


(Fn-l/2,m  "  Fn-1 /2 ,m)  *  ('  n. 


m+1/2  n,r+l 


( V . 28 ) 


(^n , n- 1 / 2  ' n , m- 1  ’ 


Note  that  this  difference  can  be  expressed  as  an  array  of  fluxes  between  adja¬ 
cent  grid  points  and  is  the  condition  required  for  flux  totreition.  We  next 

H  I 

develop  the  flux  expressions  for  the  higher  (F  )  and  lower  (FI  order  schemes. 
In  order  to  aid  in  notation,  we  make  the  following  definition  for  an  arbitrary 
variable,  F  . 
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38.  For  the  higher  order  scheme  we  employ  the  FTCS  scheme  wiitten  in 
Equation  V.ll  in  which  the  factorization  terms  developed  in  the  multiopera- 
tionai  method  are  not  shown.  Equation  V.ll  may  he  written  in  the  form  of  V.27, 
where  the  total  fluxes  are  presented  as  the  sum  of  advert ive  and  diffusive 


f 1 uxes . 


39.  From  Equation  V.ll  one  then  obtains  for  the  advert ive  fluxes: 


F  A  =  k+1/2  ,  , 

ln+l/2,m  '"n+l/2,m“  m"  *1 


SH  4  S; 


i  dk:;/j 

n+]fm 
n+1 ,  m  — 


+  ^k+l  /2  A  (v.30) 


_H A  _  k+1/2  .  / SH  +  Sk 

Fn  ,m+l/2  Un,m+l/2L  tU2)nAci2  \  2 


n  ,nr+l 


k+1/2 
n  ,m+l 


♦  c'72' 

v  n ,  n  J  / 


(V. 31) 


Note  the  subscript  A  denotes  advection.  The  diffusive  fluxes  are  then  given 


by  the  following  relations. 
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Note  the  subscript  0  denotes  diffusion. 

40.  For  the  lower  order  scheme,  the  FTUS  scheme  written  in  Equation  V.25 
is  employed.  Factorization  terms  generated  by  the  mu ] t iopera t iona 1  method  ;  re 
not  considered.  Equation  V.25  is  written  in  the  form  of  V.27.  The  total 
fluxes  are  presented  as  the  sum  of  advective  and  diffusive  fluxes. 

41.  From  Equation  V.25  one  obtains  the  following  set  of  advective 
f 1 uxes . 
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The  diffusive  fluxes  are  obtained  from  Equations  V.32  and  V.33  with  H  re¬ 
placed  by  L  . 

42.  The  ant  id i f f us i ve  fluxes  are  then  computed  as  follows. 
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(V. 39) 


In  computing  the  difference  between  the  diffusive  fluxes  (third  and  fourth 

terms  in  the  above  express i ons ) ,  note  that  the  terms  with  S  are  of  opno- 

n  ,  m 

site  signs  and  drop  out  of  the  equations. 

43.  Next  the  maximum  and  minimum  cell  values  are  determined. 
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44.  Next  the  sum  of  all  ant idi f fusive  fluxes  into  cell  (n,m), 


is  determined. 


P  =  max  (0,A  )  -  min  (0,A  ) 

n  ,m  n-1/2  ,m  ’  n+l/2,m 


+  max  (0,A  .  )  -  min  (0,A  .  )  (V.43) 
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The  maximum  allowable  mass  into  the  cell,  Q+  ,  is  then  computed  as  follows: 
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45.  Similarly,  the  sum  of  all  ant i di f f us i ve  fluxes  into  cell  (n,m), 
,  is  determined. 


P  =  max  (0,A  .„  )  -  min  (0,A  .  ) 

n,m  n+l/2,m  n-l/2,m 


+  max  (0 , A  )  -  min  (0,A  )  (V.45) 
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The  maximum  allowable  mass  to  leave  the  cell,  Q  ,  is  then  computed. 

n ,  m  ' 
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4b.  The  following  ratios  are  next  computed  for  use  in  determining  the 
uniting  coefficients. 
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The  limiting  coefficients  are  then  given  by 
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C. 
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47.  The  ant i di f f us i ve  fluxes  in  Equations  V.38  and  V.39  are  limited  by 
multiplying  by  the  limiting  coefficients  and  the  solution  is  advanced  to  the 
next  time  1  eve  1  . 
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48.  The  coding  of  the  f 1 ux-rorrect ed  transport  procedures  is  contained 
in  Subroutine  CONI'  in  the  numerical  model. 

Three  Time  I.evel  F.xp  1  t  c  i  t  scheme 

49.  In  order  to  avoid  the  averaging  of  hydrodynamic  quantities  which  is 
performed  when  employing  a  two  time  level  transport  scheme,  with  a  three  time 
level  velocity  scheme,  a  three  time  level  explicit  scheme  is  considered. 

50.  It  is  instructive  to  observe  the  form  of  the  continuity  equation 
employed  in  the  mul t i opera t iona 1  (alternating  direction)  hydrodynamic  scheme. 
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with 
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where 


At  =  time  step  length 

H*  =  water  surface  elevation  at  intermediate  time  level-' 

k-  i 

H  =  water  surface  elevation  at  time  level  k  -  ]  at  cell  (n,m) 


I 


-  or^  space  increment 

Aa  =  a  space  increment 
k+ 1  ^  ^ 

Un  m+1/2  ~  x  "  velocity  component  at  time  level  k  +  1  at  cell  (n,m) 

k-  1 

un  m+j^2  =  x  "  velocity  component  at  time  level  k  -  1  at  cell  (n,m) 

k+ 1 

Vn+l/2  m  ~  ^  ”  a2  velocity  component  at  time  level  k  +  1  at  cell  (n,m) 

k-  1 

Vn+l/2  m  =  ^  ~  a,i  velocity  component  at  time  level  k  -  1  at  cell  (n,m) 
k 

d  =  water  depth  at  time  level  k  at  cell  n,m 
n ,  m 

If  we  eliminate  the  intermediate  level  ry*  ;  e.g.,  solve  for  ry*  in  Equa¬ 
tion  V.51  and  substitute  in  V.52,  we  obtain: 
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■  It  (  II  ,  III  ) 

k  k 

Since  d  =  n  -  h  '  Equation  V.53  is  a  full  three  time  level  scheme, 
n  ,m  n  ,m  n  ,  m 

In  order  to  develop  a  three  time  level  volume  consistent  transport  scheme,  we 

k  k  k 

replace  d  in  Equation  V.53  with  d  S  ami  thereby  obtain: 
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The  scheme  in  Equation  V.54  employs  centered  space  differencing  of  the  ad¬ 
vert  ive  terms  (FTCS).  Alternatively,  one  may  employ  upwind  differencing 
(FTUS).  The  two  schemes  FTCS  and  FTUS  might  then  he  employed  in  a  f lux— 
corrected  transport  procedure. 

51.  In  order  to  analyze  the  stability  of  these  schemes,  we  consider  tin 

linear  case,  dk  =  d  ,  uk  .„  =  u  ,  vk  .„  =  v  ,  Kk  1  =  K  ,  and 

,  n,m  n,m±l/2  n±l/2,m  a  x 

V  =  K  .  Also  (,,)  =  (m2)  =  1  - 

2  mm 

52.  Let  us  consider  the  following  general  three  time  level  transport 
scheme  written  in  operator  form,  where 
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-  K  62(sk"lN)  -  K  62(sk_1)  =  0  at  (n,m) 
x  x  \  n  ,rn/  y  y  \  n ,m/ 


(V . 55 ) 


with  gf. (-1  ,  0  ,  1) 

Note:  If  g  =  o  ,  central  differencing  of  the  advective  terms  is  affected. 

If  u, v  >  o  g  =  -1  (backward  differencing)  is  employed.  If  u,v  <  o  ,  g  =  + 1 
(forward  differencing)  is  employed.  This  method  constitutes  upwind 
d i f  ferenc ing . 

53.  We  note  that  Equation  V.55  is  a  three  time  level  scheme,  thus  in 

k 

order  to  perform  a  stability  analysis  we  set  up  a  second  variable  ^  and 


63 


set  V 


then  Equation  V.55  may  be  written  in  matrix  form  as  follows 


-2At(uT  +  vT  )  +  2At[K  S2 
x  y  V  x  x 


+  K  +  1  Sk 

v  v  /  n  ,  m 


(V.56) 


54.  Equation  V.56  is  in  the  required  format  for  a  stability  analysis. 

„  ,  .  .  „k  kaAt  nyAv  mflAx 

Following  standard  practice,  assume  S  =  e  e  e 

n  ,  m 

55.  Thus  we  may  develop  the  following  supplemental  relations: 

T  (Sk  )  s  (*-  sin  PAx  -  Is  sin2  &f*)sk 
x  n ,  m  \  Ax  Ax  2  /  u ,  m 


T  (Sk  )  :  (  -J-  sin  yAy  -  sin2  ^)sk 

y  n,m  \ Ay  s  y  Ay  2  /  n,i 


,2,ck  -4  .  2  pAx 

6  ( S  )  i  -r  sin 
x  n ,m  .2  2 

Ax 


c2,ck  .  ..  -4  .2  pAx 

<5  (S  sin 

y  n,m  Ay2  2 


( V .  5  7  ) 


Thus  we  obtain  for  Equation  V.56,  the  following  matrix  equation: 


•  (  "I  I -  •  '  ''ll'  ’’  ')  (i ..  M 


>i  ■  -  •  >  hi 


I  S  I 

I  n. in  I 


( V .58) 


C  =  uAt/Ax  C  =  vAt/At  D 

x  y  x 


.  ...,ij 


In  oriter  to  achieve  stability  all  the  eigenvalues  of  A  must  be  less  than  1 
in  magnitude.  To  compute  the  eigenvalues,  the  following  chu rac t e rt i s t i c  equa- 
t i on  i s  used . 
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(A  -  AI|  =  0 


A**  +  A 


C^2i  sin  pAx  -  4g  sin2  +  C^^2i  sin  Y^y  "  4g  sin2 

(8D  sin2  +  80  sin2  = 

\  x  2  y  2  J 


(V. 59) 


-1  + 


In  considering  Equation  V.59,  let  us  introduce  the  following  supplemental 
variables : 


=  Cx(2i 


sin  |3Ax  -  4g  sin 


2  §4x 


)  +  cv(; 


2i  sin  yAy  -  4g  sin 


i  or  2  PAx  .2  yAy 

h  =  80  s l n  —  +  0  sin 

x  2  v  2 


( V . 60 ) 


Then  Equation  V.59  becomes 


A  +  2a  +  (b  -  1)  =  0 


(V.61) 


Using  the  quadratic  equation  formula: 


A 


1,2 


(V.62) 


56.  To  find  restrictions  on  C 


and  0  such  that 


|  A  j  <  1  appears  to  be  difficult  for  this  case. 

*  »  ^ 

57.  Instead,  we  consider  the  ranges  of  C  ,  C  ,  0  ,  and  0  for 

x  y  x  y 

the  Mississippi  Sound  case.  Consider  the  following  maximum  velocity  and 
dispersion  conditions  in  Mississippi  Sound:  u  =  v  =  3  fps , 

=  100  ft2/sec.  Since  on  the  global  grid  the  minimum  cell  dimensions  are 
Ax  =  Ay  =  3500  ft  and  a  time  step  At  =  360  sec  is  employed,  bounds  on  the 


supplemental  relations  in  Equation  V.53  become  C 


C  <  0.3 
v 


and 


D  ,  D  <  0.003  . 
x  y 

58.  Let  us  introduce  8  =  2n/mAx  and  y  =  2n/nAv  and  consider  the 

m  n 

wave  numbers  n  and  m  to  vary  from  2-9  over  3  Jog  cycles.  We  compute  the 
eigenvalue  for  each  set  of  wave  numbers  a  and  m  as  follows: 
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a  =  C  I2i  sin  —  -  4g  sin  -  1  +  C  ( 2 i  sin  —  -  4g  sin  -  ) 
x\  m  &  m/  y  \  n  6  n/ 

2ti  .2  n\ 

b  -  81 D  sin  -  +  D  sin  -  I 
\  x  m  v  n  / 


(V. 63) 


“ !  V  (!)  - 6 


59.  The  computer  program  employing  Equation  V.62  shown  in  Table  V-3  was 
used  to  compute  the  eigenvalues  for  the1  general  transport  scheme.  The  results 
are  shown  in  Table  V-4  for  the  cases  considered.  Case  III  represents  the  con¬ 
ditions  to  be  considered  within  Mississippi  Sound.  Since  the  upward  scheme  is 
unstable,  it  may  not  be  employed  as  the  lower  order  scheme  in  the  flux- 
corrected  transport  method. 

60.  Although  it  is  possible  to  flux-correct  the  three  time  level  center¬ 
ing  of  the  advective  terms  scheme,  this  was  not  undertaken  in  this  study. 

This  three  level  scheme  with  the  advective  terms  centered  in  time  was  coded  as 
Subroutine  CONCE  within  the  model,  thereby  allowing  the  model  user  an  alternate 
scheme  to  employ  in  the  transport  calculations.  This  scheme  may  be  employed 
for  transport  simulations  in  which  the  constituent  levels  are  reasonably  uni¬ 
form.  For  sharp  front  problems,  the  Flux-Corrected  Transport  Scheme  should  be 
used . 


Dispersion  Coefficient  Formula ti on 

61.  To  close  the  numerical  approximations  to  the  two-dimensional  depth 
averaged  transport  equation,  relations  for  the  effective  dispersion  coef¬ 
ficients  must  be  developed  in  terms  of  flow  field  properties. 

62.  The  effective  dispersion  coefficients  are  assumed  to  have  the  fol¬ 
lowing  forms: 


~  r  ulh 

K-  -  C  Vg  — _  — 
x  x  C 


D  ;  K*  =  C  Vg  +  D 

x  y  y  C  y 


(V . 64) 


whe  re 


K*  ,  K* 
x  y 


=  effective  dispersion  coefficients  in  the  x  and  y  directions 

-  acceleration  due  to  gravity 

-  velocity  components  in  the  x  and  y  directions,  respectively 


Table  V-3 

Computer  Program  for  Determination  of  the  Eigenvalues 
for  the  General  Three  Time  Level  Explicit  Scheme 


PROGRAM  EIGEN 
PARAMETER  (ND=5 , N6=6 ) 

COMPLEX  A , B , LM1 , LM2 
DATA  PI/3. 141592654/ 

READ(ND , 1 ) I  CYC , CX,CY , DX , DY 
FORMAT ( 1 5, AF 10.0) 

I F ( I  CYC . LE . 0 ) STOP 

WR  ITE(N6 , 250) 1  CYC , CX , CY ,DX , DY 

FORMAT ( 1H 1  )  , 56X ,  -'  THREE  «  1  ME  LEVE L  TRANSPORT  SCHEMES* ,////, 

1  25X  ,  1 5  WAVE  NUMBER  CYCLES--,/, 

2  25X ,  G8 . 3  ,  X  COURANT  NUMBER'--- ,  /  , 

2  25X,G8 . 3  Y  COURANT  NUMBER-’-,/, 

2  25X ,  G8 . 3  X  DIFFUSION  NUMBER  ••,/, 

2  25X,G8,3m*  Y  DIFFUSION  NUMBER---,/) 

IF (CY+DY ) 3 , 3 , A 

ICYC1=I 

L I M 1  =  1 

LIM2=1 

GO  TO  5 

ICYC1=ICYC 

LIM1=2 

L1M2=9 

G=-  1  . 

WRITE(N6,56)G 

FORMAT ( 1 H 1 , 20X , *G=* , F6 . 3 , / ) 

DO  200  1  =  1,  ICYC 1 
DO  200  J=LIM1 , LIM2 
N= J*  1 0---*  ( I  -  1 ) 

DO  200  K= 1,1  CYC 
DO  200  L=2 , 9 
M=  L*  1 0**  ( K- 1 ) 

A=CMPLX  ( -  A  ,  *G*  ( CX--S  1  N  ( P  [  /  M  )-•■  S I N  ( P  I  /  M ) +C  Y-'-S  I N  (  P I  /N  )*S  1 N  ( P I  /  N  ) ) , 

1  2 ,*CX*S I N ( 2 . *P J /M )+CY* S 1 N ( 2*P I /N ) ) ) 

B=CMPI.X  (  1  .  -8  .  DX’’- S I N  ( PI  /M  )*S  I N  (P I  /M )  -'-DY'-’-S  1  N(  P 1  /  N  )-'-S  I N  (P I  /  N ) )  ,0. 

LM I  =CSQRT (  (A/2.  )-'  (A/2.  )+B) 

LM2=-A/2 . -LM1 
LMl=-A/2 ,+LMl 
AMAG 1 =CABS ( LM 1 ) 

AMAG2=CABS ( LM2 ) 

V/R  I TE  ( N6 . 5  5  )  N  ,  M ,  AMAG  1  ,  AMAG2 
FORMAT ( 20X , 2  J  6 , 9X , EG , 3 , 1  AX , F6 . 3 ) 

IFCG.GT. .5 )GO  TO  5000 
G=G+ 1 . 

GO  TO  5005 


Eigenvalue  Analysis  Results  for  the  General 
Three  Time  Leve 1  Explicit  Scheme 


Magnitude  of  Largest 

Advect  i  ve  Term  _ Eigenvalue _ 

Case  I:  C  =  C  =  0.  D  =  £>  =  0.125  (Diffusion  Only) 

_ x _  y _ x  _ y _ 

Centered  (g  =  0)  <1 

Upwind  (g  =  -  1 )  <  1 

Forwa  rd  ( g  =  +  1 )  <3 

Case  II:  C  =  C  =0.5  D  -  D  =0.125 

_ x _ y _  x _ y 


Centered  (g  =  0)  1 

Upwind  (g  =  -  1 )  >1 

Forward  (g  =  +1)  >1 

Case  III:  C  =  C  =0.3  D  =  D  =  0.003 

_ x _ y _ x _ y _ 


Centered  (g  =  0) 
Upwind  (g  =  -1) 


<1 

>1 


h  =  water  depth 
C  =  Chezy  coefficient 

C^  ,  0^  =  dispersion  factors  in  the  x  and  v  directions,  respectively 

D  ,  D  =  dispersion  offsets  due  to  wind  effects  in  the  x  and  y 
x  y 

directions,  respectively,  (I)  ,  D  *  0) 

x  y 

63.  Elder  (1959)  has  determined  the  longitudinal  and  lateral  dispersion 
coefficients  in  open  channel  flow  experiments  to  he  given  by  the  following 
relations 


K  =5.93  hu*  Kt  „  =  0.23  hu*  (V.65) 

L  la 

where 

K  =  lateral  dispersion  coefficient 

Ln 

K  =  longitudinal  dispersion  coefficient 

L 

h  =  water  depth  (hydraulic  radius) 
u*  =  shear  (friction)  velocity 
For  open  channel  flow,  the  following  relations  hold: 

u*  =  \lghS  u  =  CjhS^  (V.66) 

where 

u*  =  friction  velocity 
g  2  acceleration  of  gravity 
h  =  water  depth 

=  slope  of  energy  grade  line 
C  =  Chezy  coefficient 
As  a  result,  we  obtain: 


“*  -  ^  C 

where 

u  =  velocity 
g  =  gravity 
C  =  Chezy  coefficient 
Therefore,  Equation  V.65  becomes 


h  =  5  f  kla 


r  uh 


(V.  67) 


Cm  c.  q  \ 
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64.  Taylor  (1954)  has  conducted  pipe  flow  experiments  to  determine  the 
longitudinal  dispersion  coefficient.  By  assuming  the  hydraulic  radius  as  half 
the  pipe  radius  in  the  pipe  experiments  and  equal  to  the  water  depth  m  a  uni¬ 
form  steady  flow  open  channel,  the  coefficient  in  the  longitudinal  dispersion 
coefficient  equation  (V.68)  is  determined  to  he  20.2  rather  than  5.93. 

65.  F.r  one-dimensional  flow  in  the  x  direction,  C  f.  (5.93,  20.2) 

x 

and  C  -  0.23  .  In  order  to  extrapolate  these  results  to  a  two-dimensional 
flow  as  occurs  in  Mississippi  Sound  ,  C  ,  D  ,  D  are  specified  for 

each  grid  cell.  The  cel!  fate  conditions  for  each  cell  are  examined  indepen¬ 
dent  Iv  m  each  coordinate  direction.  For  a  no-flow  tell  face  condition,  C 

x 

.a til  IJ  or  C  and  1)  are  set  to  zero.  (No  dispersion  may  occur  across  a 
x  v  y 

solid  boundary.)  For  a  (low  cell  face  condition  on  the  u-velocity  cell  face 
l see  Figure  A - 1  ) ,  the  second  digit  in  the  u-velocity  cell  face  flag  is  ex¬ 
amined.  If  advection  is  not  allowed,  is  reduced  by  a  user  specified  fac¬ 

tor,  F  .  For  a  flow  cell  face  condition  on  the  v-velocity  cell  face  (see 
Figure  A-  I  )  ,  analogous  ptocedures  are  employed. 

06 .  In  order  to  calibrate  the  effective  dispersion  coefficients,  C^_  , 

F  .  I  ,  It  specified  on  a  cel  1 -by-ce 1 1  basis,  and  F  are  adjusted  until 

x  y  v 

simulated  salinity  levels  correspond  to  measured  salinity  levels.  Based  upon 

experimental  results,  ('  ,  C  f.  (5.93,  20.2  )  F  -  0.23/5.93  =  0.0388  .  The 

x  y 

wind  effect  terms  D  and  II  .ire  more  difficult  to  determine.  Leendertse 

x  V  > 

(1970,  1971)  in  a  simul  it  ion  of  Jamaica  Bav  suggested  D  ,  D  >;  (25,  45)  ft*”/ 

X  V 

sec.  Since  the  velocity  magnitudes  u  and  v  in  Kquation  V.64  will  in- 

i  cease  with  increasing  wind  speeds,  the  effective  dispersion  coefficients  will 

i  n<  i  e  ise  without  considering  0  and  1'  .  F..|  this  reason,  I)  and  I) 

x  v  x  y 

w  i  !  !  he  set  to  w  r>  > ,  and  t  ,  F  ,  and  I  will  he  ad  |  us  t  ed  during  the  (..ill- 


PART  VI:  DEVELOPMENT  OF  A  GLOBAL  GRID  FOR  MISSISSIPPI  SOUND 


67.  The  hydrodynamic  salinity  model  developed  in  this  study  employs  a 
variable  (exponentially  stretched)  computational  grid.  In  transforming  real 
space  to  computational  space,  each  coordinate  axis  is  mapped  independently 
using  Program  MAPIT.  This  enables  the  specification  of  partition  points  along 
each  axis.  These  partition  points  coincide  with  the  location  of  barrier  island 
passes  and  the  shoreline. 

68.  The  grid  spacing  along  the  southern  and  eastern  extents  of  the 
global  grid  as  shown  in  Figure  VI-1  was  set  to  correspond  to  the  15-min  lati¬ 
tude  and  longitude  uniform  grid  employed  in  the  GTM.  In  this  manner,  inter¬ 
polation  in  only  one  space  dimension  is  required  in  developing  the  global 
grid  boundary  conditions.  This  consideration  established  an  upper  bound  on 
the  grid  spacing  in  each  direction.  The  lower  bound  in  grid  spacing  was  ob¬ 
tained  such  that  grid  cell  aspect  ratios  would  be  no  larger  than  20  and  ex¬ 
plicit  time-step  limitations  would  not  be  severe.  Considerable  engineering 
judgment  and  consultation  with  the  Mobile  District  was  exercised  in  developing 
the  global  grid. 

69.  The  grid  developed  extends  in  geographical  extent  from  the  western 
end  of  Lake  Borgne  to  Santa  Rosa  Island  in  the  west-east  orientation  and  from 
the  Mississippi  River  Delta  to  above  Mobile  Bay  in  the  north-south  orientation. 
The  east -west  extent  employs  116  lines  and  the  north-south  60  lines,  resulting 
in  a  grid  of  6789  (115  *  59)  computational  cells  as  shown  in  Figure  VI-2.  Min¬ 
imum  spatial  resolution  of  approximately  3500-4000  ft  is  obtained  within  the 
[lasses  into  Mississippi  Sound  and  within  the  Sound  itself.  Depths  within  Mis¬ 
sissippi  Sound  are  relatively  shallow  (10-20  ft),  except  in  the  navigation 
channels,  which  are  normally  maintained  at  30-35  ft.  As  a  result,  the  gravity 
wave  speed  within  the  sound  is  <38  fps ,  resulting  in  an  explicit  time-step  lim¬ 
itation  of  approx ima te 1 y  100  sec.  All  numerical  simulations  employ  a  360-sec 

( 0  nun)  time  step,  resulting  in  a  maximum  spatial  Courant  number  of  less  than 
U  within  the  Sound.  The  mapping  for  the  horizontal  (east-west)  orientation  is 
presented  in  Table  VI-1.  The  mapping  for  the  vertical  (north-south)  orienta¬ 
tion  is  given  in  Table  VI-2.  The  grid  obtained  is  presented  in  Figure  Vl-2. 

In  Tables  VI-1  and  VI-2,  the  real  space  values  are  in  map  inches  as  measured 
from  the  maps  shown  in  Table  VI-3  below. 
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Horizontal  Mapping  Results 


1:250000 
1:250000 
1  .-250000 
1 : 250000 
1:250000 


Ini t i  a  l  Depth  Assignment 

70.  Program  TGRID  was  employed  to  plot  snbgrids  corresponding  to  grid 
areas  shown  on  the  individual  nautical  charts  as  shown  in  Table  Vl-4.  Each 
subgrid  was  next  overlaid  on  its  corresponding  nautical  chart.  For  most  grid 
cells,  the  assigned  depth  represented  an  average  depth  over  the  cell.  In 
order  to  represent  flow  restrictions  in  some  cells,  the  minimum  depth  in  the 
cell  was  more  heavily  considered  in  the  averaging  process.  The  assigned 
depths  are  presented  in  Table  VI-5,  which  were  directly  output  from  WIFM.  In 
interpreting  Table  VI-5,  all  water  depths  are  in  feet,  preceded  by  a  minus 
sign,  and  are  with  respect  to  local  mean  sea  level  (LMSL),  which  was  taken 

1  ft  above  MLLW.  (Land  elevations  were  assigned  a  value  of  +10.) 

I n corporation  of  Hydrographic  Su r v e y  1  n f o rm a t_i ori 

71.  In  order  to  properly  model  circulation  within  Mississippi  Sound,  the 
depth  field  must  be  accurately  specified  along  transects  across  barrier  island 
passes  and  bay  entrances.  Raytheon  Ocean  Systems  was  contracted  by  the  Mobile 
District  to  obtain  soundings  at  50-ft  intervals  for  the  transects  shown  in  Fig¬ 
ure  1 1 -5.  Shallow  water  areas  and  land  obstructions  were  avoided  and  soundings 
were  plotted  in  state  (Louisiana,  Mississippi,  and  Alabama)  coordinate  systems 
at  a  scale  of  1:8000. 

72.  The  method  of  incorporating  the  hydrographic  survey  data  into  the 
previously  developed  depth  field  is  developed  in  the  context  of  a  journal  as 
presented  below. 

Journal  setup 

73.  Plots  of  transects  were  identified  and  labeled  with  respect  to  the 
general  location  map  (Figure  f I - 5 ) .  Requested  transects,  noted  on  small  craft 
charts  11367,  72,  74,  and  78,  were  contrasted  to  actual  transects  by  checking 
beginning  and  ending  points  in  state  coordinates  of  transects  on  the  plots  and 
on  the  charts.  Iri  cases  where  differences  did  exist,  actual  transects  were 
added  to  the  charts.  Latitude  and  longitude  of  end  points  of  the  surveys  were 
then  interpolated  from  the  small  craft  (harts. 

74.  A  journal  was  prepared  (Table  Vl-6)  identifying  assigned  labels  for 
the  various  transects,  giving  la t i tude/ 1 ougi t ude  and  state  coordinates  of  the 
end  points  of  each  transect,  and  identifying  pertinent  nautical  charts. 
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Table  VI-4 


Progra in  TGRID  Subg r  1  d  P  1  o ts 


’ 'm 

Grid 

I nd i cps 

M 

N 

Naut i ca ] 

Cha  rt 

Scale 

■■  i 

i 

1  5- 75 

1-50 

1  137  1 

22nd 

F.d  . 

4/80 

1  :  80000 

■  -i 

•'  ‘  At 

i 

36-  5  6 

1-30 

1137  1 

22nd 

Ed  . 

4/80 

1 : 80000 

11363 

20th 

Ed. 

1/80 

1 :80000 

# 

1  1  364 

2  3  r  d 

F.d  . 

2/80 

1  :  80000 

37-56 

1  -  50 

1  1  561 

40  th 

F.d  . 

ON 

'•'J 

oc 

1  :  80000 

1  1363 

20th 

Ed  . 

1/80 

1  :  80000 

1  1364 

23  rd 

Ed. 

2/80 

1  :  80000 

• 

15-30 

51-78 

1  1373 

25th 

F.d. 

6/80 

1:80000 

57-59 

31-115 

1  1  360 

24 1  h 

F.d  . 

2/80 

1  .-456394 

37-56 

29-45 

11363 

20th 

Ed  . 

1/80 

1:80000 

• 

1  -  <M 

79-109 

11376 

33  rd 

F.d. 

9/79 

1  : 80000 

15  -  56 

1  10-  1  15 

11382 

23  rd 

Ed . 

1  1/77 

1  : 80000 

Z  I  - 


hie  Survey  Journal 


Procedure  for  transferring 
grid  lines  to  sounding  plots 

75.  The  various  documents  used  and  their  scales  are:  (a)  small  craft 
charts  -  1:40000;  (b)  nautical  charts  -  1:80000;  (c)  sounding  plots  -  1:8000; 
and  (d)  subgrids  -  1:80000.  All  measurements  were  taken  using  a  scale  divided 
into  sixtieths  (1/60)  of  an  inch. 

76.  The  transects  as  actually  surveyed  were  transferred  to  the  nautical 
charts  from  the  small  craft  charts  using  la t i tude/ longi tude  in  combination 
with  common  land  feature  identification,  the  2:1  scale  relationship  serving  as 
a  control  for  accurate  placement. 

77.  The  nautical  charts  were  overlaid  with  the  corresponding  subgrid  to 

* 

identify  the  affected  cells.  A  1:10  scale  factor  was  used  to  transfer  the 
grid  lines  to  the  sounding  plots.  Each  cell  was  then  identified  and  labeled 
by  its  grid  coordinates.  As  each  nautical  chart  contained  several  transects, 
it  was  possible  to  identify  in  terms  of  state  coordinates  at  least  one  lateral 
grid  line  to  serve  as  a  control  for  consistent  lateral  line  placement.  Ver¬ 
tical  line  placement  was  controlled  by  visual  comp  rison  of  the  chart  and 
sounding  plot. 


Comparison  of  Previously  Assigned  Depths  to  Surveyed  Depth 

78.  The  nautical  chart  with  subgrid  overlay  and  depth  grid  (Table  VI-5) 
was  used  to  refine  cell  identification  with  regard  to  model  features  to  ac¬ 
count  for  related  cells  and  special  cases. 

79.  The  overlaid  nautical  charts  had  been  used  in  assigning  initial 
cell  depths,  thus  their  use  in  comparing  survey  depths  with  previously  as¬ 
signed  depths  would  serve  as  a  control  so  that  the  whole  cell  could  be  viewed 
in  conjunction  with  the  transect . 

80.  The  depth  grid  (Table  VT-5)  was  used  to  list  individual  cell  depths 
along  each  transect.  Then  each  sounding  plot  was  examined  to  determine  a  rep¬ 
resentative  depth  range  in  each  cell.  The  representative  range  was  compared 
with  the  previously  assigned  depth  together  with  the  specific  cell  on  the 
overlaid  nautical  chart  to  determine  if  a  cell  depth  needed  to  be  corrected. 

A  table  was  prepared  to  show  the  results  of  these  comparisons  (Table  VI-7). 


Table  Vl-7 

Comparison  of  Initial  Cell  to 
Survey  Range  Depths  (ft) 


Initial  Depth 

Revised  Depth 

Transect  Grid  Cell 

(MLLW) 

Depth  Range  Survey 

(MLLW) 

1A 

(15,39) 

7 

7.0-31.6-  5.7  (P) 

+ 

IB 

(15,38) 

9 

11 . 1-12.4-13. 1 

1 1 

IB 

(16,37) 

9 

7.  1-12.4 

+ 

1C 

( 15 , 34)* 

10 

9.0-35.0 

+ 

1C 

(16,34)* 

14 

9.0-35.0 

20 

1C 

(16,35) 

1 1 

9.5-16.0-14.2 

12 

2 

(21,24) 

8 

4.9-  8.2 

+ 

3 

(25,34) 

4 

9.6-15.0-  9.0  (P) 

+ 

3 

(26,33) 

10 

18.4-  7.0-  9.6 

+ 

3 

(26,32) 

19 

17. 1-22.0-18.4 

+ 

3 

(26,31) 

13 

13.7-21.0-16.7 

18 

3 

(26,30) 

6 

11.0-13.7  (P) 

+ 

4 

(26,34) 

10 

6.6-13.4-10.6 

1 1 

4 

(27,34) 

13 

10.8-16.0-13.0-17.0 

14 

4 

(28,33) 

20 

6.3-30.0 

+ 

5 

(31,29)* 

8 

9.5-11.2  (P) 

+ 

5 

(31,30)* 

12 

11.4-14.7 

+ 

5 

(32,30) 

1 1 

14.5-16.9 

14.5 

5 

(33,31) 

17 

18.0-24.3 

+ 

5 

(34,31) 

27 

24.4-36.8-32.7  (P) 

+ 

6 

(33,32) 

20 

24. 1-32.5 

22.5 

6 

(34,32) 

24 

32.2-21.4 

+ 

7 

(38,30) 

16 

3.9-  9.3-  2.8 

6 

7 

(38,31) 

10 

3.9-  9.3-  2.8 

6 

(39,30)** 

2 

3.9-  9.3-  2.8 

6 

(39 ,31)** 

18 

3.9-  9.3-  2.8 

6 

8T 

(42,16) 

4 

3.9-  7.8 

+ 

8t 

(43 ,16) 

4 

7.3-16.7-  5.5 

8 

10A 

(42,29) 

10 

7.2-13.6 

+ 

(43,29) 

10 

7.8-  9.7 

+ 

(44,29) 

10 

8.6-26.3 

12 

(45,29) 

16 

25.5-  7.0-16.8 

+ 

(Continued) 

Note : 

MLLW  -  Model 

Depth  +1 ;  (P) 

-  Partial  transect  of  cell;  + 

-  Indicates  no 

revision  in  previously  assigned  model  depth. 

* 

Indicates  cells  considered 

together . 

** 

Related  cells 

t 

Special  case, 

transect  8; 

survey  halved. 

(Sheet  1  of  3) 


Table  VI-7  (Continued) 


T  ransect 

Grid  Cel  1 

Initial  Depth 
(MLLW) 

10A 

(46,29) 

18 

(Cont ' d ) 

(47,29) 

8 

(48,29) 

1 1 

(49,29) 

11 

10B&C 

(42,28) 

1 1 

(43,28) 

6 

(44,28) 

12 

(45,28) 

10 

1 0B&C 

(46,28) 

6 

10C 

(47,28) 

7 

(48,28) 

4 

(49,28) 

18 

11A 

(59,31) 

20 

(60,31) 

18 

(61,31) 

12 

(62,31) 

16 

1  IB 

(59,30) 

10 

(60,30) 

10 

(61,31) 

12 

(62,31) 

16 

1 1C 

(59,30) 

10 

(60,30) 

10 

(61,30) 

12 

(62,30) 

1 1 

12A 

(71,31) 

9 

(72,31) 

7 

(73,31) 

6 

(74,31) 

5 

(75,31) 

10 

(76,31) 

18 

12B 

(71,31) 

9 

(72,31) 

7 

(73,31) 

6 

(74,31) 

5 

(75,30) 

12 

Revised  Depth 


Depth  Range  Survey _  (MLLW) 


17.0-19.0-  9.8-  5.8  16 

5.8- 10.7  + 

11.0-23.1  12 

11.3- 23.3  + 

4.1-  6.5  (B)  5 

5.9- 14.6  (C) 

7.2- 18.9-  9.4  (B)  9.5 

5.9-  4.5-16.9-13.4  (C) 

17.1- 34.3-  7.1  + 

8.5- 20.7-13.5  (C) 

16.2-  6.8-10.1  (B) 

8. 1-  8. 7  7.5 

8.1- 10.2  (B)  8 

8.2-  8.5  (C) 

6. 1-  8. 1  + 

6.2-  4.1-  6.6  5.5 

6.9- 32.7-17.6  + 

13.6-16.4  (P)  + 

16.4- 9.6  + 

6.9- 10.8  + 

6.5- 46.5  + 

4.0-  8.0-  5.0  6 

5.8- 13.9-10.0  + 

6.9- 10.8  + 

6.5-46.5  + 

4.0-  8.0-  5.0  6 

5.8- 13.9-10.0  + 

11.9-12.9-11.9  + 

5. 1- 11.4  13.9-30.7  + 

3.2- 11.5  13.7-16.4  + 

7.3- 11.5  10 

8.9- 17.0  8 

8.3- 10.9  10 

11.5- 19.9-  9.0  13 

11.4- 20.8  (P)  + 

3.2- 11.5  13.7-16.4  + 

7.3- 11.5  10 

8.9-17.0  8 

8.3- 10.9  10 

10.4- 19.0  + 


(Continued) 
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Table  VI-7  (Concluded) 


Initial  Depth 

Revised  Di 

T  caused 

Grid  Cell 

(MLLW) 

Depth  Range  Survey 

(MLI.W) 

12A 

(71,31) 

9 

3.2-11.5  13.7-16.4 

+ 

(72,31) 

7 

7.3-11.5 

10 

(73,30) 

12 

17.6-17.8  9.0-10.6 

10 

(74,30) 

12 

9.8-13.0 

+ 

(75,30) 

12 

10.4-19.0 

+ 

1  2D 

(70,31) 

b 

7.0-  8.1 

+ 

(71,31) 

9 

3.2-11.5  13.7-16.4 

+ 

(72,30) 

15 

14.9-17. • 

1  7 

(73,30) 

12 

17.1,-17.8  9.0-10.6 

10 

(74,29) 

1  7 

16  '  - i /  .  5 

t 

(7.3,29) 

16 

7.7-18.4  i  1 1  ! 

+ 

]  5 

(82,23) 

4 

j  .  7-  0  .  -t  (Pi 

+ 

(82,24) 

b 

6.3-  7.5 

1 

(82,25) 

6 

7.2-  8.3 

7 

/ 

(82,26) 

7 

7.9-  8.5 

8 

(82,27) 

8 

8.4-10.0 

+ 

(82,28) 

6 

3.9-11  .0 

+ 

14 

(86,25) 

4 

.8.5-14.8 

7 

15 

(90,28) 

6 

9.7-13.4  ( P ) 

10 

(91,29) 

8 

9.6-12.6 

10 

(92,29) 

21 

12.6-  9.9  14.1-41.2 

16 

(93,30) 

19 

20.6-42 . 2 

5  3 

(91,32) 

7 

8.6-15.0-10.0 

10 

(92,32) 

34 

10.5-  4.0-14.8 

+ 

(93,32) 

8 

14.0-  9.0-15.8 

9 . 5 

17 

(91 ,34) 

1  7 

23.5-20.0-42.6 

20 

(92,34) 

22 

44.5-17.7-20.8 

+ 

(93,34) 

21 

20.0-31.5  (P) 

+ 

Barrier  Island 


Conf i gurat ion 

81.  All  harrier  islands  were  located  on  cell  faces  with  a  land  eleva¬ 
tion  of  10  ft.  The  barriers  were  modeled  as  exposed  barriers;  no  overtopping 
occurred  at  any  barriers  in  any  of  the  simulations.  The  barrier'  islands  are 
given  in  Table  VI-8.  The  orientation  labeled  1  indicates  a  u-taie  barrier, 
while  the  orientation  labeled  2  represents  the  barrier  obstruction  to  be  along 
the  v-face  of  the  cell  as  shown  in  Figure  A- I . 

Flow  inputs 

82.  The  major  flow  inputs  (see  Table  VI-9)  and  their  locations  in  the 
grid  will  he  considered.  Average  daily  flow  values  as  obtained  from  the 
United  States  Geological  Survey  (USGS )  are  employed  in  the  numerical  model. 

CaJ_i_b  ra  1 1 on/ V'e r  l  f  i  ca  t  i  on  S  t  a  t  i  ons 

83.  In  order  to  compare  simulated  water  surface  elevations  with  the  ob¬ 
served  values,  the  location  of  the  water  surface  elevation  stations  must  be 
located  on  the  grid.  Results  are  shown  in  Table  V 1  - 1 0 . 

84.  The  placements  of  the  meteorological  and  velocity/salinity  stations 
are  given  in  Tables  VI-11  and  VI-12,  respectively. 


Table  Vl-8 

Bar r ier  Co n f igu rat  l on 


Lora  t i on 

•  i ent  at i on 

N 

M 

1 

41 

17 

2 

41 

18 

1 

42 

18 

1 

43 

18 

1 

44 

18 

1 

45 

18 

■ ) 

45 

i‘< 

1 

34 

3  1 

1 

1 

15 

3  b 

51 

i  i 

] 

n6 

i 

1 

>  / 

.30 

! 

4  i 

28 

1 

r>0 

28 

1 

51 

28 

1 

52 

28 

1 

53 

28 

1 

55 

29 

1 

56 

29 

') 

56 

30 

1 

5  7 

30 

I 

38 

30 

1 

58 

24 

5  7 

24 

I 

74 

1  7 

1 

6  i 

31 

1 

b4 

31 

I 

65 

31 

') 

65 

32 

1 

66 

32 

1 

67 

32 

1 

69 

31 

') 

78 

19 

2 

78 

20 

2 

78 

21 

2 

78 

22 

2 

41 

41 

1 

43 

43 

2 

44 

45 

1 

42 

50 

( Cont  i  nne<l ) 
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Table  VI-8  (Continued) 


Location 


No. 

Orientation 

N 

M 

41 

1 

41 

49 

42 

2 

40 

49 

43 

1 

39 

49 

44 

2 

38 

59 

45 

2 

15 

33 

46 

1 

30 

29 

47 

2 

30 

28 

48 

2 

30 

29 

49 

2 

17 

43 

50 

1 

19 

43 

5  1 

2 

19 

45 

52 

2 

19 

46 

53 

2 

19 

4  7 

54 

2 

19 

48 

55 

1 

2  1 

44 

56 

2 

21 

44 

57 

2 

22 

41 

58 

2 

O  '/ 

42 

59 

2 

2  2 

45 

60 

2 

22 

46 

61 

0 

2  l 

47 

62 

21 

48 

63 

i 

21 

48 

64 

i 

22 

47 

65 

i 

23 

47 

66 

i 

24 

47 

67 

i 

22 

46 

68 

i 

23 

49 

69 

i 

24 

49 

70 

i 

25 

49 

71 

i 

24 

45 

72 

i 

25 

35 

73 

2 

24 

37 

74 

2 

24 

43 

75 

2 

24 

44 

76 

2 

23 

44 

77 

2 

23 

45 

78 

2 

2  3 

46 

79 

2 

23 

4  7 

80 

1 

76 

29 

( Coni  i  rilled  ) 
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Table  VI-8  (Concluded) 


Table  VI-9 

Flow  Inputs  in  the  Global  Grid 


Inflow 

Grid  Cell 

Mobile  River  System 

(97,3) 

East  Pascagoula  River 

(59,19) 

Pascagoula  River 

(59,17) 

Pearl  River 

03,33) 

Jourdan  and  Wolf  Rivers 

(19,20) 

Biloxi  River  System 

(34,15) 

Table  VI-10 

Water  Surface  Elevation  Stations 


Station  No 


WIFM  Grid  Coordinates 


8* 

S 

(110,28) 

12 

s 

(107,24) 

13 

s 

(  90,28) 

10 

s 

(  88,14) 

9 

s 

(  87,27) 

'll 

s 

(  92,  5) 

7 

s 

(  70,21) 

6 

s 

(  67,32) 

'5 

s 

(  61,21) 

'4 

s 

(  55,29) 

’21 

s 

(  47,17) 

'20 

s 

(  43,15) 

’15 

s 

(  43,17) 

16 

s 

(  35,18) 

’19 

s 

(  34,15),  (35,15) 

’14 

s 

(  35,31) 

'17 

s 

(  24,23) 

18 

s 

(  21,22) 

'22 

s 

-- 

'1 

s 

(  32,59) 

'2 

s 

(  26,57) 

3 

s 

(  22,51),  (23,51) 

22 

DSP 

(  67,49) 

’23 

DSP 

(  89,49),  (90,49) 

’24 

DSP 

(108,49) 

S  =  Standard  Gage;  DSP  =  Deep  Sea  Pressure 
Gage. 

Refer  to  Figure  FI-3  for  station  locations. 


Table  VI-11 

Meteorological  Stations 


Station  No. 

Ml 

M2 

M3 

M4 

M5 


WIFM  Grid  Coordinates 

(22,23) 

(34,31) 

(42,17) 

(55,29) 

(89,28) 


Table  VI-12 

Velocity/Salinity  Continuous  Stations 


Station  No. 

Type 

WIFM  Grid  Coordinates 

VI 

V 

(15,39)* 

V2 

V 

(16,37),  (16,38) 

V3 

V 

(16,35) 

V4 

V 

(27,32),  (26,31) 

V5 

V 

(29,22) 

VS6 

VS 

(28,27),  (27,27) 

V  7 

V 

(32,25) 

V8 

V 

(34,32) 

V9 

V 

(45,29) 

VS10 

VS 

(49,29),  (49,28) 

Vll 

V 

(49 ,23) 

VS  12 

VS 

(54,27) 

V13 

V 

(62,27) 

VI 4 

V 

(60,30) 

VS  15 

VS 

(62,31),  (62,32) 

V16 

V 

(67,28),  (68,28) 

VI 7 

V 

(72,30),  (72,31) 

VS  18 

VS 

(75,30) 

V19 

V 

(76,25) 

VS20 

VS 

(87,25) 

V21 

V 

(93,30) 

PART  VII:  SELECTION  OF  THE  CALIBRATION  AND 
VERIFICATION  PERIODS 


85.  In  order  to  understand  the  general  behavior  of  the  hydrodynamics 
and  salinity  distribution  in  Mississippi  Sound  over  the  sampling  period,  river 
inflows,  meteorological  station  wind  information,  and  salinity  transect  data 
were  tabulated.  A  brief  discussion  of  the  significant  findings  and  the  tables 
themselves  are  presented  in  turn  for  each  data  group.  Periods  for  numerical 
study  are  selected  based  upon  these  findings  and  on  modeling  requirements. 

The  selection  process  is  presented  in  detail  in  the  final  section. 

HiverJ  nf  lqiws 

86.  United  States  Geological  Survey  {USGS)  daily  average  stream! lows 
were  tabulated  from  Apr i 1 -September  1980  for  the  major  streams  to  be  considered 
in  the  numerical  model  as  shown  in  Table  VII— 1 .  Flow  for  the  Biloxi  River 
System  consisted  of  flows  in  the  Tchoutacabouf fa  and  Biloxi  Rivers  and  Bayou 
Bernard.  Total  drainage  areas  for  these  three  systems  totaled  588.68  square 
miles  (Bettandorff  1972).  Daily  mean  flows  for  the  Biloxi  River  at  Wortham, 
Miss.,  (Station  02481000j  were  obtained  from  the  USGS.  The  drainage  area  at 
Wortham,  Miss.,  was  given  as  96.11  square  miles.  The  total  flow  for  the  Biloxi 
River  System  was  estimated  by  multiplying  the  Wortham,  Miss.,  gage  reading  by 
6.1251  =  588.68/96.11.  The  flow  for  the  Jourdan  and  Wolf  System  was  obtained 
similarly.  The  drainage  areas  totaled  759.20  square  miles  for  the  two  rivers 
at  Mississippi  Sound  (Bettandorff  1972).  Daily  mean  flows  at  Station  02481510 
on  the  Wolf  River  near  Landon,  Miss.,  (drainage  area  .308.28  square  miles)  were 
multiplied  by  2.4627  =  759.20/308.28  to  estimate  the  total  inflow  to  the 
Sound . 

87.  The  total  drainage  area  of  the  Pascagoula  River  System  as  given  by 
the  USGS  is  9498  square  miles.  Daily  mean  flows  at  Station  02479000  at 
Merrill,  Miss.,  gage  (drainage  area  6590  square  miles)  were  multiplied  by 
1.4413  =  9498/6590  to  estimate  the  total  inflow  into  Mississippi  Sound.  Based 
upon  USGS  Water  Supply  Paper  No.  1763  (Harvey  et  al.  1965),  this  flow  was  di¬ 
vided  between  the  West  Pascagoula  and  Pascagoula  Rivers  as  follows.  The  West 
Pascagoula  River  received  0.6  of  the  flow  with  the  remaining  0.4  ol  the  flow 
being  assigned  to  the  Pascagoula  River. 
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Table  VII-1 

Daily  Discharges  (USGS)  into  Mississippi  Sound,  in  Cfs 


Apri 1-September  1980 


River  Systems 


Pearl 

Jourdan 

_w°n 

Bi loxi 

West 

Pascagoula 

Pascagoula 

Mobile 

122,929 

8,422 

4,753 

85,476 

57,651 

369,150 

117,286 

8,472 

7,656 

93,395 

62,263 

385,200 

116,633 

13,816 

12,985 

99,448 

66,299 

354,170 

111,403 

8,644 

8 , 085 

95,124 

63,416 

372,360 

102,904 

5,566 

3,957 

85,093 

56,729 

395,900 

93,489 

4,137 

2,437 

73,678 

49,119 

396,970 

85,252 

4,014 

2,499 

64,252 

42,835 

395,900 

80,152 

4,039 

2,536 

56,988 

37,992 

383,060 

78,976 

4,014 

2,413 

49,464 

32,976 

325,280 

72,176 

4,261 

1,611 

42,114 

28,076 

300,349 

67,600 

3,817 

1,255 

35,715 

23,810 

329,560 

66,423 

9,087 

11,944 

34,936 

23,291 

315,650 

71,784 

31,769 

38,955 

45,054 

30,036 

297,353 

87,082 

31,276 

33,627 

55,172 

36,781 

239,573 

105,519 

24,405 

10,045 

62,868 

41,912 

247,384 

110,095 

7,782 

4,269 

77,051 

51,367 

247,170 

105,519 

4,531 

2,854 

94,259 

62,840 

212,609 

102,250 

3,546 

2,144 

98,583 

65,722 

209,720 

101,596 

3,325 

1,960 

91 ,665 

61,110 

225,235 

97,673 

2,980 

1,838 

80,683 

53,788 

261,080 

89,959 

2,487 

1,531 

69,959 

46,640 

278,200 

74,268 

2,066 

1,103 

60,793 

40,529 

294,250 

77,668 

1,822 

919 

49,378 

32,919 

328,490 

74,268 

1 ,625 

796 

39,693 

26,462 

339,190 

72,176 

1,478 

735 

31 ,910 

21,273 

243,470 

78,976 

4,728 

7,963 

28,796 

19,198 

347,750 

82,114 

8,373 

6,125 

32,342 

21 ,561 

342,400 

70,607 

6,009 

3,063 

34,072 

22,714 

322,070 

63,416 

3,497 

1  ,838 

33,985 

22,657 

304,950 

59,362 

2,266 

1,225 

33,120 

22,080 

271,887 

57,532 

1  ,899 

919 

32,342 

21,561 

246,742 

54,786 

1 ,675 

735 

29,748 

19,832 

234,437 

51,779 

2,194 

1 ,531 

27,067 

18,045 

218,280 

46,548 

1,877 

980 

24,040 

16,027 

205,012 

42,364 

1,428 

735 

19,890 

13,260 

186,822 

38,703 

1 ,133 

613 

17,468 

1 1 ,645 

161,784 

(Continued) 
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Table  VII-1  (Continued) 


River 

Systems 

Date 

Pearl 

Jourdan 

Wolf 

Biloxi 

West 

Pascagoula 

Pascagoula 

Mobile 

lay 

(Continued) 

7 

35,826 

1  ,034 

551 

15,912 

10,608 

160,179 

8 

33,473 

1  ,010 

521 

14,269 

9,512 

190,032 

9 

30,204 

985 

502 

12,712 

8,475 

110,531 

10 

23,366 

948 

490 

11,501 

7,668 

86,991 

1  1 

20 , 1 36 

886 

465 

10,377 

6,918 

80,143 

12 

16,867 

637 

441 

9,512 

6,342 

65,056 

13 

1 4 , 644 

788 

416 

8,648 

5,765 

68,694 

14 

12,369 

7  76 

392 

8,025 

5,350 

70,620 

13 

10,604 

788 

368 

8,189 

5,460 

65,591 

lb 

13,598 

2,857 

1,225 

9,512 

6,342 

55,961 

1  7 

28,766 

1 6  ,  +0  1 

15,313 

21,879 

14,586 

42,265 

18 

45,241 

16.254 

9,188 

36,234 

24,156 

43,763 

19 

57,793 

23,790 

27,563 

43,152 

28,768 

41,409 

20 

73,876 

36 ,694 

24,500 

52,837 

35,225 

44,619 

21 

65,769 

30,784 

12,250 

61,831 

41 ,220 

53,179 

22 

59,362 

17,608 

4,288 

61 ,571 

41,047 

121,445 

23 

53,609 

6,009 

2,756 

60,447 

40,298 

199,662 

24 

48,510 

4,113 

1,838 

57,853 

38,569 

214,749 

23 

46,287 

2,931 

1,347 

54,913 

36 , 608 

228,124 

26 

44,587 

2,327 

1 ,041 

53,269 

35,513 

263,220 

2  7 

42,364 

1  ,963 

796 

53,010 

35,340 

280,340 

28 

37,526 

1  ,6  75 

674 

49,983 

33,322 

295,320 

29 

31  ,381 

1  ,478 

551 

43,584 

29,056 

294,250 

30 

24,974 

1  , 305 

490 

37,271 

24,848 

276,060 

31 

20,921 

1  ,  157 

453 

30,786 

20,524 

278,200 

June 

1 

18,698 

1,07  l 

416 

21,879 

14,586 

273,920 

2 

16,998 

960 

367 

14,009 

9,339 

260,010 

3 

14,383 

874 

331 

10,291 

6,860 

136,318 

4 

11,716 

808 

306 

8,994 

5,995 

111,815 

3 

9,624 

739 

282 

8,189 

5,460 

95,444 

b 

7,924 

690 

263 

7,498 

4,998 

76,933 

7 

6,982 

640 

251 

6,961 

4,641 

66,233 

8 

6,904 

566 

239 

6,425 

4,283 

57,673 

9 

7,165 

591 

220 

6,399 

4,266 

52,216 

10 

6,825 

566 

233 

7,160 

4,774 

52,216 

(Continued) 
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Table  VII-1  (Continued) 


River 

Systems 

Date 

Pearl 

Jourdan 

Wolf 

Biloxi 

West 

Pascagoula 

Pascagoula 

Mobile 

June 

(Continued) 

1 1 

6,015 

566 

208 

6,797 

4,531 

50,183 

12 

5,570 

520 

190 

7,281 

4,854 

37,022 

13 

5,269 

475 

172 

7,169 

4,779 

33,919 

14 

5,086 

446 

159 

6,685 

4,456 

23,273 

15 

4,890 

419 

147 

6,053 

4,036 

23,112 

16 

4,773 

406 

135 

5,318 

3,546 

19,945 

17 

4,655 

389 

129 

4,817 

3,211 

17,302 

18 

4,576 

387 

123 

4,566 

3,044 

16,820 

19 

4,642 

419 

153 

4,479 

2,986 

18,115 

20 

4,668 

431 

135 

4,436 

2,957 

17,923 

21 

4,576 

441 

123 

4,454 

2,969 

21,956 

22 

4,472 

414 

110 

4,687 

3,125 

21,828 

23 

4,472 

458 

116 

5,085 

3,390 

23,925 

24 

4,655 

677 

214 

4,869 

3,246 

27,360 

25 

4,838 

857 

184 

4,211 

2,808 

25,894 

26 

4,579 

823 

153 

4,635 

3,090 

30,281 

27 

4,576 

754 

123 

5,422 

3,615 

37,236 

28 

4,943 

549 

110 

5,967 

3,978 

24,396 

29 

5,871 

468 

98 

6,884 

4,589 

29,853 

30 

7,924 

468 

92 

6,754 

4,503 

63,986 

July 

1 

9,022 

539 

123 

5,578 

3,718 

73,509 

2 

9,218 

529 

98 

4,817 

3,211 

57,994 

3 

9,676 

396 

86 

4,255 

2,836 

37,450 

4 

10,434 

342 

73 

3,943 

2,629 

29,746 

5 

11,519 

313 

67 

3,701 

2,467 

21,775 

6 

12,212 

296 

61 

3,494 

2,329 

22,342 

7 

11,964 

278 

59 

3,364 

2,243 

28,783 

8 

10,447 

276 

56 

3,563 

2,375 

29,960 

9 

7,963 

320 

86 

4,099 

2,733 

23,439 

10 

6,642 

362 

80 

3,753 

2,502 

21,689 

1  1 

5,975 

335 

67 

3,442 

2,295 

18,768 

12 

5,243 

288 

61 

3,191 

2,127 

19,314 

13 

4,603 

264 

58 

3,035 

2,024 

23,198 

14 

4,250 

244 

53 

2,888 

1,926 

26,108 

15 

4,067 

236 

51 

2,759 

1,839 

26,483 

16 

3,923 

229 

50 

2,672 

1,781 

27,595 

(Continued) 
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Table  VI I  —  1  (Continued) 


River 

Systems 

Date 

Pearl 

Jourdan 

Wolf 

Biloxi 

West 

Pascagoula 

Pascagoula 

Mobile 

July 

(Continued) 

17 

3,949 

229 

50 

2,629 

1,753 

18,821 

18 

3,923 

288 

51 

2,689 

1,793 

18,211 

19 

4,001 

328 

61 

2,854 

1,902 

16,221 

20 

4,406 

406 

80 

3,277 

2,185 

18,425 

21 

5,034 

1,953 

484 

3,727 

2,485 

16,275 

22 

5,544 

2,931 

1,323 

5,906 

3,938 

15,857 

23 

6,015 

1,832 

1,023 

7,212 

4,808 

17,441 

24 

6,629 

1,271 

606 

7,722 

5,148 

15,665 

25 

7,061 

798 

331 

7,964 

5,310 

15,697 

26 

7,845 

589 

208 

6,711 

4,474 

14,541 

27 

r^. 

oc 

r-» 

475 

184 

5,768 

3,845 

15,943 

28 

7,715 

426 

178 

5,223 

3,482 

17,559 

29 

8,028 

589 

196 

5,748 

3,165 

19,442 

30 

6,211 

525 

159 

4 , 661 

3,107 

17,366 

31 

5,479 

480 

141 

4,782 

3,188 

21,817 

August 

1 

5,243 

406 

110 

4,514 

3,009 

23,208 

2 

4,851 

340 

98 

3,935 

2,623 

21,871 

3 

4,406 

305 

92 

3,468 

2,312 

23,273 

4 

4,119 

313 

129 

3,200 

2,133 

21,421 

5 

3,962 

414 

172 

3,070 

2,047 

24,781 

6 

3,844 

335 

123 

3,234 

2,156 

26,140 

7 

3,753 

325 

92 

2,958 

1 ,972 

16,628 

8 

3,635 

278 

86 

2,733 

1,822 

16,425 

9 

3,687 

268 

80 

2,560 

1,706 

14,167 

10 

3,687 

241 

73 

2,404 

1 ,603 

14,370 

1 1 

3,596 

236 

67 

2,776 

1,851 

18,693 

12 

3,517 

251 

92 

2,724 

1,816 

17,527 

13 

3,465 

310 

135 

2,473 

1  ,649 

13,439 

14 

3,347 

340 

220 

2,352 

1  ,568 

14,830 

15 

3,282 

362 

263 

2,447 

1  ,632 

15,215 

16 

3,269 

394 

245 

2,292 

1  ,528 

12,594 

17 

3,334 

397 

190 

2,153 

1  ,436 

12,487 

18 

3,308 

328 

129 

2,058 

1,372 

14,049 

19 

3,217 

278 

153 

1  ,946 

1  ,297 

12 ,b69 

20 

3,138 

239 

104 

1,851 

1  ,234 

13,568 

21 

3,073 

256 

86 

1  ,807 

1  ,205 

13, 150 

( Cont inued ) 
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Table  VII-1  (Concluded) 


88.  The  development  of  the  Mobile  River  System  was  treated  as  described 
by  Schroeder  (1979).  The  total  inflow  to  Mobile  Bay  was  assumed  to  be  equal 
to  1.07  times  the  sum  of  the  gage  readings  for  the  Tombigbee  River  at  Coffee- 
ville,  Ala.,  Station  02429761  and  the  Alabama  River  at  Clairborne,  Ala.,  Sta¬ 
tion  0 2429500.  To  allow  for  the  travel  time,  the  resulting  inflow  time  series 
was  lagged  by  5  days  to  estimate  the  inflow  time  series  at  Mobile  Bay. 

89.  The  total  drainage  area  reported  by  the  USGS  for  the  Pearl  River  at 
its  mouth  is  8669  square  miles.  To  estimate  flow  into  Mississippi  Sound, 
daily  mean  values  at  Station  02489500  near  Bogalusa,  La.  (drainage  area  6630 
square  miles)  were  multiplied  by  1.3075  =  8669/6630. 

90.  River  inflows  were  high  in  the  spring  months  (April  and  May)  and 
receded  during  June  and  were  relatively  low  in  July,  August,  and  September. 
Total  average  daily  inflows  for  all  six  river  systems  were  518,000  cfs, 

258,000  cfs,  77,500  cfs,  39,150  cfs,  23,900  cfs,  and  16,900  cfs  for  April,  May, 
June,  July,  August,  and  September,  respectively. 


Meteorological  Station  Wind  Data 


91.  Daily  maximum  hourly  averaged  wind  speeds  and  daily  maximum  (2-sec 
gust)  wind  speeds  were  tabulated  over  the  survey  period  at  Station  4  on  Horn 
Island  as  shown  in  Table  VII-2.  For  periods  in  which  data  at  Station  4  were 
not  available,  data  were  tabulated  for  Station  2  at  Ship  Island  as  shown  in 
Table  VII-3  tc  provide  a  record  of  wind  information  as  complete  as  possible. 

The  range  of  daily  maximum  hourly  averaged  wind  speed  was  5.6  (July  10)  to  31.7 
(April  13)  mph.  The  range  of  daily  maximum  winds  speed  was  10.7  (July  10)  to 
45.4  (April  13)  mph.  Typical  daily  maximum  hourly  averaged  windspeeds  were  20, 
12,  10,  15,  13,  8-28,  and  8-24  mph  for  April,  May,  June,  July,  August,  Sep¬ 
tember,  and  October  1980.  April  and  early  May  constituted  a  period  of  rela¬ 
tively  high  winds.  Fall  storms  occurred  in  September  and  October.  The  gages 
were  removed  during  Hurricane  Allen  (7  August-17  August). 


Sal i nity  Transect  Data 


92.  Throughout  the  180-day  survey  period  at  approximately  3-week  inter¬ 
vals,  salinity  transects  were  obtained  at  locations  shown  in  Figure  1 1 - 4 .  Sa¬ 
linity  values  were  measured  at  5-ft-depth  intervals.  In  order  to  characterize 
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Table  VII-2 


Wind  Characteristics  for  Meteorological  Station  4,  Horn  Island 


Maximum 

Maximum 

Hourly 

Maximum 

Hourly 

Maximum 

Julian 

Average  Speed 

Gust 

.Julian  Average  Speed  Gust 

Date 

Day 

mph 

mph 

Date 

Da>'_ 

mph 

mph 

May  8 

129 

12.6 

20.0 

16 

168 

7.6 

15.4 

9 

130 

24.0 

31.4 

17 

169 

8.3 

16.0 

10 

131 

10.6 

18.7 

18 

170 

11.4 

21.4 

1  1 

132 

11.5 

22.0 

19 

171 

9.9 

22  .  7 

12 

133 

10.8 

18.7 

20 

172 

10.8 

16.0 

13 

134 

10.8 

20.0 

21 

173 

10.1 

18.7 

14 

135 

9.9 

18.0 

22 

174 

10.0 

20. 7 

15 

136 

18.3 

38.1 

23 

175 

10.0 

25.4 

16 

137 

31.9 

43.4 

24 

176 

10.0 

32.7 

17 

138 

26.7 

38.7 

25 

177 

11.3 

23.4 

18 

139 

10.0 

16.0 

26 

178 

15.0 

36.1 

19 

140 

13.4 

26.7 

27 

179 

9.2 

21.4 

20 

141 

17.6 

31.4 

28 

180 

9.3 

16.7 

21 

142 

11.9 

29.4 

29 

181 

9.6 

18.7 

22 

143 

14.0 

21.4 

30 

182 

13.4 

21 .4 

23 

144 

9.5 

17.4 

Jul  1 

183 

14.3 

28.7 

24 

145 

9.8 

18.0 

2 

184 

11.3 

19.4 

25 

146 

10.2 

18.0 

3 

185 

11.2 

20.7 

26 

147 

10.3 

34.7 

4 

186 

8.3 

13.4 

27 

148 

12.1 

17.4 

5 

187 

8.7 

15.4 

28 

149 

9.7 

16.7 

6 

188 

8.8 

15.4 

29 

150 

13.2 

21 .4 

7 

189 

21  .  1 

46.8 

30 

151 

9.6 

17.4 

8 

190 

19.2 

34.7 

31 

152 

8.8 

15.4 

9 

191 

7.6 

14.8 

10 

192 

5.6 

10.7 

Jun  1 

153 

8.2 

14.0 

Out 

of  Service 

2 

154 

9. 7 

17.4 

31 

213 

8.0 

14.0 

3 

155 

6.  1 

12.0 

4 

156 

6.2 

11.4 

Aug  1 

214 

7.7 

14.0 

5 

157 

13.4 

19.4 

2 

215 

8.8 

17.4 

6 

158 

6.8 

15.4 

3 

216 

9.  1 

28. 1 

7 

159 

9.2 

18.7 

4 

217 

8.4 

16.7 

8 

160 

11.3 

19.4 

5 

218 

8.3 

15.4 

9 

161 

20.9 

29.4 

6 

219 

13.4 

18.7 

10 

162 

17.6 

22.7 

Out  of 

Servi ce 

(Hurricane 

A1 1 en) 

18 

231 

6.3 

12.0 

1  1 

163 

11.3 

14.7 

19 

232 

7.9 

12.0 

12 

164 

16.  1 

20.7 

20 

233 

8.6 

12.7 

13 

165 

11.4 

16.7 

21 

234 

8.8 

20.7 

14 

1 66 

7  .  1 

13.4 

22 

235 

12.3 

22.7 

15 

167 

7.2 

13.4 

(Continued) 
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Table  VI [-2  (Concluded) 


Maximum 

Maximum 

Hourly 

Maximum 

Hourly 

Maximum 

Ju  1  i  an 

Average  Speed 

Gust 

Julian 

Average  Speed 

Gust 

Date 

Day 

mph 

mph 

Date 

Day 

mph 

mph 

Aug  23 

236 

19  .  1 

25.4 

0(  t  1 

275 

18.3 

24.0 

24 

237 

18.4 

23.4 

2 

276 

10.4 

14.0 

25 

238 

10.5 

20.0 

!  3 

277 

22.6 

30. 1 

26 

239 

12.7 

16.7 

-♦ 

278 

19.  1 

26.  1 

27 

240 

13.3 

18.0 

s 

2  79 

9.6 

29 . 4 

28 

241 

13.8 

22.0 

0 

280 

23.7 

30.  1 

29 

242 

8.8 

16.7 

7 

281 

17.2 

21.4 

30 

243 

7.4 

1 0 . 0 

8 

282 

8.8 

11.4 

11 

244 

8.  ' 

16.7 

9 

283 

6.2 

11.4 

10 

284 

10.9 

1  4 . 0 

Sep  1 

245 

1  0 . 4 

io.O 

246 

1  7  .  i 

25.4 

]  1 

285 

12.7 

19.4 

3 

247 

24.0 

12.  1 

1  2 

286 

21  .  1 

27.4 

4 

248 

28.0 

38 .  1 

13 

287 

17.0 

21.4 

5 

249 

20.0 

32.7 

14 

288 

9.9 

12.7 

15 

289 

15.1 

19.4 

6 

250 

12.0 

18. 7 

7 

251 

12.2 

15.4 

16 

290 

15.1 

26.7 

8 

252 

7.9 

17.7 

17 

291 

10.2 

18.0 

9 

253 

8.3 

14.0 

18 

292 

7.9 

15.4 

10 

254 

8.6 

14.7 

19 

293 

18.5 

22.7 

20 

294 

24.0 

32.7 

1  1 

255 

14.5 

18.0 

12 

256 

10.0 

18.7 

21 

295 

18.8 

24.7 

13 

257 

9.6 

14.7 

22 

296 

II  .5 

14.7 

14 

258 

12.3 

18.7 

23 

29  7 

13.3 

20.0 

15 

259 

8.  1 

14.0 

24 

298 

20.  7 

34.  1 

25 

299 

23.6 

3b.  I 

16 

260 

8.0 

15.4 

1  7 

261 

10.3 

32.  7 

26 

300 

11.8 

18.0 

18 

262 

15.6 

18.4 

27 

301 

15.8 

28.  1 

19 

263 

1.3.3 

23.4 

28 

30? 

17.0 

31.4 

20 

264 

8.9 

1<1. 4 

29 

23.7 

30. 7 

30 

.>u4 

21.7 

28  .  1 

21 

265 

8.8 

15.4 

31 

305 

20.8 

26.  1 

22 

266 

7 . 2 

13.4 

23 

267 

9.  1 

14.7 

Nov  1 

.306 

11.2 

14.7 

24 

268 

6 . 6 

lb.  7 

2vV 

307 

7.2 

10.7 

25 

269 

7.5 

15.4 

3* 

308 

7.3 

15.4 

26 

270 

17.2 

22.0 

4* 

309 

19.5 

20.  7 

27 

271 

19.3 

24. 7 

5'-' 

3  1 0 

15.8 

31.4 

28 

272 

17.  1 

22.0 

6* 

311 

12.5 

16.0 

29 

273 

9.8 

19.4 

7* 

312 

7.7 

20 . 0 

30 

274 

20.2 

32.7 

Gage  possibly  ma  1  furirt  ion  i  ng  . 
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Table  VII-3 

Wind  Characteristics  for  Meteorological  Station  2,  Ship  Island 


Maximum 

Maximum 

Hourly 

Maximum 

Hourly 

Maximum 

Jul  lan 

Average  Speed 

Gust 

Julian 

Average  Speed 

Gust 

Date 

Day 

mph 

mph  i 

Date 

Day 

mph 

mph 

Apr  10 

101 

13.  1 

i 

! 

16. 7 

Jul  9 

191 

10.3 

14.0 

1  1 

102 

21.9 

28. 1 

10 

192 

10.5 

14.0 

12 

103 

30.5 

44.8 

1  I 

193 

14.3 

18.7 

13 

104 

31.7 

45.4  | 

12 

194 

19.5 

24.0 

14 

105 

29.5 

40.  1 

1 

13 

195 

18.7 

24.0 

15 

106 

23.  1 

30.7  i 

14 

196 

18.8 

22.7 

16 

107 

10.3 

14.0 

15 

197 

16.6 

20.7 

17 

108 

10.4 

14.7 

16 

198 

14.9 

21.4 

18 

109 

18.4 

34 . 7 

17 

199 

15.1 

23.4 

19 

1 10 

26.6 

37.4 

18 

200 

16.7 

38.1 

20 

111 

13.0 

16.7  ! 

19 

201 

12.5 

28.7 

21 

1 12 

14.7 

20.0 

20 

202 

22.3 

28.1 

22 

113 

18.0 

24.7 

21 

203 

20.3 

27.4 

23 

114 

16. 1 

20.7 

22 

204 

15.9 

24.7 

24 

115 

17.6 

22.7 

23 

205 

16.6 

33.4 

25 

116 

20.8 

28.7 

24" 

206 

23.6 

26.1 

25" 

207 

20.5 

23.4 

26 

117 

25.8 

44. 1 

27 

1  18 

24.  1 

32.1 

26" 

208 

25.3 

24.7 

28 

119 

26.0 

28.1 

27" 

209 

28.9 

30.7 

29 

120 

21.5 

28.1 

28" 

210 

29.6 

31.4 

30 

121 

18.9 

24.0 

29" 

211 

34.  1 

17.4 

30 

212 

14.8 

18.7 

May  1 

122 

14.8 

19.4 

31 

213 

12.4 

18.7 

2 

123 

15.8 

20.0 

3 

124 

15.8 

28.7 

Aug  1 

214 

11.6 

15.1 

4 

125 

15.5 

20.7 

2 

215 

17.1 

30.1 

5 

126 

14.8 

18.7 

3 

216 

21.7 

36.7 

4 

217 

15.8 

20.0 

6 

127 

18.0 

24.0 

5 

218 

17.3 

22.0 

7 

128 

18.5 

24.0 

6 

219 

14-3 

20.0 

8 

129 

17.7 

29.4 

Out 

of  Service 

Gage  possibly  malfunctioning. 


the  general  horizontal  salinity  distribution,  middepth  salinity  conditions  for 
11  stations  are  shown  in  Table  V1I-4  with  the  range  over  depths  listed  under¬ 
neath  the  middepth  value.  Stations  T26 ,  T28,  T30,  T36 ,  and  T42  characterize 
the  western  portion  of  the  Sound.  Stations  T2,  T52,  and  T54  represent  the 
central  portion  of  the  Sound,  while  stations  T60,  T70,  and  T80  portray  condi¬ 
tions  in  the  eastern  section  of  the  Sound.  In  referencing  the  table,  if  one 
goes  across  the  table  at  any  one  date,  the  horizontal  salinity  distribution  is 
obtained.  If  one  goes  down  a  column,  the  change  in  conditions  over  time  are 
represented  for  each  station,  group,  or  for  the  entire  Sound.  In  April  and 
May  1980,  the  entire  Sound,  particularly  the  western  end  and  the  entrance  to 
Mobile  Bay,  exhibited  very  low  salinity.  As  June  progressed,  salinity  levels 
increased  overall.  This  trend  continued  through  July.  By  August,  conditions 
in  the  Sound  had  stabilized  to  a  summer  pattern.  Note  that  the  8/22-23  tran¬ 
sect  values  correspond  very  closely  to  the  9/24-25  values. 

93.  Horizontal  salinity  gradients  are  largest  in  the  Spring.  However, 
local  gradients  may  be  found  throughout  the  period  in  the  highly  transient 
Lake  Borgne  (T26)  and  Mobile  Bay  (T80)  entrances. 

Cal ib ration  an d  Ver if i cation  Periods 

94.  It  is  desirable  to  calibrate  and  verify  the  hydrodynamic  parameters 
(bottom  stress)  and  salinity  parameters  (dispersion  coefficients)  separately 
but  over  the  same  periods.  The  following  criteria  are  used  to  select  the  cali 
bration  and  verification  periods: 

a.  Streamflow  must  remain  relatively  constant.  Storm  or  flood 
periods  exhibiting  time-varying  flow  characteristics  are  dis¬ 
carded  from  consideration. 

b.  In  order  to  minimize  wind  effects,  periods  for  which  wind 
speeds  are  minimum  are  to  be  considered. 

e.  Since  continuous  salinity  levels  are  available  at  only  six  sta¬ 
tions  (VS  stations  in  Figure  1 1  —  1 ) ,  it  will  be  necessary  to  em¬ 
ploy  salinity  transect  values  in  defining  the  initial  condi¬ 
tions.  Therefore,  these  periods  must  begin  in  a  transect  day. 

95.  By  considering  river  inflows  (Table  VII-1),  wind  conditions 
(Tables  VII-2  and  VI 1-3),  and  salinity  transect  conditions  in  Table  V 1 1  — 4 ,  a 
general  understanding  of  the  dynamics  of  the  Sound  may  be  obtained. 

96.  In  the  numerical  study  of  the  Sound,  it  is  necessary  (for  rea¬ 
sons  of  computational  cost)  to  consider  a  5-  to  6-day  period.  This  constraint 
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Table  VI 1-4 

>  resent  at  Ive  Spatial  Sa  1  inity  Tran  sort  Pat  a  ,  1  ‘jh  0  (  p  p  t  ) 


coupled  with  the  previous  three  criteria  is  sufficient  to  determine  the  pe¬ 
riods  for  numerical  study. 

97.  Consider  Table  VII-4;  the  4/28-4/29  period  is  eliminated  due  to 
excessive  wind  and  highly  transient  flood  flow  conditions.  The  S /2 1 —5/22  pe¬ 
riod,  although  relatively  calm,  exhibits  transient  flow  conditions  and  is  dis¬ 
carded.  The  7/23-24,  8/22-23,  and  9/8-9  periods  exhibit  excessive  wind  and 
are  eliminated.  The  9/2-3  transect  is  incomplete  and  therefore  this  period  is 
not  considered. 

98.  The  6-day  period,  9/20-9/25,  was  selected  for  further  study  as  a 

potential  calibration  period  for  both  hydrodynamics  and  salinity  mechanisms. 
From  the  hydrodynamic  perspective,  the  astronomic  tide  characteristics  at 
Pascagoula  as  shown  by  Outlaw  (1983,  Plate  42)  indicate  that  during  9/23  and 
9/24  the  semidiurnal  constituents  M^  and  assume  increasing  importance 

as  the  tidal  range  is  reduced  on  entering  a  neap  tide  period.  From  the 
salinity  perspective,  values  are  relatively  constant  in  time  over  this  period. 
Therefore,  salinity  transect  values  obtained  over  the  2-day  period,  20-21  Sep¬ 
tember  may  be  used  to  define  initial  conditions.  In  addition,  the  24-25  Sep¬ 
tember  values  may  be  used  as  values  to  compare  with  the  numerical  results. 
Maximum  daily  hourly  averaged  wind  speeds  are  less  than  9.0  mph  over  the  en¬ 
tire  period.  Therefore,  wind  effects  should  be  small. 

99.  The  5-day  period,  6/12-6/16,  was  selected  for  further  study  as  a 
potential  verification  period  for  both  hydrodynamic  and  salinity  mechanisms. 
From  the  hydrodynamic  viewpoint,  the  astronomic  tide  characteristics  at 
Pascagoula  as  shown  in  Outlaw  (1983,  Plate  39)  indicate  that  the  diurnal  con¬ 
stituents  0^  ,  Kj  ,  Pj  dominate  and  determine  the  character  of  the  tide. 
The  tide  range  is  characteristic  of  a  Spring  tide.  From  the  salinity  view¬ 
point,  this  period  exhibits  a  relatively  large  horizontal  gradient.  Vertical 
salinity  gradients  are  significant  in  the  middle  and  eastern  station  groups  in 
Table  VII- 4.  The  maximum  daily  hourly  averaged  wind  speed  declines  from  16.1 
to  11.4  mph  from  6/12-6/13  and  to  less  than  8.0  mph  for  the  6/14-6/16  period. 
The  vertical  salinity  gradients  indicate  that  the  middle  and  eastern  sections 
of  Mississippi  Sound  are  partially  mixed,  while  the  western  section  is  rea¬ 
sonably  well  mixed.  The  vertically  integrated  two-dimensional  modeling  con¬ 
cept  employed  for  salinity  may  not  be  applicable  during  this  period.  As  a  re¬ 
sult,  salinity  was  not  considered.  The  hydrodynamic  mechanisms  were  verified 
using  this  Spring  tide  period. 


PART  VIII:  CALIBRATION  PERIOD  20-25  SEPTEMBER  1980 


100.  In  order  to  obtain  a  more  detailed  understanding  of  the  dynamics 
of  the  Sound  prior  to  simulating  this  period,  hourly  values  of  wind  speed  and 
direction  and  salinity  are  considered  in  turn  below.  A  complete  set  of  sa¬ 
linity  transect  values  are  presented  in  order  to  investigate  the  degree  of 
stratification.  Based  upon  the  results  of  the  harmonic  analysis,  the  pre¬ 
dicted  values  of  water  surface  elevation  and  currents  are  developed  and  com¬ 
pared  versus  the  unfiltered  (raw)  data. 

Wind  Information 

101.  Hourly  wind  speeds  and  directions  are  as  shown  in  Table  VI 1 1  —  1 . 

At  most  stations  average  hourly  wind  speeds  are  below  10  mph,  with  daily  maxi¬ 
mum  wind  speeds  (2-sec  gust)  below  20  mph.  Wind  direction  and  velocity  are 
relatively  uniform  spatially  at  each  hour  over  the  period.  The  wind  effects 
should  be  relatively  small  on  water  surface  displacements  over  this  period. 

Instantaneous  Salinity  Information 

102.  Hourly  values  of  salinity  meters  nearest  middepth  were  considered. 
Middepth  locations  were  selected  as  being  the  most  appropriate  to  compare  with 
the  vertical  averaged  model  results.  The  range  for  each  day  at  each  station 
as  located  in  Figure  1 1  —  1  is  given  in  Table  VIII-2.  Salinity  levels  are  nor¬ 
mally  within  a  range  of  approximately  1-2  parts  per  thousand  for  all  days. 

The  maximum  spatial  difference  over  the  Sound  is  approximately  5-6  parts  per 
thousand . 


Sa 1 inity  Transect  Information 

103.  Representative  middepth  salinity  values  are  shown  in  Figure  VI 1 1 - 1 . 
The  salinity  range  over  depth  is  given  to  indicate  the  degree  of  stratifica¬ 
tion.  Stratification  effects  even  within  the  navigation  channels  were  not 
significant  during  this  period.  Therefore,  the  well-mixed  assumption  appeared 
to  be  valid  during  this  period. 

104.  In  order  to  investigate  the  reliability  of  the  instantaneous 
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Table  VI 1 1  —  1  (Continued) 


Next  larger  value. 


Table  VI 1 1 - 1  (Continued) 


(Continued) 


24-hr  24-hr  24-hr  24-h 

ate/  Direc-  Max  Di rec-  Max  Direc-  Max  Direc-  Max 
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Station  15B  was  used. 

Some  spurious  data  points  were  discovered. 


September  20-21  salinity  transect  da 


salinity  and  temperature  data  obtained  using  the  deployed  conductivity- 
temperature  sensors,  the  salinity  and  temperature  transect  values  obtained 
using  a  field  salinometer  were  compared  at  common  times  and  spatial  locations 
as  shown  in  Table  VIII-3. 

105.  In  most  instances,  temperature  values  were  extremely  close.  How¬ 
ever,  the  salinity  values  obtained  from  the  conductivity-temperature  sensors 
tended  to  be  2-3  parts  per  thousand  lower  than  the  field  salinometer  values. 
Raytheon  Ocean  Systems  Company  (Parker  1981)  indicated  that  there  were  bio¬ 
logical  fouling  problems  with  the  conductivity  probe.  Antifouling  paints  can¬ 
not  be  used  on  the  sensor  because  the  metal  ions  released  by  such  paints  affect 
the  conductivity  readings.  Therefore,  in  analyzing  the  salinity  data,  the  in¬ 
stantaneous  salinity  data  should  be  considered  less  reliable  then  the  salinity 
transect  data.  Since  salinity  transect  data  were  available  for  20-21  September 
and  24-25  September  these  data  were  used  to  establish  initial  conditions  and 

as  calibration  data,  respectively. 

Water  Surface  Elevation  Information 

106.  Program  TIDE  (Appendix  B)  was  developed  to  predict  the  tidal  eleva¬ 
tions  presented  in  Part  III  after  Schureman  (1940)  based  upon  the  constituent 
amplitudes  and  phases.  Astronomic  data  are  calculated  in  order  to  determine 
constituent  node  factors  and  equilibrium  arguments  at  the  start  of  the  pre¬ 
diction  period.  Based  upon  specified  amplitudes  and  local  epochs  for  the 

0^  ,  ,  Pj  ,  M^  ,  anc*  "*2  constituents >  as  given  previously  in 

Tables  1II-2  and  1 1 1  —  3 ,  the  program  reconstructs  the  water  surface  elevation 
for  each  station.  Subroutine  TAPE  accesses  the  edited  hourly  elevations  as 
provided  by  NOAA  and  the  filtered  elevations  input  to  the  harmonic  analysis  so 
that  these  elevations  may  be  printed  next  to  the  reconstructed  elevation  at 
each  hour  of  the  prediction  period. 

107.  Table  VIII-4  presents  the  astronomic  data  with  reference  to  Green¬ 
wich  Mean  Time  (GMT).  By  employing  local  epochs,  the  time  reference  is  trans¬ 
formed  directly  to  Central  Standard  Time  (CST).  Consider  the  Pascagoula  sta¬ 
tion,  874-1196;  tidal  constituents  are  given  in  Table  VI 1 1 —5 .  Since  local 
epochs  are  used,  the  water  surface  elevations  given  in  Table  VI 1 1-6  are  in  CST 
beginning  at  20  September,  hour  00.  Elevations  are  given  for  the  first  2  days 
of  the  5-day  period.  The  unfiltered  (UF)  and  filtered  (F)  columns  present  the 
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Tahir  VI ! 1 -3 


Comparison  of  20-21  September  Salinity  Tr.insm  t  Hat  a  with 
Instantaneous  Salinity  Data 


Location  Watt'r  Measurement 


Slat  ton 

l)ate-Jul 

i  .in  Day 

CST 

Time 

Lat i tude 

_ deg 

Long  i  tmie 
deg 

Depth 
f  t 

Dept  h 
ft 

S.i  ]  nut  y 

ppt 

Tempe  ra l u  re 
°C 

T  42 

20 
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Table  VI [1-4 

Verification  Period,  Astronomic  Tide  Conditions, 
2 0  Sept  era b er  1980,  Hour  0000 


As  t  ronom  ic 


T crms  at  Greenwich : 
Term 


Longitude  of  lunar  node 

Mean  longitude  of  sun 

Longitude  of  solar  perigee 

Mean  longitude  of  moon 

Longitude  of  lunar  perigee 

Intersection  of  lunar  orbit  and  equator 

Longitude  in  the  lunar  orbit 

Longitude  in  the  celestial  orbit 

V  ''rime 

V  double  prime 
Solar  hour  angle 


Degrees 

L37. 96955 
179. 044; u 
>82.60883 
303.55682 
18. 7900 
19.91312 
9.44233 
10.1 c2°9 
h  .  >,'/  _ 00 
1  2 . 8()6  1 ) 
1  80.00000 


Node  Factors  and 

Tidal 
Con- 
st ituent 

v: 


N  od  e 
Factor 

1  .  (2  f 


Equilibrium  Argument 

(V  +  U) 


1.668 
.tee 
l .  T 

i  -  c  :■  ; 

l  • '  f 
l  - 1  9  7 
.735 
.“7.; 

.  6  T-; 
l .  •.  o  t 
l . '  f.  p 
.  6  7  > 

l .  r  r  c 

1  •  t  2  b 
1  ■ 4  2  “ 

•  C  t  • 

.“73 

“TTIT" 

1.117 
1  .  1 5  7 


127.3? 
?<=  ’.01 


329.32 
183.66 
26 ft.77 
353.t  r 

r'-.n 

24?.. 2 
228.23 
246.4  ^ 

2  e c . 3  : 

1  .3.56 
7b. 4ft 
1.13 
9r.9f 
_  2 5  .'j.4_^_ 

164.33 
21 “.21 
316.7ft 

“3^755“ 

76.22 

l'>9.55 


Components 
V _ U_ 

9. 98  358.58 


8.73 

355.74 


357.16 


C."0 

358.53 


33  .-‘6 
35°. 58 


3.64 

349.35 


156.16 

255.23 


325.48 

1  a3 . 81 

3  c  5 . 0  ° 
17<3  .  04 

~2*7?7T? 

247.11 

237.67 

T7TT1? 

1  23.  *6 
76.44 

-55?7T!' 
90.96 
249_.  C2 
~T  66.46 

2  lc  .  7ft 
312.  °1 

~55V.1« 
93.99 
11C. 9« 


341.12 

8.73 


354  .4  7 
3.-3 


354.31 

358.58 


Table  VIII-5 

Tidal  Constituents  at  Pascagoula  (8 74- 1196) 


NGVD  -  +0.58  Ft  Above  MSL 
Amjj  Jj_t  ude  ,  f  l 


L  oca  LSr  ol  i 


unfiltered  (edited  NOAA)  and  filtered  (harmonic  analysis  input)  data,  re¬ 
spectively.  All  elevations  are  with  respect  to  NGVD  (1929).  Information 
analogous  to  that  shown  in  Table  VIII-6  for  station  874-1196  is  developed  by 
the  program  for  all  other  stations  shown  in  Table  I  1-3. 

108.  In  interpreting  the  elevation  table,  it  should  be  noted  that  pre¬ 
dicted  levels  are  based  upon  the  five  major  tidal  constituents  < 0 ^  , 

Pj  ,  ,  and  S^),  whereas  the  unfiltered  levels  contain  all  tidal  energy. 

Reid  and  Whitaker  (1981)  note  that  these  five  tidal  constituents  contain  ap¬ 
proximately  95  percent  of  the  tidal  energy. 

109.  The  unfiltered  levels  contain  meteorological  input  in  addition  to 
all  the  tidal  energy.  In  examining  Table  VIII-6,  the  filtered  and  recon¬ 
structed  levels  are  usually  very  close  (less  than  0.1  ft).  Unfiltered  levels 
are  0.4-0. 6  ft  higher  at  each  hour.  The  winds  shown  in  Table  VIII-1  do  not 
appear  to  be  of  sufficient  strength  to  account  for  all  of  the  difference.  In 
addition,  this  pattern  in  elevations  holds  at  all  other  tidal  stations  as  well 
The  majority  of  the  0.4-  to  0.6-ft  difference  between  unfiltered  and  recon¬ 
structed  levels  is  believed  to  constitute  a  measure  of  the  forerunner  and  far 
field  surge  of  tropical  storm  Hermine,  which  passed  through  the  Carribean  over 
the  Yucatan  Peninsula  and  through  the  lower  Gulf  of  Mexico  during  20-25  Sep¬ 
tember  as  shown  in  the  1980  supplement  to  Cry  (1977).  The  magnitude  of  this 
effect  as  noted  from  the  data  is  less  than  0.6  ft,  which  is  within  the  range 
of  magnitudes  of  these  effects  observed  for  other  tropical  storms  and 
hurricanes . 


Current  Component  Information 

110.  Program  TIDE  may  also  be  used  to  reconstruct  tidal  currents  as 
well.  Tidal  currents  are  represented  in  E-W  and  N-S  component  form,  with  the 
E  and  N  directions  considered  positive.  The  same  astronomic  data  as  given  in 
Table  VI 11-4  are  employed.  The  time  reference  is  again  CST  since  local  epochs 
as  shown  in  Table  VI  1 1-7  for  station  V10-S  are  used  for  each  current  component 

111.  The  unfiltered  (edited),  filtered  (harmonic  analysis  input),  and 
reconstructed  currents  for  the  first  day  of  the  5-day  period  are  given  for  sta 
tion  V10-S  in  Table  VIII-8.  Since  a  vertically  integrated  two-dimensional 
modeling  approach  is  being  employed,  the  meter  nearest  middepth  was  considered 
for  stations  shown  in  Table  III-5,  for  which  meters  were  placed  at  various 
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depths.  The  subsets  of  stations  considered  are  indicated  by  asterisks  in 
Table  III-5. 

112.  Predicted  currents  are  reconstructed  from  the  five  major  tidal  con¬ 
stituents  (Oj  ,  Kj  ,  Pj  ,  ^2  ’  ant*  ^2^’  ^1  *n  contains 

both  Kj  and  P^  energy.  In  order  to  separate  and  P^  ,  at  least 

182  days  of  record  are  necessary.  Since  current  meters  were  removed  during 
the  shrimp  season  and  during  Hurricane  Allen,  a  complete  18^-day  record  was 
not  aval lable  for  any  of  the  current  stations.  As  a  result  the  con¬ 

stituent  reported  in  Tables  I1I-6  through  1 1 1 - 9  contains  P^  energy  and  is 
valid  only  during  the  analysis  period.  Therefore,  the  total  tidal  current 
signal  may  be  predicted  at  a  given  station  only  during  the  period  given  in 
Table  MI-5. 
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PART  IX:  VERIFICATION  PERIOD  12-16  IL'NK  1980 


113.  I  n  order  to  obtain  ,i  more  detailed  mule  rs  t  .1  ml  1  ng  of  the  behavior 
ot  the  Sound  prior  to  s  i  11111 1 .1 1  1  ng  this  per  toil,  hourly  values  of  wind  speed  and 
direr  tint)  and  salinity  are  lonsidered.  A  complete  set  of  salinity  transect 
values  is  also  presented  in  order  to  investigate  str.it  1  I  ic.it  ini  effeits  Sound- 
wide.  Based  upon  the  results  of  the  harmonic  analysis,  the  predicted  values 
ot  water  surface  elevation  and  currents  are  developed  and  compart’d  with  the 
Lin  filtered  (  raw  ;  data. 


Wind  Information 

’.  i  a  .  „  ,  mi  peeds  and  directions  are  as  shown  in  Table  1 X  -  1  .  At 

most  stations  average  hourly  wind  speeds  are  below  I  f>  mph ,  while  daily  maximum 
wind  speeds  (2-sec  gust)  are  below  20  mph.  Wind  direction  and  velocity’  are 
spatially  uniform  at  any  one  hour  over  the  period. 

Instantaneous  Sa 1  imi ty  I nformat ion 

lib.  Hourly  values  of  salinity  at  meters  nearest  middepth  appear  to  ex¬ 
hibit  tidal  variations  at  most  stations.  In  order  to  investigate  the  temporal 
variations,  the  maximum  and  minimum  values  are  given  in  Table  I X - 2 ,  such  that 
their  difference  (W)  is  maximum  for  the  day  given.  Maximum  daily  variations 
in  salinity  range  from  4.18  ppt  at  station  6S  to  11.0  ppt  at  station  20M.  The 
maximum  spatial  difference  in  salinity  over  the  entire  Sound  is  greater  than 
13  ppt  as  shown  i  ri  Table  I X— 2 .  This  period  exhibits  highly  variable  salinity 
conditions  both  in  space  and  in  time. 


Salinity  Transect  Information 

116.  Representative  middepth  salinity  values  are  shown  in  Figure  I X  —  1 . 
The  range  over  depth  is  given  to  indicate  the  degree  of  stratification. 
Stratification  effects  are  quite  significant  over  the  mid  and  eastern  section 
of  the  Sound,  even  outside  navigation  channels  during  this  period.  The  well- 
mixed  assumption  does  not  hold  over  these  areas  during  this  period.  As  a  re¬ 
sult,  a  two-dimensional  vertically  averaged  approach  may  not  hold  for  salinity 
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12  JUNE  1980  13  JUNE  1980 


117.  It  should  also  be  mentioned  that  the  biological  fouling  problem 
associated  with  the  conductivity  probe  was  also  evident  during  this  period. 
Instantaneous  salinity  values  were  consistently  2-3  ppt  lower  than  correspond¬ 
ing  salinity  transect  values.  As  a  result,  salinity  transect  values  should  be 
considered  to  be  more  reliable  than  the  continuous  salinity  levels  determined 
at  the  conductivity  temperature  (CT)  stations.  However,  the  time  and  spatial 
variability  characteristic  of  the  continuous  salinity  data  shown  in  Table  IX-2 
is  considered  valid.  As  a  result,  the  initial  salinity  conditions  cannot  be 
specified  using  the  transect  data,  since  these  data  were  obtained  over  2  days 
and  do  not  constitute  a  snapshot  of  the  system.  As  a  result  salinity  will  not 
be  simulated  during  this  period. 


Water  Surface  Elevation  and  Current 
Component  Information 


118.  The  astronomic  data,  constituent  node  factor,  and  equilibrium  argu¬ 
ments  are  shown  in  Table  IX-3  for  conditions  at  the  start  of  the  verification 
period.  The  major  constituents  at  Pascagoula  station  874-1196  are  as  pre¬ 
viously  given.  The  water  levels  are  presented  for  this  station  in  Table  1X-4 
from  20  September,  hr  0,  CST,  for  the  first  2  days  of  the  verification  period. 

119.  In  interpreting  the  elevation  table,  predicted  levels  are  based  on 

the  five  major  constituents  (0^  ,  ,  and  S^).  whereas  the  fil¬ 

tered  levels  contain  all  tidal  energy.  The  unfiltered  levels  contain  meteoro¬ 
logical  input  in  addition  to  all  the  tidal  energy.  In  examining  Table  IX-4, 
the  filtered  and  predicted  levels  are  usually  within  0.1  ft.  Unfiltered 
levels  are  within  0.3  ft  of  filtered  levels.  This  difference  is  consistent 
with  the  wind  information  shown  in  Table  IX— 1 .  No  tropical  storms  or  hurri¬ 
canes  influenced  Mississippi  Sound  during  this  period.  Current  data  for  sta¬ 
tion  10S  were  not  available  for  this  period;  as  a  result  station  V10S  was  not 
considered . 


Table  IX-3 


li 


r. 


Verification  Period  Astronomical  Tide  Conditions,  6  June  1980,  H our  0000 


a.  Astronomic  Terms  at  Greenwich: 


Term 


Longitude  of  lunar  node 
Mean  longitude  of  sun 
Longitude  of  solar  perigee 
Mean  longitude  of  moon 
Longitude  of  lunar  perigee 
Intersection  of  lunar  orbit  and  equator 
Longitude  in  the  lunar  orbit 
Longitude  in  the  celestial  orbit 

V  prime 

V  double  prime 
Solar  hour  angle 


Degrees 

143.26493 

80.47970 

282.60412 

65.91713 

7.65024 

19.55813 

8.58509 

9.21961 

6.02663 

11.20029 

180.00000 


PART  X:  GLOBAL  GRID  HYDRODYNAMIC  CALIBRATION  AND  VERIFICATION 

AND  SYSTEM  MODIFICATION 


120.  In  this  part,  hydrodynamic  simulation  results  over  the  global  grid 
developed  in  Part  VI  are  presented.  Bottom  friction  mechanics  ("Manning’s  n 
versus  stilled  water  depth)  are  calibrated  considering  the  period  20-24  Septem¬ 
ber  1980.  The  period  12-14  June  1980  was  employed  to  verify  the  bottom  fric¬ 
tion.  In  order  to  consider  the  spatial  extent  of  system  induced  changes  in  the 
hydrodynamics,  a  hypothetical  system  modification  in  the  vicinity  of  Sand  Is¬ 
land  was  made  and  a  hydrodynamic  simulation  performed  over  the  period  20-24  Sep¬ 
tember.  Salinity  was  not  considered  in  any  of  these  three  simulations.  As  a 
result,  meteorological  effects  (wind,  surge  setup,  pressure  anomaly)  were  not 
considered.  The  simulations  corresponded  to  solely  astronomic  conditions. 

121.  Program  TIDE  (Appendix  B)  was  incorporated  in  WIFM  as  a  set  of  sub¬ 

routines  in  order  to  develop  water  surface  elevations  at  the  boundary  of  the 
global  grid.  Harmonic  constants  (mean  amplitude  and  Greenwich  epoch)  developed 
by  the  GTM  were  input.  Astronomic  arguments,  node  factors,  and  equilibrium  ar¬ 
guments  are  computed  at  the  start  of  the  simulation  period.  Contributions  from 
each  of  the  five  constituents  (0^  ,  ,  P^  ,  M^  ,  and  S2)  are  summed  to 

produce  the  predicted  or  reconstructed  tidal  elevation  at  the  boundary.  Since 
GTM  grid  spacing  did  not  correspond  to  global  grid  spacings,  it  was  necessary 
to  employ  an  interpolation  procedure.  The  amplitudes  (cm)  and  phases  (“Green¬ 
wich)  are  shown  in  Table  X-l  for  the  14  tidal  signals  located  on  the  GTM  grid 
as  shown  in  Figure  VI-1.  The  convective  acceleration  and  eddy  dispersion  terms 
in  the  motion  equations  were  set  to  zero  on  the  solid  and  open  boundaries  of 
the  global  grid,  thereby  linearizing  the  motion  equations  at  locations  around 
the  boundary.  The  same  approach  was  employed  around  barriers.  Subroutine 
ADVBAR  (Appendix  D)  was  developed  in  order  to  perform  this  task  for  the  model 
user  automatically.  For  global  grid  boundary  cells  (115,  22-59),  values  were 
interpolated  linearly  on  cel  1 -centered  distances  between  the  appropriate  sig¬ 
nals  1-5.  For  global  grid  boundary  cells  (31-114,  59)  values  were  interpolated 
linearly  on  ce 1 1 -centered  distances  between  the  appropriate  signals  6-14. 

122.  Zero  offset  was  assumed  so  that  the  reconstructed  tidal  signal 
oscillated  about  local  mean  sea  level  (LMSL),  which  was  selected  as  the  model 
datum.  Based  upon  Table  1 1 1  —  1 ,  1  ft  was  added  to  all  soundings  to  convert 
their  reference  Gulf  Coast  Low  Water  Datum  (MLLW)  to  the  model  datum  (LMSL). 
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Table  X-l 

Gulf  Tide  Model  Bounda ry  Inputs 


Const i t  ue,it 

Amplitude , 

cm  Phase,  °G 

Amp  1  . 

Mult  .  Fa c  t  u  t 

Signal 

1,  GTM 

Lat  1 ,  Globa  1 

Grid  Cell 

(115,  59) 

°1 

12.4 

18.  1 

1.15 

P, 

3.6 

24.2 

1  .  1<) 

K1 

12.6 

25.4 

1.12 

"2 

1 .8 

136.4 

0 . 72 

S2 

0.7 

120.6 

1  00 

S igna 1 

2,  GTM 

Lat  2,  Global 

Grid  Cell 

( 1 15,56) 

°1 

12.4 

18.3 

,.13 

P 

‘  1 

3.6 

24.5 

:  .  19 

K1 

12.7 

25.7 

1.12 

M2 

1.8 

136.8 

0.72 

S2 

0.8 

0.72 

1 .00 

S i gna  1  3 

,_  GTM 

Lat  3,  Global 

Grid  Cell 

1 1  l 5a  .50  ) 

°1 

12.4 

18.5 

1.15 

P, 

3.6 

24.8 

1.19 

K1 

12.7 

25.9 

1.12 

M2 

1 .8 

137.5 

0.  72 

S2 

0.8 

120. 1 

1  . 00 

S i gna 1  4 

Lat  4l  Global 

Grid  Cell 

(.115,  37) 

°, 

12.5 

18.6 

1 . 00 

P, 

3.6 

25.0 

1 .00 

K1 

12.8 

26.  1 

1 . 00 

^2 

1 .8 

137.2 

1 .00 

S2 

0.8 

119.2 

1 .00 

Signal  5 

,  GTM 

Lat  4,  Global 

Grid  Cell 

(115,  22) 

°1 

12.5 

18.6 

1 .00 

P1 

3.6 

25.0 

1 .00 

K1 

12.8 

26. 1 

1.00 

«2 

1.8 

137.2 

1 .00 

S2 

0.8 

119.2 

1 .00 

(Cont inued) 
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Table  X-l  (Continued) 


Constituent 

Amplitude,  cm 

Phase,  °G  Ampl . 

Mult.  Factor 

Signal  6,  GTM  Long.  1, 

Global  Grid  Cell  (31 ,  59) 

°! 

12 . 3 

17.8 

1 .00 

P1 

3.6 

22.4 

1 .00 

K1 

12.5 

25.  1 

1 .00 

M2 

1 .4 

147.2 

1 .00 

S2 

0.6 

124.8 

1 .00 

Signal  7,  GTM  Long.  1 

,  Global  Grid  Cell  (42,  59) 

0 

1 

12.3 

17.8 

1 .00 

P1 

3.6 

22.4 

1 .00 

K1 

12.5 

25.1 

1.00 

N, 

1.4 

147.2 

1.00 

S 

t 

0.6 

124.8 

1.00 

Signal  8,  GTM  Long.  2, 

Global  Grid  Cell  (57,  59) 

0 

1 

12.3 

17.8 

1.08 

p, 

3.6 

22.6 

1.02 

S 

12.5 

25.1 

1.08 

1 .4 

144.3 

0.76 

S2 

0.7 

123.3 

0.84 

Signal  9,  GTM  Long.  3, 

Global  Grid  Cell  (73,  59) 

°, 

12.3 

17.8 

1.08 

P. 

3.6 

22.8 

1 .20 

K1 

12.5 

25.1 

1.08 

M2 

1 .4 

141 . 3 

0.76 

*2 

0.7 

121  .7 

0.84 

S igna 1  10,  GTM  Long.  4 

i_  Glob  a  1  Grid  Cel  1  (87,  59) 

°. 

12.3 

17.8 

1  .  10 

Pl 

3.6 

23.0 

1 .24 

K1 

12.5 

25. 1 

1  .  12 

»2 

1  .5 

139.2 

0.76 

S„ 

0.7 

120.7 

0.95 

2 

(Continued) 
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Table  X-l  (Concluded) 


onstituent 


Amp  1 i t ude ,  cm  Phase ,  °G 


Amp  1 .  Mu  It. _ Factor 


°1 

P1 

K1 

M2 

S2 

°1 

P1 

K1 

M2 

S2 

°! 

P1 

K1 

M2 

S2 

°1 

P1 

K1 

M2 

S2 


Signa 1  1 1 ,  GTM  Long .  5  ,  Global  Grid  Cell  ( 10  5,  59) 

12.3  17.9  1.10 

3.6  23.2  1.24 

12.5  25.1  1.12 

1.5  137.5  0.76 

0.7  119.9  0.95 

Signal  12,  GTM  Long . _ 6,  Globa  1  Grid  Ceil  ( 11 0 59) 

12.3  17.9  1.15 

3.6  23.4  1.19 

12.5  25.2  1.12 

1.6  136.2  0.72 

0.7  119.4  1.00 

Signal  13,  GTM  Long.  7,  Global  Grid  Cel  1  (112,  59) 

12.3  17.9  1.15 

3.6  2.3.6  1.19 

12.5  25.2  1.12 

1.6  133.5  0.72 

0.7  119.4  1.00 

Signal  14,  GTM  Long.  8,  Globa  1  G r i d  Cell  (1 15,  58) 

12.3  18.0  1.15 

3.6  23.8  1.19 

12.5  25.2  1.12 

1.7  135.2  0.72 

0.7  119.7  1.00 


(3  of  3  sheets) 


123.  All  water  surface  elevations  and  velocity  components  were  initial¬ 
ized  to  zero.  A  cubic  polynomial  feathering  routine  (Appendix  C)  was  developed 
t.i  smoothly  transition  the  boundary  elevations  from  zero  to  their  reconstructed 
1  i'v"  i .  The  convective  acceleration  and  eddy  dispersion  terms  in  the  motion 
'•<pi.it  ions  wore  set  to  zero  on  the  solid  and  open  boundaries  of  the  global  grid, 
’hereby  i  i  him.  icing  the  motion  equations  at  Ioc.it  i  ms  around  the  boundary.  The 
ciine  approach  was  employed  a  round  barriers.  Subroutine  ADVRAK  (Appendix  1) ) 
wa.<  developed  in  order  to  perform  this  task  for  t  li  •  model  use--  automatically. 

124  hilling  tile  calibration  process,  it  was  nocessaiv  to  adjust  the 
amp  1 l t u. b-  "t  tidal  constituents  from  the  values  tabulated  in  fable  X-I  .  The 
GTM  amplitudes  were  multiplied  by  the  factors  shown  in  the  last  lolumn. 

Phases,  as  repented  by  the  GTM ,  were  used  I  i  redly.  In  order  to  develop  the 
appropriate  amplitude  multiplication  factors,  the  harmonic  analysis  results 
obtained  by  'hit  law  (198  5)  and  those  reported  by  Reid  and  Whitaker  (1981)  for 
the  GTM  were  i  impanel  at  t  fie  three  deep  sea  pressure  gages  as  shown  in 
fable  X-2.  As  can  be  observed,  the  Outlaw  (  1  98  5)  .•■suits  were  approximately 
lc-13  percent  higher  tor  the  three  major  diurnal  constituents.  Pressure  sta¬ 
tion  22  diurnal  amplitude  factors  as  reported  in  Table  X-2  were  used  to  modify 
the  diurnal  constituents  of  signals  8  and  9  in  Table  X-l.  Pressure  station  23 
diurnal  amplitude  factors  of  Table  X-2  were  used  to  modify  the  diurnal  con¬ 
stituents  c>  f  signals  10  and  11  in  Table  X-l,  while  diurnal  factors  in  Table  X-2 
for  pressure  station  24  were  used  to  inodifv  the  diurnal  constituents  of  sig¬ 
nals  1-4,  and  12-14  in  Table  X-l.  Since  simulation  results  employing  the  GTM 
results  at  the  boundary  are  compared  with  the  tide  reconstructed  from  the  0^  , 

K  ,  R  ,  M  ,  and  S  ,  constituents  given  by  "itlaw  (198  5)  in  Tables  IT  1-12 
and  1  1  1  -  } ,  tii<-  use  of  t  he  mu  i  t  i  p  I  i  eat  ion  factors  t  <•  m.iK*-  the  GTM  results  con¬ 
sistent  with,  outlaw’s  (  1  98  51  results  is  warranted. 

129.  Initial  attempts  at  employing  this  same  procedure  for  t he  semi¬ 
diurnal  factors  resulted  in  simulated  tide  exceeding  i evens t rue t ed  tides  dur¬ 
ing  the  last  2  days  of  the  r>-dav  (20-24  Sep  1980)  cal  i  brat  i  on  period,  when  the 
semidiurnal  components  were  most  important  in  determining  the  structure  ot  the 
t i de .  In  order  to  obtain  the  best  fit  between  simulated  and  reconstructed 
tides  during  this  period  of  t  tic-  simulation,  it  was  necessary  to  reduce  the  am¬ 
plitudes  of  these  constituents  as  shown  in  Table  X-l  slightly  below  the  levels 
reported  by  Reid  and  Whitaker  (1981).  The  M  constituent  was  effectively 
reduced  by  29  percent,  while  the  S  component  was  essentially  unchanged.  In 


Table  X-2 

Comp a^r i son  of  0 ut  1  aw  (1983)  Analy spp_ pnd 
Re id  and  Whitaker  (1981)  GTM  Results 


Amp  1 i tude  , 

ft 

Phase^ 

°G 

Outj 

aw  Reid 

Factor 

Out  law 

Reid 

Different 

Deepi  Sea 

Pressure 

Station  22 

0.4b 

0.427 

1.0785 

299. 1 

295.36 

3.7 

0.47 

0 . 436 

1.0770 

309.2 

296.45 

12.7 

0. 15 

0. 125 

1.2029 

305.3 

294.65 

10.6 

O 

o 

0.059 

1.5228 

320.6 

332.72 

-12.1 

0 . 05 

0.030 

1.6949 

328.6 

306.90 

21 . 7 

Deep  Sea 

Pressure 

Station  23 

0.46 

0.420 

1.0955 

297.5 

294.86 

2.6 

0.48 

0.430 

1.1168 

308.1 

295.75 

12.3 

0.  15 

0.121 

1.2356 

304.3 

294.15 

10.  1 

0.09 

0.059 

1.5228 

312.0 

328.22 

-16.2 

0.05 

0.026 

1.9084 

321.9 

302.90 

19.0 

Deep  Sea 

Pressure 

Station  24 

0.47 

0.41 

1  .  1463 

296.4 

294.86 

1 .5 

0.47 

0.420 

1.1193 

309.9 

295.65 

13.2 

0.  14 

0.118 

1 . 1854 

300.5 

295.25 

5.2 

0 . 09 

0 . 0b2 

1 . 4446 

311.0 

324.62 

-13.6 

-he  immediate  nearshore  regions,  the  CTM  results  were  used  directly. 

126.  The  depth  versus  Manning's  ri  relationship  shown  in  Table  X - >  was 
al ibrated.  Simulation  results  employing  the  calibrated  relationship  are  pre- 
;ented  in  Figures  X- 1  through  X- 1  1  for  water  surface  e  levut  i  -  n  :.id  in  Fig- 
lies  X-  12  through  X-iK  for  u  and  v  current  component  s  1  >»  hi  starting 

it  20  Septembri  hour  0000  CS 1  .  Results  are  presented  fie  ..  >  •  e  •  n  .  e  eastern 
h.ssissi  pp  •  Sound.  Simulated  water  surface  elevation.  <.  <  >  r ;  <  p  aid  more  closely 
■  |  i  *  ■  1  i  c  t  e .  i  (  ret  oust  ructed)  values  east  of  Gulfport  -i;L  ugh  simulation  re¬ 
al  1 1  s  a  r<-  ar  •  ept  ah  1  e  in  western  Mississippi  Sound,  :  j  • .  ,r  eim-nt  :•  con  1  d  be  male 
in.  'udirg  the  effects  of  the  Lake  Fontrha:  train  sy.-lrr  lie  .  miniated 
iej.t  J  i%  e  raged  current  components  reveal  the  same  sti  ;<  .  e  ;  i  e.  ,  ebb  and 

Table  X-3 

Calibrated  Depth  Versus  Manning's  u  Relationship 

Manning's  n 

Category  1.  Depth  Range,  ft"  Roughness  Factor 


nl 

lL 

1 

0 

- 

s 

0.022 

2 

5 

- 

10 

0.021 

1 

10 

- 

20 

0.020 

4 

20 

- 

2S 

0.01  9 

S 

2r> 

- 

10 

0.01  # 

6 

10 

- 

IS 

0.017 

7 

IS 

- 

40 

0  .  o  !  !• 

y 

4i : 

- 

SO 

1  *  .  0  1  s 

c, 

St 

- 

100 

:  ;  i  i  -i 

10 

!  00 

- 

200 

0  r  i  a 

1  1 

200 

- 

SOO 

C  .  0  12 

i  _/ 

S00 

- 

1000 

;  i.O  11 

1  1 

1000 

- 

sooo 

(i  .  i  ,  Ip 

If  the  cell  water  depth  with  respect  to  mo -lei  •:  itun 
is  gi eater  than  or  equal  to  and  less  than  1  ^ 

the  category  L.  Manning's  n  value  is  assign- 
to  t  fia  t  cell. 


Figure  X- 1 .  Water  surface  elevations,  simulated  and  predicted,  at 
station  874-7437,  20-24  September  1980 
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Figure  X-14.  Velocity  components,  simulated  and  predicted, 
station  V7-S,  20-24  September  1980 
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Veloeity  components,  simulate<i  and  predicted,  at 
station  Vll-M,  20-24  September  1980 


flood  pattern,  as  the  predicted  (reconstructed)  middepth  velocity  stations. 
Since  the  and  tidal  current  components  could  not  be  separated  and 

the  calibration  period  is  outside  the  period  for  which  and  was  deter¬ 

mined  for  all  stations  (refer  to  Table  ill-S),  the  niiddepth  velocity  station 
predicted  (reconstructed)  data  should  be  interpreted  only  in  terms  of  general 
st  ructure . 

127.  In  order  to  verity  the  adjustments  to  the  GTM  constituent  ampli¬ 
tudes  and  the  calibrated  depth  versus  Manning'.',  n  relation,  a  second  period 
12-14  Jun  1)80  was  considered.  During  this  time,  the  structure  of  the  tide 
was  determined  by  the  three  major  diurnal  constituents.  Since  the  character 
e!  t  fie  tide  remained  relatively  constant,  only  5  days  were  simulated.  Ini¬ 
tially  ail  water  surface  elevations  and  current  components  were  set  to  zero. 
Alter  the  first  day,  simulated  levels  corresponded  quite  favorably  to  pre¬ 
dicted  levels.  Plots  of  simulated  versus  predicted  (reconstructed)  water  sur¬ 
face  elevations  are  shown  in  Figures  X-19  through  X-29  for  the  72-hr  period 
starting  12  June  hour  0000  CST .  The  simulated  depth  averaged  currents  are 
compared  with  predicted  (reconstructed)  middepth  values  in  Figures  X-30 
through  X-36,  Simulated  water  surface  elevations  correspond  more  closely  to 
predicted  (reconstructed)  values  east  of  Gulfport,  as  in  the  case  of  the  cali¬ 
bration  period. 

Summary  of  calibra¬ 
tion  ap  d  verifi  c  ja  t  i  o  n 

128.  The  calibration  and  verification  periods  comprise  two  different 

tidal  regimes.  During  the  calibration  period,  the  character  of  the  tide  is 
changing  from  daily  to  semi  daily  and  the  range  is  i ess  than  1.2  ft  (neap  tide). 
Du,  mg  this  period  t !«-  s<  midi  i.-iic  1  mi  <<  ,mp.  >:  c'ct  ?  bminate.  During  the  veri¬ 
fication  period,  t  fie  i  ii.ii  ail  i  ■  o'  t  n<-  :  h  :  t  e  i  at  n.  »• .  v  constant  and  is  daily 
wiili  a  range  greater  1 1:  in  1  . 9  r:  spi  in.;  i  .  m  )  .  During  this  period  the  diurnal 

tidal  components  dominate. 

129.  Since  t  fie  simulated  tidal  i eve  Is  correspond  well  to  the  predicted 
(reconstructed)  results,  the  model  may  fie  confidently  used  to  predict  tidal 
elevations  and  currents  over  the  complete  tidal  ( lunar  1  cycle. 

Hypothetical  Sand  Island  Re¬ 
gional  Dreclge  Material  Disposal  Site 

110.  Consider  a  hypothetical  regional  dredge  material  disposal  site  in 
the  vicinity  of  Sand  Island.  The  water  depths  with  respect  to  local  mean  sea 


Figure  X-21 


Water  surface  elevations,  simulated  and  predicted 
at  station  874-6819,  12-14  June  1980 


Figure  X-22.  Water  surface  elevations,  simulated  and  predicted, 
at  station  874-4986,  12-14  June  1980 
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gure  X-23.  Water  surface  elevations,  simulated  and  predicted 
at  station  874-3735,  12-14  June  1980 


12  21  n  IB  50  72  81  S6  '.CB  U 

HGUfrS 


Figure  X-24.  Water  surface  elevations,  simulated  and  predicted 
at  st.it  ion  874-3081,  12-14  June  1 9  80 


Figure  X-25 .  Water  surface  elevations,  simulated  and  predicted, 
at  station  874-2221,  12-14  June  1980 


igure  X-29.  Water  surface  elevations,  simulated  and  predicted, 
at  station  873-5587,  12-14  June  1980 
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Velocity  components,  simulated  and  predicted,  at 
station  V13-S,  20-24  September  1980 


level  are  modified  as  shown  in  Table  X-4  in  the  vicinity  of  Sand  Island  to  re¬ 
flect  accumulated  dumping  of  dredged  materials.  The  regional  site  encompasses 
an  area  of  approximately  4.35  square  miles.  The  modified  depths  over  this 
area  reflect  a  disposal  of  approximately  1.27  billion  cubic  feet  (47.03  mil¬ 
lion  cubic  yards)  of  dredged  material.  To  place  this  figure  in  perspective, 
over  the  8-year  period  (1970-1979)  an  average  of  13.7  million  cubic  yards  per 
year  were  removed  from  Mississippi  Sound  (USACE,  Mobile  District  1979). 

131.  The  20-24  September  1980  period  as  employed  in  the  hydrodynamic 
calibration  was  simulated  assuming  the  hypothetical  Sand  Island  regional  dis¬ 
posal  site  was  in  place.  Simulated  water  levels  under  this  alternative  at 

11  stations  were  compared  with  the  calibration  simulation  levels  previously 
determined.  Water  levels  within  the  Sound  are  unaltered  by  this  hypothetical 
offshore  disposal  site  even  at  grid  cell  (87,  27). 

132.  The  detailed  flow  patterns  in  the  vicinity  of  Sand  Island  are 
shown  in  Tables  X-5  and  X-6  for  the  present  conditions  at  hr  72  and  120,  re¬ 
spectively.  The  flow  patterns  in  the  vicinity  of  the  hypothetical  site  are 
shown  in  Tables  X-7  and  X-8  at  hr  72  and  120,  respectively.  Comparing  Tables 
X-5  and  X-f>  and  Tables  X-7  and  X-8,  one  notes  that  the  detailed  flow  structure 
is  changed  only  in  the  immediate  vicinity  of  the  hypothetical  dredge  material 
disposal  site.  Changes  in  flow  structure  induced  by  system  modifications  in 
Mississippi  Sound  and  adjacent  waters  similar  to  the  hypothetical  Sand  Island 
site  will  be  felt  only  in  the  immediate  vicinity  and  will  not  propagate  to  the 
boundary.  Therefore,  these  modifications  may  be  effectively  studied  using  the 
global  grid  developed  in  this  study. 


Table  X-4 


Globa  1  Grid  Alternative 
(Sand  Island  Complex) 

Original  Depth  Alternative 


Grid  Cell 

ft* 

Depth , 

(88,  31) 

13.0 

5 

(88,  32) 

30.0 

5 

m 

00 

11.0 

5 

(89,  32) 

19.0 

5 

(90,  31) 

12.0 

5 

(90,  32) 

8.0 

5 

With  respect 

to  local  mean 

sea  leve 1 . 
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SUPPLEMENTAL  NOTES  FOR  TABLES  X-5  THROUGH  X-8 


IN  ORDER  TO  INTERPRET  TABLES  X-5  THROUGH  X-8  CONSIDER  THE  GRID 
NOTATION  IN  THE  FOLLOWING  FIGURE 


M  X.  U 


THEREFORE  <N  TABLE  X--  rOP  COWpU  T  A  ' .  ON  A,  cEi  .  '>5  ...  ThE  Y 

VELOCITY  COMPONENT  IS  -0  i  t  F ...  Whil  f  THE  X  VE;  OC  .  1  V  COMPONENT 
iS  5  F«S 


OBSERVE  IN  TABLES  X  '  A  NO  x-M  rHfc  »v  K'T  HE  Tu  AL  ORE  P'  TE  DISPOSAL 
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NOTE  IN  ’HE  i'.'-'l  T--  •ICN--  A;  ■  f  .  T  If  S  NORMA,.  I  THISLINE 
/.  i  E  /f.Rf ) 
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Table  X-7 

Simulated  Global  Grid  Velocity  Components  (fps  x  ]Q)  for 
the  Sand  Island  Alternative  at  Hour  72 
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PART  XI:  REFINED  GRID  METHODOLOGY  AND  APPLICATION 


133.  The  global  grid  previously  developed  is  to  be  used  to  provide 
boundary  information  for  a  set  of  refined  grids.  Each  refined  grid  is  to  be 
located  around  a  potential  channel  alteration  project.  Detailed  circulation 
in  the  vicinity  of  the  proposed  channel  modifications  is  to  be  studied  within 
the  refined  grid.  The  assumption  underlying  this  approach  is  that  alternative 
channel  alignments  and/or  deepening  effects  will  be  localized  and  will  affect 
circulation  only  within  tfie  area  covered  on  the  refined  grid. 

1  14.  The  width  of  the  principal  navigation  channels  within  Mississippi 
Sound  is  from  2 lid  to  330  ft  .  It  is  desired  to  use  space  steps  on  the  same 
otd<  r  i n  the  ref i ned  grid.  Since  the  maximum  depth  even  under  proposed  deep¬ 
ening  alternatives  is  approximately  40  ft,  the  gravity  wave  speed  in  the  re¬ 
fund  grid  area  is  less  than  38  fps .  Therefore,  the  explicit  time  step  limit 
is  less  than  9  sec.  We  attempt  to  employ  an  implicit  time  step  of  60  sec. 

i  3.3 .  Consider  a  refined  grid  system  such  that  when  a  given  time  period 
is  simulated,  the  computational  cost  incurred  will  not  exceed  the  amount  re¬ 
quired  to  simulate  the  same  time  period  on  the  global  grid.  Observe  on  the 
6783-cell  global  grid  a  time  step  of  360  sec  was  utilized.  Thus  to  advance 
'he  simulation  one  time  step  on  the  global  grid,  six  time  steps  will  need  to 
fie  performed  on  the  refined  gr  id. 

136.  The  cost  involved  in  inverting  a  matrix  by  the  Thomas  algorithm  as 
employed  in  WJFM  is  approximately  proportional  to  its  rank  cubed.  Estimate 
t  he  rank  c.  f  the  matrix  to  be  inverted  foi  a  given  grid  setup  as  the  square  ron 
v-t.  the  nnmbei  of  cells  i  t.  the  grid,  N  ,  and  obtain  the  following  relationship 
ur'niling  a  unit  pr<>|  rti  tta  1  ,  ;  v  factor: 

W  -  (N  )  1/1  (X!  .  1  ) 


who  r  i 

W  -•  .  I  i  nv.  'red  in  inverting  the  matrix  for  a  given  grid  system 
N  -  number  of  rolls  in  the  grid  system 
The  total  s  i  mn  1  a  t  j  o,i,  rrst,  (’  ,  is  then  given  by  the  following  relation 


C  --  MW 


(XI  .2) 


where 


C  =  total  simulation  cost 

M  =  number  of  time  steps  in  the  simulation 
W  =  cost  involved  in  inverting  the  grid  system  matrix 
Let  us  equalize  simulation  cost  in  considering  a  given  time  period  on  each 
grid.  In  so  doing,  an  upper  bound  on  the  number  of  grid  cells  to  employ  in 
the  refined  grid  is  obtained  as  follows 


/N  \3/2 

MgW8  =  MrWr  =  (6M8)^j  W8  (XI.  3) 

Therefore  if  approximately  5000  cells  in  the  global  grid  are  involved  in  the 

=  1  and  N  ~  1500  .  Thus  a  refined  grid  involving 
approximately  1500  cells  in  the  computations  must  be  constructed. 

137.  Let  us  now  consider  the  general  orientation  of  the  principal  navi¬ 
gation  channels  to  the  axes  in  the  global  grid  as  shown  in  Figure  XI-1.  As 
may  be  observed,  the  navigation  channels  do  not  run  parallel,  throughout  their 
extent,  to  the  global  grid  axis.  If  the  refined  grid  orientation  corresponds 
to  that  of  the  global  grid,  stair-stepping  is  normally  employed  to  resolve  the 
channel  sections  nonparallel  to  the  grid  systems.  This  approach,  however, 
would  necessitate  the  reduction  of  the  space  step  by  more  than  100  percent. 
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Figure  XI-1.  General  orientation  of  major  navigation 
channels  in  Mississippi  Sound 
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As  a  result  less  than  750  cells  could  be  used  in  the  refined  grid  system  to 
maintain  a  simulation  cost  for  a  given  time  equal  to  that  of  the  global  grid. 
This  number  of  cells  would  probably  be  insufficient  to  resolve  the  channel 
properly  in  the  refined  grid.  Therefore,  it  is  necessary  to  idealize  the 
channel  systems  somewhat  within  the  refined  grid.  In  order  to  demonstrate  the 
refined  grid  approach  on  a  typical  navigation  channel,  the  Pascagoula  Channel 
was  selected  for  further  study.  A  refined  grid  was  constructed  in  the  ideal¬ 
ized  manner  shown  in  Figure  XI-2.  By  utilizing  this  channel  idealization  it 
was  possible  to  employ  the  same  orientation  in  the  refined  grid  as  in  the 
global  grid  and  use  two  dense  vertical  bands  in  the  vi  inity  of  the  vertical 
channel  sections  and  a  single  dense  band  of  cells  in  the  vicinity  of  the 
horizontal  channel  section.  The  details  of  the  development  of  this  refined 
grid  system  in  the  vicinity  of  the  Pascagoula  Channel  are  now  presented. 
Pascagoula  Channel  refined  j^rid 

138.  Nautical  Chart  11374  at  a  scale  of  1:40000  was  used  in  mapping  the 
refined  grid. 


N-S  M-S 

Figure  XI-/.  Pascagoula  Channel  configuration 

139.  The  mapping  for  the  horizontal  (east-west)  orientation  is  pre¬ 
sented  below  in  I  able  XI- 1. 
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Table  XI - 1 


East-West  Orientation 


NAPPING 

C0NPIETE i 

1« 

REGIONS 

REG 

LPR 

REAL 

SPACE 

ALPHA 

SPACE 

GIVEN 

X  PR  I  HE 

CALC 

X  PRINE 

A 

1 

C 

1 

2 

1 .9990 

1 

1 .50000 

1.50900 

- . 5000006*00 

.  150909E«91 

.1000006*01 

2 

7 

4.0000 

3 

1.50009 

1.50009 

- . 6450486*01 

.6511586*01 

.4306026*00 

3 

2 

n . leee 

19 

.75000 

.75573 

. 1S9673E+02 

-.  173763E*93 

-.1552676*01 

4 

2 

12.3908 

12 

.50909 

.47451 

. 1299946*02 

-. 4278716*09 

-.8141936*01 

s 

2 

12.8900 

14 

.19090 

.11594 

.1413026*02 

-.3329746*92 

-.122919E*9i 

6 

2 

13.0999 

ie 

.  19090 

.08619 

.1291646*02 

.1134766-20 

.  164994E«02 

7 

2 

13.5909 

18 

.59990 

.53493 

.1212366*02 

.2275256-08 

.6995856*01 

8 

S 

15.9099 

20 

1.90999 

1.00613 

-.1016456*04 

.9728956*03 

.  195090E-01 

9 

4 

19.5990 

25 

.89900 

.80841 

.1556826*92 

.2563456-06 

.5140256*01 

ie 

0 

24.9900 

29 

1.59090 

1.49453 

0. 

0. 

0. 

The  mapping  for  the  vertical  (north-south)  orientation  is  given  in  Table  XI-2. 


Table  XI-2 

North-South  Orientation 


NAPPING 

COMPLETE i 

IS 

REGIONS 

REG 

LPR 

REAL 

SPACE 

ALPHA 

SPACE 

GIVEN 

X  PRINE 

CALC 

X  PRINE 

A 

8 

C 

1 

4 

1.0999 

1 

2.10909 

2.10999 

-.1571 66E*9 1 

.  257166E*91 

.8165926*00 

2 

2 

8.0000 

5 

1.56990  , 

1 . 56323, 

•337351E*02 

-.41 9579E  *92 

-.3037156*09 

3 

2 

10.5009 

7 

1.00000 

1.09812 

. 7685856*92 

- .816147E*02 

- . 1063446*00 

4 

3 

12.2509 

9 

.75000 

.76342 

.270150E*02 

-.41 0470E  *02 

-.465341E*00 

5 

2 

14.1009 

12 

.50000' 

.50982 

.  796248E*01 

. 5  38587E*00 

. 979292E*00 

6 

2 

15.1999 

14 

.59000 

.49922 

.1S7612E*02 

- . 862203E  *12 

-. 195706E*02 

7 

2 

15.6909 

16 

.10000 

.10648 

. 45  4964E  *02 

- . 349252E*02 

- .57 1 698E-01 

B 

3 

15.8009 

18 

.10000 

.09402 

. 1S6492E*02 

. 1 Z2977E- 1 4 

. 1 1 2259E*92 

9 

5 

16.5009 

21 

.  45009 

.4S479 

. 1 096416*02 

.  702285E-01 

.1483946*91 

10 

3 

18.9909 

26 

. S9900 

.59431 

.21 4908E  *02 

- .  375466E  *08 

-.506097E+01 

11 

3 

29.0000 

29 

.25909 

.26017 

.298723E*02 

- . 3889946*1 3 

-.8649696+01 

12 

9 

20.5099 

32 

.10999 

.19063 

.1711616*02 

. 1 25073E*00 

.951S71E+09 

13 

3 

21.4000 

41 

.19909 

.09943 

•212172E*02 

.  1 99S64E-36 

.222989E+02 

14 

6 

22.1000 

44 

.40999 

.44739 

.291 31 6E*92 

. 721 794E-16 

.100097E+02 

15 

9 

27.2000 

50 

1.49990 

1.41377 

9. 

0. 

9. 

The  grid  obtained  contains  49  *  28  =  1372  cells  and  is  presented  in  Fig¬ 
ure  XI-.*.  In  the  above  tables,  the  real  space  values  are  in  map  inches. 

140.  Program  'I  OR  11)  was  employed  to  plot  the  relined  grid  at  a  scale  of 
1:40000  so  that  it  would  be  overlaid  on  Chart  11374.  The  average  depth  over 
tin  cell  was  assigned  as  the  cel]  depth  in  most  cases.  The  hydrographic  survey 
information  was  also  used  to  modify  channel  depths  and  chart  depths  in  the  vi¬ 
cinity  of  the  passes.  The  assigned  depths  are  presented  in  Table  XI-3,  which 
was  directly  output  from  the  model.  All  water  depths  are  in  feet  preceded  by 

a  minus  sign  and  are  with  respect  to  local  mean  sea  level.  All  land  is  repre¬ 
sented  as  +10  ft . 

Ba mt i e r  I s 1  a nd  conf igurat ion 

141.  All  barrier  islands  were  located  on  cell  faces  and  assigned  ^n 


GULF  OF  MEXICO 


(font  in 


elevation  of  10  ft  above  local  mean  sea  level.  The  barriers  were  typed  as  ex 
posed  since  no  overtopping  occurred.  The  barrier  island  configuration  is 
shown  in  Table  XI-4  below.  The  orientation  numbered  1  corresponds  to  a  barrier 
on  the  u-cell  face,  while  orientation  2  represents  a  v-face  barrier. 


Table  Xl-4 

Ba  r  r  i _er  Con  f  i  gu  ration 


Loca  t i on 

No . 

Orientation 

N 

M 

1 

2 

7 

2 

2 

2 

7 

i 

2 

8 

- 

i  :$ 

•.) 

; '  i 

f. 

i 

10 

1  i) 

I 

*-♦ 

2  a 

8 

i 

9 

i 

24 

10 

i 

'1 

2  5 

1  1 

! 

10 

25 

12 

1 

1 1 

25 

13 

1 

12 

25 

14 

1 

13 

25 

15 

1 

14 

25 

16 

1 

15 

25 

17 

1 

16 

25 

18 

1 

30 

25 

19 

1 

31 

25 

20 

1 

43 

26 

21 

1 

44 

26 

22 

1 

45 

26 

23 

1 

26 

26 

24 

I 

47 

26 

25 

1 

48 

26 

Mow  inputs  and  calibration  stations 

142.  Two  f low  inputs  for  the  Pascagoula  River  System  are  considered. 

The  Pascagoula  Rivet  is  assigned  as  input  to  cell  (8,1)  and  the  Past  Pascagoula 
River  is  assign'd  as  input  to  cell  (16,1).  Average  daily  USGS  flows  corrected 
by  drainage-  area  ratios  were  used  to  specify  input  flows  in  c f s . 

14'!.  llie  location  of  tidal  stations,  velocity  stations,  and  salinity 
transect  stations  on  the  refined  grid  system  are  shown  in  Table  Xl-5. 
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Table  XI-5 


Calibration  Station  Locations  on  the  Refined  Grid 


_ Tide  Station 

T4 

T5 

T6 

T7 

Velocity  Sta t  ion 

V 1  2 
V]  3 
V 14 
V15 
VI 6 

Salinity  Transect  Station 

T54 

T64 

T62 

T66 

T60 

T68 

T56 


WIFM  Grid  Coordinates 


(4,24) 

(25.6) 

(49,27) 

(49.8) 

WIFM  Grid  Coordinates 

(3.20) 

(3.3.20) 

(24.25) 

(30,27), (33, 26) 
(49,22) ,(49,23) 

WIFM  Grid  Coordinates 

(8,22) 

(17.6) 

(24.9) 

(31.7) 

(33,16) 

(49,19) 

(33.26) 


1  < > ! > . 1 1  grid  boundary  interface 

144.  A  subgrid  of  the  global  grid  encompassing  the  Pascagoula  Channel 
.(cm  was  plot  teil  at  a  scale  of  1:40000.  The  global  grid  subgrid  was  over- 
:  :  • !  the  relined  grid  and  boundary  cells  in  the  global  grid  surrounding  the 
:  ea  ,ii  e  unsigned  t<>  the  nearest  corresponding  cell  center  of  the  refined 
i  1  1  tie  assignments  are  shown  i  ri  Table  Xl-0.  During  execution  of  the 

•  i  ‘  gild  simulation  water  surface  elevations  are  written  to  logical  unit 
under  .  and  salinities  to  logical  unit  number  35.  During  execution  of  the 
tine. I  grid  simulation  these  two  files  are  accessed,  and  linear  in  tune  and 
.  .  :  -  .  en  i  e  i  e.|  line. ii  in  distance  interpolations  are  performed  to  determine  the 
.  t  , ned  giid  water  surface  elevation  and  salinities  at  the  boundary. 
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Table  Xl-6 

Global  Grid  Cell-Refined  Grid  Cell  Boundary  Assignment 


PART  XII:  REFINED  GRID  HYDRODYNAMIC  CALIBRATION  AND 
CHANNEL  MODIFICATION 


145.  In  this  part,  two  hydrodynamic  simulations  over  the  Pascagoula 
Channel  refined  grid  developed  in  Part  XI  are  considered.  Water  surface  ele¬ 
vations  computed  on  the  global  grid  for  the  20-24  September  1980  calibration 
period  were  accessed  to  define  the  boundary  conditions  along  the  refined  grid. 
The  convective  acceleration  and  eddy  dispersion  (advective  terms)  in  the  motion 
equation  were  not  considered  along  the  open  and  closed  boundaries  and  around 
barriers  (Appendix  B)  as  well  as  in  the  two  dense  vertical  bands  used  to  de¬ 
scribe  the  Pascagoula  Channel  and  in  the  one  dense  horizontal  band  used  to 
represent  a  section  of  the  Pascagoula  Channel.  The  sudden  discontinuity  in 
the  depth  profile  outside  the  navigation  channel  was  thought  to  perhaps  cause 
problems  in  the  advective  calculations.  For  this  reason,  advection  was  not 
considered  in  the  vicinity  of  the  navigation  channels. 

146.  In  simulation  one,  the  same  Manning’s  n  versus  water  depth  rela¬ 
tions  considered  in  the  global  grid  hydrodynamic  calibration  are  used.  The 
computed  (A)  versus  predicted  (reconstructed)  water  surface  elevations  are 
presented  in  Figures  XII-1  through  XII-4  for  the  120-hr  period  starting  20  Sep¬ 
tember  hour  0000  CST.  Simulated  currents  are  shown  in  Figures  XII-5  through 
XII-11.  No  predicted  (reconstructed)  currents  are  shown,  since  for  this  pe¬ 
riod,  due  to  the  K^-P^  separation  problem,  they  cannot  be  determined.  The 
simulated  water  surface  elevations  match  very  closely  the  predicted  (recon¬ 
structed)  levels.  From  this  we  may  infer  that  the  global  grid  boundary  eleva¬ 
tions  are  accurate  and  that  the  L-shaped  representation  of  the  Pascagoula  Chan¬ 
nel  is  sufficient  to  describe  water  surface  elevations  within  the  overall  chan¬ 
nel  system. 

Channel  modification 

147.  In  simulation  two,  the  Pascagoula  Channel  was  considered  to  be  in¬ 
creased  in  width  from  330  to  660  ft.  The  depth  of  the  widened  channel  was 
maintained  at  39  ft  with  respect  to  local  mean  sea  level  (model  datum).  The 
changes  in  depth  required  to  update  the  original  channel  width  are  shown  in 
Table  XII-1.  Hydrodynamics  was  considered  over  the  same  5-day  period  used  in 
simulation  one.  No  wind  effects  were  considered  and  salinity  patterns  were 
not  simulated. 

148.  Simulated  water  levels  for  the  channel  modification  are  presented 
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Figure  XII-1.  Simulated  versus  reconstructed  water  surface  elevations 
at  station  874-0199,  20-24  September  1980  (330-ft  channel) 
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Figure  XII-2.  Simulated  versus  reconstructed  water  surface  elevations 
at  station  874-0405,  20-24  September  1980  (330-ft  channel) 
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Figure  XII-3.  Simulated  versus  reconstructed  water  surface  elevations 
at  station  874-1196,  20-24  September  1980  (330-ft  channel) 


Figure  XII-4.  Simulated  versus  reconstructed  water  surface  elevations 
at  station  874-2221,  20-24  September  1980  (330-ft  channel) 
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Figure  XII-10.  Simulated  current  velocities  at  station  35  20, 
20-24  September  1980  (330-ft  channel) 
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Table  XII-1 

Pascagoula  Channel  Widening 


Grid  Cell  Range 

(16,  1-15) 
(17,  2-15) 
(18-34,  14) 
(18-34,  15) 
(35,  15-28) 


Channel  Depth,  ft 

38.0 

38.0 

38.0 

38.0 

38.0 


versus  predicted  (reconstructed)  levels  for  the  unmodified  system  in  Fig¬ 
ures  XII-12  through  XI I  —  1 5 .  Water  surface  elevations  with  the  widened  channel 
corresponded  to  the  predicted  (reconstructed)  levels  for  the  unmodified  system 
and  hence  to  the  simulation  one  results. 

149.  The  simulated  current  structure  at  hr  72  of  the  5-day  simulation 
period  is  presented  in  Table  XII-2  for  the  330-ft  channel  (simulation  one)  and 
in  Table  XI 1-3  for  the  660-ft  channel  (simulation  two).  In  considering 
Tables  XII-2  and  XI 1-3,  the  reader  should  consult  the  note  on  page  174.  Ob¬ 
serve  the  outline  of  the  channel  is  indicated  as  a  set  of  solid  lines  in  both 
tables.  The  impact  of  widening  the  channel  is  felt  locally  in  the  vicinity  of 
the  channel  (2-3  cells  in  width),  and  generally  the  changes  are  extremely 
small. 

150.  Based  upon  these  results,  the  nested  (refined  grid)  approach  is 
appropriate  for  studying  or  identifying  the  extent  of  channel  alterations  on 
the  local  circulation  pattern.  To  study  circulation  patterns  within  the  chan¬ 
nel  area,  a  three-dimensional  or  two-dimensional  laterally  averaged  model  is 
needed . 
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Figure  XII-12.  Computed  and  reconstructed  water  surface  elevations 
at  station  874-2221,  20-24  September  1980  (660-ft  channel) 


Figure  XI I  —  1 3 .  Computed  and  reconstructed  water  surface  elevations 
at  station  874-1196,  20-24  September  1980  (660-ft  channel) 
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Figure  XII-14.  Computed  and  reconstructed  water  surface  elevations 
at  station  874-0405,  20-24  September  1980  (660-ft  channel) 
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PART  XIII:  SALINITY  CALIBRATION  ON  THE  GLOBAL  AND 
REFINED  GRID  SYSTEMS 


151.  The  simulation  of  salinity  requires  the  consideration  of  meteoro¬ 
logical  as  well  as  astronomical  effects.  Atmospheric  pressure  anomalies  were 
not  considered  in  this  study.  Predominate  meteorological  effects  in  most 
cases  are  due  to  wind.  Therefore,  the  wind  pattern  must  be  studied  for  the 
period  of  consideration.  The  wind  pattern  is  applied  to  develop  a  surface 
shear  stress,  which  affects  the  hydrodynamic  computations.  Wind  effects  were 
not  considered  in  the  GTM  simulations  used  to  specify  elevations  along  the 
boundary  of  the  global  grid.  In  performing  combined  global  and  refined  grid 
simulations,  the  same  wind  pattern  applied  over  the  global  grid  is  applied 
over  the  refined  grid  to  develop  the  water  surface  shear  stress  effects  over 
this  grid.  During  the  global  grid  simulation,  water  surface  elevations  and 
salinities  are  saved  at  various  points  within  the  global  grid,  which  encompass 
the  refined  grid  area.  During  the  simulation  on  the  refined  grid,  the  pre¬ 
viously  saved  water  surface  elevations  and  salinities  are  accessed  to  develop 
the  appropriate  boundary  conditions. 

Wind  stress  relations 

152.  The  following  form  of  wind  stress  is  employed 

t  =  Cnp  V2  (XIII. 1) 

w  Da  w 

where 

T  =  wind  stress  on  the  water  surface 
w 

Cp  =  dimensionless  drag  coefficient 

V  =  wind  speed  at  the  anemometer  level 
w  r 

p  =  air  density 
a 

153.  Several  relationships  are  available  for  the  determination  of  . 

Wilson  (1962)  employs  =  0.0015  for  light  wind  conditions.  Van  Dorn  (1953) 
employs  =  0.000923  for  winds  less  than  or  equal  to  14  knots. 

154.  Garret  (1977)  suggests  the  following  power  law  relation 

Cn  x  io3  =  0.51V0’46  (XIII. 2) 

L)  W 
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where  V  is  measured  at  the  anemometer  level  in  knots.  (Note:  1  knot 
w 

=  1.6884  fps  =  0.5146  mps . )  Evaluating  XIII-2  at  14  knots,  one  obtains 
=  0.00126  . 

155.  We  observe  the  units  of  stress  in  Equation  XIII-1  above  are  force 

2 

per  unit  area  (M/LT  ).  In  the  motion  equations  in  the  x  and  y  directions, 

2 

the  units  are  in  terms  of  acceleration  (L/T  ).  Thus,  in  applying  Equation 
XIII-1  to  a  motion  equation,  one  divides  the  stress  by  the  product  of  water 
density  and  depth;  e.g.,  (L/T^)  =  (M/LT^)/ (M/L^)  .  Formally, 


a 

w 


CDPaVw  _  K  2 
pD  D  w 


(XIII. 3) 


where 


a  =  acceleration  induced  by  wind  stress 

D  =  water  depth 

p  =  air  density 
a 

p  =  wafer  density 

V  =  wind  speed  at  the  anemometer  level 
w  r 

K  =  dimensionless  coefficient 


If  we  use  C  =  0.001  ,  we  obtain  K  =  1.191  *  10  ^  .  If  we  employ  British 

u  2 
units  (slug,  ft,  sec),  and  V  is  given  in  mph  ,  we  multiply  K  by  (1.4667) 

-6  w 

to  obtain  2.562  *  10  .  This  is  the  coefficient  employed  in  the  computer 

code . 

Dispersion  coe f f ic i ents 


156.  The  effective  dispersion  coefficients  are  assumed  to  have  the  form 
given  in  Equation  V.64.  The  dispersion  offsets  due  to  wind  are  assumed  equal 
to  zero. 


Globa  1  Grid  Salinity  Simulations 


12-16  June  1980 

157.  This  period  was  originally  considered  as  a  verification  of  the 
salinity  computations.  Initial  conditions  in  Mississippi  Sound  were  developed 
on  the  global  grid  based  on  salinity  transect  data.  In  order  to  verify  the 
effective  dispersion  coefficients,  model  results  were  to  be  compared  with  con¬ 
tinuous  salinity  data  at  the  end  of  the  5-day  period.  Upon  further  study  of 
the  continuous  salinity  data  and  comparison  with  the  salinity  transect  data 
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at  common  time  and  spatial  location,  it  appeared  that  the  continuous  data  were 
always  lower  than  the  transect  data.  Raytheon  Ocean  Systems  acknowledged  a 
biological  fouling  problem  with  the  conductivity-temperature  probe.  Continu¬ 
ous  salinity  values  always  increased  drastically  after  the  meters  were  ser¬ 
viced.  As  a  result  of  these  data  uncertainties,  it  was  decided  not  to  con¬ 
sider  this  period  in  the  salinity  study. 

20-24  September  1980 

158.  This  period  was  selected  for  calibration  of  the  dispersion  coef¬ 
ficients  in  the  salinity  algorithm.  Salinity  transect  values  were  available 
on  20  and  21  September.  These  values  were  located  on  the  global  grid  and  two 
rectangular  areas  were  set  up  in  which  salinity  values  were  visually  inter¬ 
polated  from  the  located  transect  values.  National  Marine  Fisheries  (1981) 
data  were  obtained  for  cruises  #106  (April  1980)  and  #112  (November  1980)  of 
the  OREGON  II.  These  data  provided  a  general  understanding  of  salinity  pat¬ 
terns  in  the  vicinity  of  the  Mississippi  Delta.  A  deep  sea  vertically  aver¬ 
aged  value  of  36  parts  per  thousand  was  employed. 

159.  Initial  conditions  were  assigned  in  a  three-step  process  as  shown 
in  Table  XI 1 1  —  1 .  In  step  one,  values  were  assigned  based  on  cell  water  depth. 
In  step  two,  salinity  values  were  specified  within  Mississippi  Sound  based  on 
salinity  transect  data.  In  step  three,  initial  salinity  values  within  Lake 
Borgne  were  specified  in  a  zone  format.  In  this  process  each  succeeding  step 
overrides  the  previous  step. 

160.  Salinity  boundary  conditions  which  remained  constant  over  time  are 
shown  in  Table  XI 1 1  —  2 .  A  cell-centered  spatial  interpolation  analogous  to  that 
employed  for  water  surface  elevations  was  used  to  determine  salinity  values 
along  the  seaward  boundary. 

161.  Meteorological  effects  were  treated  in  the  following  manner.  The 
surge  setup  due  to  tropical  storm  Hermine  was  approximately  in  the  range 

0. 1-0.5  ft.  The  effect  causes  an  increase  or  shift  in  the  nodal  datum  cor¬ 
responding  to  local  mean  sea  level  equal  to  the  setup,  which  is  less  than 
0.5  ft.  Since  soundings  are  not  known  to  this  accuracy  in  most  areas  of  the 
Sound,  the  surge  setup  effect  is  negligible  and  was  not  considered.  However, 
if  one  attempts  to  compare  simulated  water  levels  with  the  unfiltered  (raw) 
elevation  data,  then  the  datum  of  simulated  levels  must  be  adjusted  by  this 
surge  setup.  Wind  data  over  the  period  are  presented  in  Table  X 1 1 1  —  3 .  Based 
upon  a  study  of  the  spatial  and  temporal  wind  variations  shown  in  Table  XI 1 1  —  3 , 
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Table  XIII-1 

Initial  Salinity  Conditions  on 
the  Global  Grid 


Water  Depth 
ft 

Initial  Salinity 
Value,  ppt 

0  - 

10 

22.0 

10  - 

20 

23.0 

20  - 

30 

25.0 

30  - 

50 

30.0 

50  - 

75 

34.0 

75  - 

100 

34 . 3 

100  - 

120 

34.5 

120  - 

200 

35.0 

200  - 

300 

35.5 

300  - 

5000 

36.0 

Salin'ty  Grid  Cel  1  by^  Grid  Cell 
Int e  rpo 1  a  t ed  L i m 1 1 s 

GlobaJ  Grid  Cell  Range 
Patch  N  _  M 

1  15-27  19  -  39 

2  28-87  15  -  32 

Sa  1  in  i  t_y  Zone  Spec  i  fjed  I_n  i  t  i  a  1  Condpt  ions 

G 1  oba  1  Cell  Range  Salinity 

Zone  _ N_  M  ppt _ 

1  1-15  33-50  15.0 
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Table  XIII-2 

Global  Grid  Boundary  Salinity  Conditions 


Tidal 

Global 

Salinity 

Signal 

Grid  Cell 

Value,  pp 

1 

(115,  58) 

36.0 

2 

(115,  56) 

36.0 

3 

(115,  50) 

36.0 

4 

(115,  37) 

36.0 

5 

(115,  22) 

34.0 

6 

(  31,  59) 

30.0 

7 

(  42,  59) 

36.0 

8 

(  57,  59) 

36.0 

9 

(  73,  59) 

36.0 

10 

(  87,  59) 

36.0 

11 

(103,  59) 

36.0 

12 

(110,  59) 

36.0 

13 

(112,  59) 

36.0 

14 

(115,  59) 

36.0 

F  resh- 

water 

Global 

Salinity 

Inflow 

Grid  Cell* 

Value,  pp 

(59,  19) 
(59,  17) 
(13,  33) 
(19,  20) 
(33,  15) 


Refer  to  Table  VI-9. 
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Wind  Data  for  20-24  September  1980 
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wind  speed  and  direction  were  spatially  averaged  at  6-hr  intervals. 

162.  A  time  step  of  360  sec  was  employed  1200  times  in  order  to  simu¬ 
late  a  120-hr  period  starting  at  hour  0000  CST,  20  September.  In  order  to  re¬ 
move  the  effects  of  transients  in  the  hydrodynamics  from  the  salinity  compu¬ 
tations,  salinity  computations  were  initiated  after  240  time  steps  (1  day  pro¬ 
totype  time).  Thus  salinity  computations  were  performed  over  the  period 
21-24  September.  The  initial  conditions  specified  over  the  western  portion  of 
the  Sound  corresponded  to  the  transect  data  taken  on  20  September.  Since 
salinity  time  variations  in  salinity  were  extremely  small  during  this  period, 
it  was  assumed  that  conditions  on  21  September  were  effectively  equal  to  con¬ 
ditions  on  20  September.  Simulation  results  at  the  end  of  the  simulation 
(hour  2400  CST  24  September)  were  compared  with  transect  data  obtained  over 
the  western  half  of  the  Sound  during  the  day  of  24  September  and  over  the 
eastern  half  of  the  Sound  during  the  day  of  25  September.  Since  time  varia¬ 
tions  in  salinity  were  extremely  small,  this  was  felt  to  be  a  valid  procedure. 
Wind  information  input  at  6-hr  intervals  was  interpolated  in  time  at  each  time 
step.  Both  salinity  schemes  were  considered  using  the  effective  dispersion 

coefficient  parameters  shown  in  Table  XIII-4.  Note  that  C  and  C  are  in 

x  y 

the  interval  of  experimentally  determined  values  (5.93,  20.2)  as  developed  in 
Part  V.  The  reduction  factor,  F  =  0.0388  ,  employed  equals  the  ratio  of  the 
lateral  to  longitudinal  dispersion  coefficient  as  determined  by  Elder  (1959) 
as  reported  in  Part  V.  Therefore,  the  values  of  the  calibrated  dispersion  co¬ 
efficients  are  within  the  range  of  or  are  equal  to  experimentally  determined 
values.  The  scheme  1  FCT  results  and  the  scheme  2  three-time  level  results 
ate  compared  with  measured  data  as  shown  in  Table  XIII-5.  In  most  regions  of 
the  Sound,  the  scheme  1  and  scheme  2  results  are  nearly  identical  and  are  in 
agreement  with  measured  salinity  values.  However,  in  the  vicinity  of  the  Up¬ 
per  Mobile  Bay  freshwater  inflow  the  results  diverge  as  shown  in  Table  XIII-6. 
The  scheme  1  results  are  nonnegative  and  exhibit  no  oscillations.  The  scheme  2 
results  exhibit  oscillations  behind  the  freshwater  front.  Thus  in  order  to 
study  sharp  front  problems,  scheme  1  is  recommended. 


Refined  Grid  Salinity  Simulation 


163.  The  20-24  September  period  was  simulated  on  the  Pascagoula  Channel 
refined  grid  using  the  Scheme  1  FCT  approach.  Boundary  salinity  and  elevations 


Table  XIII-4 

Calibrated  Dispersion  Coefficient  Parameters 


C  =  C  =10 
x  y 

D  =  D  =  0 
x  y 

Reduction  factor  (F)  =  0.0388 
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21.0 

22.3 
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25.2 
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T12 
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T4 
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25.6 

24.6 

T16 
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26.1 

T2 
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T44 

(49,21) 
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25.5 

T48 
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28.2 
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(57,28) 
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T62 
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T66 
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T60 
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28.  1 
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T58 
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29  .  1 

29.0 
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28.4 

28.2 

T56 

(62,32) 

29.7 

30.0 

30.3 
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29.4 

T68 

(67,26) 

27.9 
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27.5* 

28.  1 

28.2 

T70 

(71,28) 

28.4 

28.0 

ON 

On 

CM 
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28. 1 

T74 
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28.  1 

28.0 

-- 

28.4 

28.3 
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29.0 

28.5* 

26.3 

26.9 

T76 

(76,25) 

26.6 

27.0 

-- 

28.  1 

28.0 

T78 

(81,25) 

22.5 

22.0 

-- 

22.5 

22.7 

T80 

(86,25) 

22.9 

23.0 

22.6 

23.  1 

*  28  Sep  1980. 
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were  accessed  from  the  scheme  1  FCT  global  grid  simulation  for  this  period  at 
the  locations  shown  in  Table  XI-6.  Temporal  and  spatial  interpolation  were 
employed  to  develop  the  boundary  information  at  intermediate  cells  along  the 
refined  grid  boundary.  Initial  conditions  over  the  refined  grid  were  deter¬ 
mined  from  transect  data  collected  on  20  and  21  September  and  were  specified 
on  a  cell  by  cell  basis.  Since  salinity  variations  in  time  were  extremely 
small,  it  was  assumed  that  the  initial  conditions  were  representative  of  con¬ 
ditions  at  hour  2400  CST  on  20  September.  In  order  to  study  the  behavior  of  a 
purely  freshwater  inflow,  zero  values  of  salinity  were  specified  for  cells 
(8,1)  and  (16,1)  representing  the  West  Pascagoula  and  Pascagoula  River  inflows. 
Flows  during  the  period  were  taken  from  Table  VII-1  and  were  approximately 
1600  and  1000  cfs  for  the  West  Pascagoula  and  Pascagoula.  The  surge  setup 
generated  by  tropical  storm  Hermine  was  not  considered. 

164.  A  60-sec  time  step  was  employed  7200  times  in  order  to  simulate 
the  120-hr  period  beginning  hour  0000  CST  20  September.  In  order  to  allow 
hydrodynamic  transients  associated  with  initial  conditions  to  dissipate,  the 
salinity  computations  were  initiated  at  time  step  1440  (hour  2400  CST  20  Sep¬ 
tember).  Wind  information  employed  for  the  global  grid  simulation  at  6-hr  in¬ 
tervals  was  interpolated  in  time  at  each  time  step.  The  calibrated  effective 
dispersion  coefficient  parameters  in  Table  XIII-4  were  used  to  determine 
dispersion. 

165.  The  growth  of  the  freshwater  influence  is  represented  by  computing 
the  24-ppt  contour  after  48,  72,  96,  and  120  hr  shown  in  Figures  XIII-1  through 
XIII-4.  Initial  salinity  values  were  larger  than  24  ppt.  Therefore,  fresh¬ 
water  must  dilute  the  water  in  any  given  cell  to  obtain  a  cell  concentration 
less  than  or  equal  to  24  ppt.  Salinity  values  in  the  shaded  areas  in  Figures 
XIII-1  through  XIII-4  are  less  than  24  ppt.  Tidal  conditions  at  48,  72,  96, 
and  120  hr  correspond  to  a  mean  level  between  high  and  low  tides  as  shown  in 
Figure  XIII-5,  which  represents  the  water  surface  elevation  at  Pascagoula  ob¬ 
tained  during  the  global  grid  simulation.  At  72,  96,  and  120  hr,  the  tide  is 

a  lower  high  tide  in  a  semidiurnal  (neap)  period. 

166.  Figures  XIII-1  through  XIII-5  indicate  that  the  growth  of  the 
freshwater  influence  is  a  complicated  function  of  wind  loading  (wind  magnitude 
and  direction  history),  freshwater  flow  rate,  and  tidal  condition.  It  appears 
that  salinity  conditions  are  beginning  to  stabilize  at  the  end  of  this  5-day 
period.  However,  the  tidal  range  is  relatively  low  corresponding  to  a  neap 


48  HOURS 


LEGEND 

■  SALINITY  TRANSECT  STATION 
4-  SIMULATED  SALINITY  24  PPT  CONTOUR 


Figure  XIII-1.  Refined  grid,  24-ppt  salinity  contour,  48  hr 


tide  and  the  Pascagoula  River  system  inflows  correspond  to  a  late  summer  low 
flow  condition.  For  higher  tidal  ranges  and  inflows  it  may  require  a  signifi¬ 
cantly  longer  period  in  order  to  stabilize  the  salinity  distribution. 

167.  Simulated  velocity  levels  at  120  hr  (24  September  hour  24  CST)  are 
compared  with  transect  values  measured  on  24  and  25  September  in  Table  X 1 1 1  —  7 
at  locations  indicated  in  the  darkened  cells  in  Figures  XI 1 1  - 1  through  X1II-4. 
As  may  be  noted  from  these  figures,  the  freshwater  influence  does  not  reach 
any  of  the  transect  stations.  Therefore,  the  specification  of  inflow  salinity 
values  does  not  influence  the  simulation  results  at  these  locations.  Simu¬ 
lated  salinity  levels  at  the  transect  station  location  correspond  closely  to 
observed  values  and  further  verify  the  dispersion  coefficient  parameters 
employed . 

168.  In  order  to  study  the  behavior  of  refined  grid  simulation  results 
in  the  neighborhood  of  the  Pascagoula  Channel  inflow  on  a  cell  by  cell  basis, 
the  initial  salinity  levels  are  compared  with  simulated  levels  at  120  hours 
(24  September  hour  24  CST)  in  Table  XIII-8.  The  simulated  results  appear  to 
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Table  XIII-7 

Salinity  Conditions  on  the  Refined  Grid 


28  Sep  1980. 


24  Sep  1980 


T  ransect 
Stat i on 
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Grid  Cell 

20-21  Sep 
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1980 

Initial 

Condi tion 
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Computed 
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T54 

(8,22) 

27.3 

27.0 

27.6 

28.8 

T64 

(17,6) 

27.7 

28.0 

27.5 

28.2 

T62 

(24,9) 

28.5 

29.0 

26.8 

27.7 

T66 

(31,7) 

27.3 

27.0 

27.7 

27.9 

T60 

(33,7) 

28. 1 

27.0 

29.  1 

27 . 9 

T58 

(36,23) 

29. 1 

29.0 

29.6 

2,  .7 

TS6 

(34,26) 

29.7 

30.0 

30.3 

27.9 

T68 

(49,19) 

27.9 

28.0 

27.5* 

28.1 

exhibit  spurious  fluctuations  in  cells  (15-17,  2-6).  This  area  of  cells  is  lo 
cated  within  the  Pascagoula  River  itself  as  shown  in  Figure  XIII-4.  Although 
simulated  salinity  values  are  nonnegative,  the  behavior  of  the  simulation  in¬ 
dicates  that  the  refined  grid  cannot  be  used  effectively  to  predict  salinity 
levels  within  the  Pascagoula  River  itself.  The  numerical  difficulties  are  due 
to  the  averaging  process  employed  in  coupling  a  two  time  level  salinity  finite 
difference  scheme  with  a  three  time  level  hydrodynamic  scheme.  In  practice, 
these  difficulties  would  be  resolved  by  employing  an  even  more  refined  grid  in 
the  vicinity  of  the  Pascagoula  Channel  inflow. 
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PART  XIV:  STUDY  RESULTS,  CONCLUSIONS,  AND 
RECOMMENDATIONS 


169.  Salient  study  results  are  presented  and  the  following  conclusions 
are  drawn. 

a.  Twenty-two  tide  gages,  three  deep  sea  pressure  gages,  and 
twenty-one  oceanographic  moorings  at  rayed  w,th  current  meters 
at  various  depths  were  deployed  over  Mi  ssjppi  Sound  and  ad¬ 
jacent  waters  over  the  period  Apr  i  i -Nov  eii.t.  •  r  1980.  Outlaw 
(1989)  performed  a  harmonic  analvs's  r  u  *  I t  e  red  elevations 
and  currents  and  determined  tidal  ha:  .  -‘r  eties  for  nine 

constituents  (0  ,  K  ,  P  ,  M  ,  J  ,  0  ,  M  ,  S..  ,  and  N  )  . 

1  1  1  i  I  •  z  z 

Since  the  GTM  used  to  drive  tin  gloli.l  grid  produced  tidal 
charac t er i st i cs  for  the  five  major  tidal  constituents  CO  , 

Kj  ,  I' |  ,  M;  ,  and  S,),  tidal  rlrv  itions  at  11  stations  and 

currents  at  7  stations  throughout  the  study  area  were  recon¬ 
structed  using  Outlaw's  (198i)  r<  suits  using  these  live  major 
constituents  for  the  calibration  period  24-25  September  and 
verification  period  12-14  June  1 9 8c . .  Simulated  water  levels 
over  the  global  grid  matched  reconstructed  levels  to  within 
0.2  ft  at  stations  west  of  Pascagoula  and  to  within  0.1  ft  at 
Pascagoula  and  the  stations  to  the  east  for  both  periods. 
Simulated  current  components  exhibited  the  same  general  trend 
(flow  reversal  times)  as  the  reconstructed  currents  at  meters 
nearest  middepth. 

Based  upon  these  results,  the  prototype  data  collection 
program  and  subsequent  harmonic  analysis  provided  sufficient 
water  surface  elevation  and  current  data  to  calibrate  and 
verify  the  bottom  friction  mechanics  in  the  global  grid 
hydrodynami cs . 

b.  The  simulated  water  surface  elevations  on  the  global  grid  for 
the  calibration  period  20-25  September  were  used  to  provide 
boundary  elevations  for  a  refined  grid  simulation  for  the  same 
period.  Simulated  water  surface  elevations  agreed  with  the 
reconstructed  levels  to  within  0.1  ft  at  four  data  stations 
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wit  fnn  the  refined  grid.  An  L-shaped  idealization  was  em¬ 
ployed  t>  represent  tfie  branching  of  the  Pascagoula  Channel. 

H  im'iI  up*-n  these  results,  the  idealization  of  the  chan¬ 
nel  rep  i  i'm'i:  t  ed  the  .Ktual  system  in  sufficient  detail  to  ef¬ 
fective:'.  s  I’nnlite  reconstructed  levels  on  the  refined  grid. 

In  the  hvd  i  odyn.ini  1 1  simulation  outlined  in  a  and  b,  astronomic 
tide  conditions  were  considered  Meteorological  effects, 

Haim  1'’.  pies:  in>  .inomaiv.  and  long-wave  setup  must  he  con- 
s .  de  red  it  '!;*•  attempts  !i  lompare  simulated  hydrodynamic 
w.itei  level  and  currents  with  the  observed  f unfiltered)  data, 
tor  the  ill!  fit!  n  perioil,  the  long  wave  setup  due  to  tropical 
t  i  m  »!•  .''line  ;  estimated  t,  !<•  approximately  0.3  ft. 

It  was  necessary  to  increase  the  amplitude  of  the  diurnal 
(0  ,  Kj  ,  P  )  and  decrease  the  amplitudes  of  the  semidiurnal 

(M^  ,  S ; )  tidal  constituents  from  the  values  produced  by  the 
GTM  in  order  to  achieve  hydrodynamic  calibration  and  verifica¬ 
tion  on  the  global  grid  and  hydrodynamic  calibration  on  the 
refined  grid. 

Since  M2  and  S.;  were  dominant  during  the  calibration  pe¬ 
riod  and  ,  Kj  ,  and  P^  were  dominant  during  the  verifi¬ 

cation  period,  effectively  all  five  constituents  were  con¬ 
sidered  during  the  calibration  and  verification  process. 
Therefore,  the  hydrodynamics  for  any  period  may  be  effectively 
simulated  on  the  global  grid  and  refined  grid. 

A  hypothetical  regional  dredge  disposal  site  in  the  vicinity 
of  Sand  Island  studied  on  the  global  grid  appeared  to  alter 
the  tidal  pattern  only  in  lower  Mobile  Bay.  Since  the  tidal 
pattern  was  not  changed  near  the  boundary,  the  boundary  condi¬ 
tions  supplied  by  the  GTM  remain  valid  and  regional  disposal 
alterations  of  this  nature  may  be  effectively  studied  on  the 
g  loba 1  grid. 

The  hydrodynamic  effects  of  doubling  the  present  channel  width 
of  the  Pascagoula  Channel  as  studied  on  the  refined  grid  are 
localized  to  a  narrow  hand  of  cells  outlining  the  channel  and 


do  not  propagate  to  the  boundary.  As  a  result,  adjustments  to 
the  navigation  channels  (minor  adjustments  in  alignment,  deep¬ 
ening,  and  widening)  will  induce  changes  in  the  flow  pattern 
only  in  a  narrow  band  of  cells  surrounding  the  area  in  which 
the  adjustments  are  made.  Flow  changes  will  not  propagate  to 
the  boundary,  and  the  integrity  of  the  global  grid  supplied 
boundary  conditions  will  be  maintained.  This  finding  indi¬ 
cates  that  proposed  changes  in  navigation  channels  may  be 
studied  on  a  case  by  case  basis  by  developing  individual  re¬ 
fined  grids  around  the  appropriate  channels.  Boundary  condi¬ 
tions  for  each  of  these  refined  grids  may  be  developed  using 
the  global  grid  developed  in  this  study. 

Salinity  data  were  collected  at  40  stations  encompassing  Mis¬ 
sissippi  Sound  over  a  2-day  period  at  regular  intervals  over 
the  period  Apri 1-November  1980.  Conductivity-temperature 
readings  were  recorded  and  converted  to  salinity  on  a  contin¬ 
uous  basis  at  six  locations  over  the  entire  Sound  during  the 
same  period.  Due  to  biological  fouling  of  the  conductivity- 
temperature  probes,  continuous  salinity  data  could  not  be  con¬ 
sidered  reliable  and  were  not  used  to  calibrate  the  salinity 
component  of  the  model. 

As  a  result,  no  observed  time-varying  salinity  data  were 
available.  Transect  data  on  20-21  September  and  24-25  Sep¬ 
tember  were  used  to  establish  initial  conditions  and  provide 
conditions  to  compare  with  simulated  salinities.  Both  the  FCT 
scheme  1  and  full  three-time  local  explicit  scheme  2  were  con¬ 
sidered  on  the  global  grid.  Meteorological  effect  due  to  wind 
only  was  considered.  Within  Mississippi  Sound  simulation  re¬ 
sults  for  both  schemes  were  nearly  identical  and  corresponded 
Uj  observed  levels  to  within  1  ppt.  In  upper  Mobile  Bay, 
within  the  freshwater  influence,  scheme  2  results  oscillated 
behind  the  freshwater  front.  The  scheme  1  FCT  results  were 
nonnegative  and  nonosci 1 latory .  Although  due  to  the  nonavail- 
abiLity  of  time  varying  salinity  data,  the  global  grid  salinity 
calibration  is  not  as  rigorous  as  the  hydrodynamic,  the 


calibrated  dispersion  coefficient  parameters  are  within  experi 
mentally  determined  ranges. 

j  ■  The  simulated  water  surface  elevations  and  salinity  levels  on 
the  global  grid  for  the  calibration  period  20-24  September 
were  used  to  provide  boundary  conditions  for  a  refiuei  grid 
simulation  for  the  same  period.  Temporal  and  spatial  distribu 
tions  ol  freshwater  fronts  associated  with  the  Pascagoula 
River  system  were  studied  using  the  FCT  scheme  1  approach. 
Meteorological  effect  due  to  wind  only  was  considered.  During 
'•he  21-24  September  period,  the  freshwater  front  associated 
with  the  Pascagoula  River  inflow  (1000  cfs)  extended  along  the 
•■hannel  and  into  the  middle  of  the  Sound.  In  contrast,  the 
! reshwater  front  associated  with  the  West  Pascagoula  River  in¬ 
flow  (loOO  cfs)  was  confined  to  an  area  near  the  shoreline. 

This  front  appeared  to  be  driven  to  the  west  by  the  prevailing 
winds  which  were  from  the  south-southeast.  Simulated  salinity 
levels  corresponded  to  within  2.T  ppt  of  measured  transect 
values  at  all  eight  stations,  which  all  were  located  outside 
the  region  of  freshwater  influence. 

As  a  result  of  the  closeness  of  the  comparison,  the  cali¬ 
brated  values  of  effective  dispersion  coefficient  parameters 
were  substantiated  further. 

k.  The  growth  of  the  freshwater  fronts  associated  with  the  Pasca¬ 
goula  and  West  Pascagoula  Rivers  was  influenced  by  the  wind 
magnitude  and  direction,  the  magnitude  of  the  freshwater  in¬ 
flows,  and  the  astronomic  tide. 

As  a  result,  in  considering  alternatives  to  the  Pasca¬ 
goula  Channel,  care  must  be  taken  in  developing  representative 
conditions  for  these  phenomena. 

170.  The  following  recommendations  are  made  based  upon  study  results. 

a.  In  the  development  of  additional  refined  grids  for  study  of 

alternate  channel  systems,  an  implicit  time  step  in  the  neigh¬ 
borhood  of  60  sec  can  be  economically  employed  with  a  minimum 
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space  step  on  the  order  of  the  channel  width  and  a  grid  with 
1000-1500  computational  cells. 

In  studying  the  effects  on  salinity  of  dredge  practices  on 
either  the  global  or  refined  grid,  the  FCT  scheme  1  is  pre¬ 
ferred  over  scheme  2  for  sharp  front  problems.  However, 
scheme  1  is  on  the  order  of  three  times  slower  than  scheme  2. 
Thus,  if  fronts  are  not  sharp;  i.e.,  if  one  specifies  the  same 
value  of  salinity  for  the  flow  input  as  is  found  in  the  neigh¬ 
boring  cells  of  the  inflow,  scheme  2  may  be  effectively 
employed . 

Additional  research  is  warranted  to  flux  correct  the  three 
time  level  explicit  scheme  (scheme  2).  This  research  would 
produce  a  reasonably  fast  (approximately  50  percent  faster 
than  the  present  scheme  1)  but  accurate  transport  scheme  and 
eliminate  the  averaging  of  the  hydrodynamic  variables  used  in 
scheme  1  FCT. 

Due  to  the  dependence  of  the  salinity  distribution  associated 
with  freshwater  fronts  on  wind,  flow  rate,  and  astronomical 
tide,  it  is  recommended  that  additional  simulations  be  per¬ 
formed  in  order  to  develop  the  sensitivity  of  the  salinity 
distribution  with  respect  to  each  of  these  influences.  The 
results  of  this  sensitivity  analysis  would  provide  useful  in¬ 
formation  in  developing  design  scenarios  for  each  channel  al¬ 
teration  project. 
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APPENDIX  A:  GUIDE  TO  WIFM  MODEL  INPUT  REQUIREMENTS 


1 .  The  preparation  of  the  necessary  input  data  required  by  WIFM  is  con¬ 
sidered  in  light  of  the  need  to  develop  additional  refined  grids  in  the  vicin¬ 
ity  of  navigation  channels  (Gulfport  and  Biloxi).  This  appendix  is  not  in¬ 
tended  to  be  a  complete  user  guide  but  to  outline  in  detail  the  requirements 
necessary  in  developing  an  input  data  deck  for  a  refined  grid  system.  The 
basic  requirements  are  presented  in  Table  Al.  Each  item  in  this  table  is  dis¬ 
cussed  in  terms  of  WIFM-SAL  input  requirements  using  the  Pascagoula  Channel 
refined  grid  system  developed  in  this  study.  A  complete  description  of  model 
input  variables  and  required  data  formats  is  reported  by  Schmalz  (1984). 

Table  Al 

WIFM  Input  Requirements  for  a  Refined  Grid 


I.  Plot  global  subgrid  encompassing  the  refined  grid  region 

II.  Map  the  refined  grid  system 

III.  Specify  the  depth  field 

IV.  Specify  the  barrier  characteristics 

V.  Specify  flow  input  locations  and  discharges 

VI.  Specify  tidal  signal  input 

VII.  Specify  initial  salinity  conditions 

VIII.  Specify  wind  information 

IX.  Specify  simulation  control  variables 


Mapping  the  Global  Subgrid  of  Pascagoula  Region 

2.  Initial  development  of  the  refined  grid  required  reproducing  the  com¬ 
plete  global  grid  of  Mississippi  Sound. 

8.  Program  MAPIT  was  executed  using  the  table  of  input  variables  pre¬ 
viously  prepared  in  generating  the  global  grid.  Once  this  file  was  stored  a 
subgrid  of  the  region  selected  (in  this  case,  Pascagoula  Channel)  was  plotted 
at  the  same  scale  (1:40000;  as  that  chosen  for  the  refined  grid.  The  global 
subgrid  was  overlaid  on  the  project  map  (chart  11374)  for  the  refined  grid  in 
order  to  facilitate  commonality  of  feature  specification. 


Al 


4.  The  first  effort  to  map  a  grid  of  approximately  one  thousand  (1000) 
cells  having  three  (3)  bands  of  high  resolution  was  unac cep! ah  1 e . 

3 .  The  aspect  ratio  (cell  width  versus  length.)  exceeded  desirable 
limits.  Aspect  ratios  should  be  less  than  20. 

6.  To  reduce  the  aspect  ratio  and  total  number  of  it- 1  1  s  required,  the 
second  effort  eliminated  two  of  the  high  resolution  bands.  This  grid,  while 
only  10  54  cells  (47  x  22),  one  e  plotted  and  overlaid  chart  11374  was  clearly 
inadequate  (with  only  one  high  resolution  band)  to  mod  the  branching  channel 
i n t o  Pa scagoula . 

7.  The  third  effort  attempted  increasing  the  number  of  cells,  decreas¬ 
ing  the  region  to  be  mapped  and  adding  a  fourth  band  of  high  resolution.  This 
produced  a  grid  of  2058  cells  (49  x  42).  While  the  grid  would  allow  better 
idealization  of  the  channels  and  rivers  it  was  considered  too  large  for  eco¬ 
nomical  computation. 

8.  Using  this  third  grid  as  a  guide,  the  following  decisions  were  made 
with  regard  to  producing  the  final  grid: 

a.  Left  boundary  of  refined  grid  coincided  with  global  grid, 

N  =  53  . 

b.  Bottom  or  seaward  boundary  of  the  refined  grid  coincided  with 
global  grid  M  =  33  . 

c.  Three  high  resolution  bands  whose  cells,  0.1  map  inches  in  width 
(300  ft  in  nature),  were  positioned  to  achieve  an  idealized  chan 
nel.  Two  dense  bands  were  mapped  in  the  y  direction,  at  the 
main  channel  and  at  the  west  river,  and  one  band  was  mapped  in 
the  x  direction  at  the  channel  junction. 

d.  The  right  boundary  of  the  refined  grid  extended  to  approximately 
global  grid  N  =  67  . 

e.  The  upper  (land)  boundary  extended  to  approximately  global  grid 
M  =  18  . 

9.  The  final  grid  covered  an  area  from  near  Bellafontaine  Point  to  near 
Point  Aux  Chenes  Bay.  Tables  A2  and  A3  show  the  Y  ,  Y'  ,  and  a  and  X  , 

X'  ,  and  a  for  the  Y  direction  and  the  X  direction,  respectively.  These 
input  variables  for  Program  MAPIT  produced  a  grid  of  1372  points  (49  x  28). 

The  scale  is  1:40000  and  the  file  is  ordered  as  follows:  Y  direction 


followed  by  X  direction.  The  refined  grid  was  then  plotted  on  paper  and 
mylar  by  using  the  DOGRID  procedure. 

10.  The  stretfhing  coefficients  were  also  punched  using  the  DOGRID  pro¬ 
cedure  for  subsei|iien;  input  to  W1FM-SAI. 

!»  ot  li  f  ii1  hi  one i  .1  1 1  a  t  inii 

11  rin  dept  *  field  «.  i  !  ,  v  .  '  i  .'t-1  f  i  i .mi  nautical  chart  1  1') 7 4  (Ed  16  ) 
dated  Nov*  -iiher  :  o-o  >.  .  !o  ii’i.i'.ed  tin  idealized  v  ha  title  1  were  set 

t  •  dept  fi  -  18  1  pi  <  <  li  »  •  ’>  t  '  1  ■  i  '  rept  eseut  i  ng  1  and  were  set  to 

depth  *  I  1  1  tie  d  a  ’  all  w  »  ''  .  ■  i  '  i,.l  1  i  t  lilli . 

I  ‘  i  tie  •  a  *  n  viii*  ’  r  i  -  .  t  to  a.  't  eaii  Svstems  Co.  June/Julv  1980) 

♦  or  iiom  1  >.  I  an-i  f  oi.  n  •  •  -is  reviewed  to  ensure  that  depths 

asr  •  gned  to  t  to-  I  -  a  i .  •  i  !  h  ,  gB  ..  o  ells  in  t  tie  pass  would  reflect  the 

survey.  I  he  lep  t  tl  t  I  e  !  I  i  .  m  (  .  se  i  '  ■ .  l  ou;  It 

Ba  ii  hi  Spei  it  iv.it  ion 

I).  Barriers  weie  spei  i I icd  for  the  refined  grid  based  upon  the  island 
configuration  and  shallow  depth  areas.  For  harrier  specification  all  the  bar- 
riers  were  assumed  t  o  fie  exposed  throughout  simulation  and  barrier  elevations 
we  re  set  to  ♦ 1 0  ft. 

14.  Barriers  are  located  on  either  the  u  or  v  face  of  each  cell. 
Figure  A-l  shows  this  orientation. 


M,  X 

Figure  A1 .  Barrier  cell  face  orientation 


15.  To  code  the  barrier  location,  IDIR  is  set  to  1  for  a  u  face  bar- 
rier  and  set  to  2  for  a  v  face  barrier.  IDIR  indicates  flow  direction. 

16.  Further  to  locate  the  barrier  it  has  to  be  specified  by  its  refined 
grid  indices  and  their  range.  The  following  examples  illustrate: 


A. barrier  on 

the  u 

cell  faces 

of  cells 

(4, 

24)  (5, 

24) 

and  (6,  24) 

is  coded 

thus : 

ITYP 

INDX 

IDIR 

11 

12 

13 

M 

N 

a 

Nb 

1 

1 

1 

24 

4 

6 

11  is  the 

M  row  of  the 

u  face 

barrier  location. 

12  and  13  correspond  to 

the  range 

of  columns  (N's)  that 

the  barrier 

extends 

along  row  M 

(ID. 

17. 

Alternately, 

a  barrier  on  the  v 

cell  faces 

of  cell 

(7,  2)  and 

(7,  3)  is 

coded : 

ITYP 

INDX 

IDIR 

11 

12 

13 

N 

M 

a 

Mb 

1 

1 

2 

7 

2 

3 

II  is  the  N  column  of  the  u  face  barrier  and  12  and  13  correspond  to  the 
range  of  rows  that  the  barrier  extends  along  column  N  (II).  Barrier  location 
codes  are  the  first  set  of  input  variables  in  Card  Group  17. 


FI ow  Specification 


18.  Freshwater  inflow  locations  were  noted  for  the  Pascagoula  and  West 
Pascagoula  Rivers  in  terms  of  the  refined  grid  coordinate  system. 

19.  Location  information  is  specified  in  the  boundary  input  Card 
Group  17.  For  the  West  Pascagoula  the  input  is  specified  as  follows. 


Refined  Grid  Indices 


ITYP 

9 


1NDX 

1 


IDIR 

1 


M 

1 


N, 


N„ 


ITYP  =  9  for  a  flow  input.  IDIR  =  1  fur  a  vertical  oriented  x  axis  ori¬ 
ented  inflow.  INDX  specifies  the  discharge  array.  Since  this  inflow  is  the 
first  inflow  specified,  INDX  =  1  . 

20.  Note  in  simulating  salinity,  salinity  values  associated  with  the 
given  inflow  must  also  be  specified  in  Card  Group  21. 


Tidal  Signal  Input  Specification 

21.  To  specify  the  locations  for  the  tidal  signal  inputs,  where  the 
global  and  refined  grid  interface,  the  cell  centers  of  the  global  grid  cells 
on  the  refined  grid  boundaries  must  be  referenced  in  their  respective  refined 
grid  indices.  To  accomplish  this,  all  cells  in  the  global  grid  where  tidal 
inputs  would  be  made  to  the  refined  grid  were  identified  and  numbered.  Forty 
(40)  global  grid  cells  lie  along  refined  grid  (water)  boundaries.  To  iden¬ 
tify  their  centers  in  optimum  refined  grid  indices,  the  refined  grid  was 
overlaid  the  global  grid.  Table  A4  lists  the  tidal  signal  inputs  and  their 
N,  M  indices  in  both  the  global  and  refined  grid.  The  refined  grid  M  indices 
are  the  IGX  array;  the  N  indices  are  the  IGY  array  (Card  Group  3a). 

22.  The  forty  cells  in  the  table  are  assigned  boundary  values  based 
upon  the  results  obtained  at  their  corresponding  locations  in  the  global  grid. 
For  intermediate  refined  grid  points  interpolation  is  required.  To  specify 
interpolation  along  the  x-axis  for  cells  (1,  6)  through  (1,  8)  in  the  refined 
grid  the  following  data  are  required. 

II  12  13  14  15  16 

Refined  Grid  Boundary  Sweep  Tide  Inputs 

ITYF  INDX  IDIR  N  M  M£  Orientation  Srgna^  Signalg 

8  1  2  1  6  8  0  1  2 

Note  14  =  0  for  a  lower  sweep  boundary  and  IDIR  =  2  for  x  interpolation 
ITYP  =  8  for  a  tidal  elevation  boundary.  Observe  in  Table  A4  that  signals  1 
and  2  correspond  to  cells  (1,  6)  and  (1,  8)  in  the  refined  grid.  Thus  values 
at  refined  cell  (1,  7)  will  be  interpolated  from  the  values  of  tidal  signal 
inputs  1  and  2  (based  upon  x  distance). 

23.  To  specify  interpolation  along  the  y-axis  for  cells  (5,28) 
through  (8,28)  in  the  refined  grid  the  following  data  are  required. 

II  12  13  14  15  16 

Refined  Grid  Boundary  Sweep  Tide  Inputs _ 

ITYF  INDX  IDIR  M  N^  Orientation  Signal^  Signalg 


ffl 


Table  A4 


Global  Grid  Cell-Refined  Grid  Cell _ Boundary  Assignment 


Tidal 

N,M 

Indices 

|  Tidal 

N,M  Indices 

Signal  No. 

Global  Grid 

Refined  Grid 

Signal 

No. 

Global 

Grid 

Refined  Grid 

1 

53,  21 

1, 

6 

21 

67, 

29 

49, 

24 

2 

53,  22 

1, 

8 

22 

67, 

30 

49, 

25 

3 

53,  23 

1, 

9 

23 

67, 

31 

49, 

26 

4 

53,  24 

1, 

10 

24 

67, 

32 

49, 

27 

5 

53,  25 

1, 

14 

25 

67, 

33 

49, 

28 

6 

53,  26 

1, 

18 

26 

53, 

33 

1, 

28 

7 

53,  27 

1, 

20 

27 

54, 

33 

2, 

28 

8 

53,  28 

1, 

22 

28 

55, 

33 

4, 

28 

9 

53,  29 

1, 

24 

29 

56, 

33 

5, 

28 

10 

53,  30 

1, 

25 

30 

57, 

33 

8, 

28 

11 

53,  31 

1, 

26 

31 

58, 

33 

H, 

28 

12 

53,  32 

1, 

27 

32 

59, 

33 

17, 

28 

13 

53,  33 

1, 

28 

33 

60, 

33 

23, 

28 

14 

67,  22 

49, 

8 

34 

61, 

33 

26, 

28 

15 

67,  23 

49, 

9 

35 

62, 

33 

31, 

28 

16 

67,  24 

49, 

10 

36 

63, 

33 

42, 

28 

17 

67,  25 

49, 

14 

37 

64, 

33 

45, 

28 

18 

67,  26 

49, 

18 

38 

65, 

33 

47, 

28 

19 

67,  27 

49, 

20 

39 

66 , 

33 

48, 

28 

20 

67,  28 

49, 

22 

40 

67, 

33 

49, 

28 

Freshwater 
Flow  Inputs 

1 

West  Pascagoula 

57, 

19 

8, 

1 

2 

Pascagoula 

59, 

19 

16, 

1 

Reference  to  Card  Group  3a  ■ 

-  Tidal  signal  input 

locations 

1 .  Each  M 

index  of  the 

ref ined 

grid  indices  is 

an 

element  1 

of  the 

IGX 

array 

2.  Each  N 

index  of  the 

refined 

grid  indices  is 

an 

element 

of  the 

I  GY 

array 

Reference  to  Card  Group  17. 

1.  The  tidal  signal  numbers  are  15  and  16. 

2.  The  N,  M  indices  locate  where  the  tidal  signals  or  flows  are  input 
on  the  refined  grid. 


A7 


•ijtfvi-.w 
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Note  14  =  1  for  an  upper  sweep  boundary  and  IDIR  =  1  for  y  interpolation 
ITYP  =  8  for  a  tidal  boundary.  Observe  in  Table  A4  that  signals  29  and  10  cor 
respond  to  cells  (5,  28)  and  (8,  28)  in  the  refined  grid.  Thus  values  at  re¬ 
fined  cells  (6,  28)  and  (7,  28)  will  be  interpolated  from  the  values  of  tidal 
signals  1  and  2  (based  upon  y  distance). 

Salinity  Initio 1  Cond i t i on  Spec i f i cat  ion 

24.  The  calibration  period  (20-24  September  1980)  was  selected  for  simu 
lating  salinity  on  the  refined  grid. 

25.  Initial  conditions  for  each  cell  were  specified  by  transferring  the 
calibration  period  initial  salinity  values  from  the  global  subgrid  to  the  re¬ 
fined  grid.  To  facilitate  this  task,  the  global  subgrid  (scale  1:80000)  was 
reproduced  and  the  initial  salinity  values  were  recorded  thereon.  Then  the 
refined  grid  was  reproduced  at  scale  1:80000  and  overlaid  on  the  subgrid. 

26.  Initial  salinity  conditions  for  the  global  subgrid  had  been  speci¬ 
fied  based  upon  salinity  transect  surveys  taken  during  20-21  September  1980. 

27.  See  Table  A5  for  transect  station  location  summary. 

28.  Salinity  initialization  is  specified  in  Card  Group  13b. 

Specification  of  Wind  Information 

29.  The  user  has  the  option  of  considering  wind  effects  on  the  hydro¬ 
dynamics  and  the  salinity.  Wind  speeds  are  input  in  miles  per  hour  and  wind 
directions  are  in  the  meteorological  or  "from"  convention  as  opposed  to  the 
oceanographic  or  "to"  convention  (0°  from  the  N  ,  90°  from  the  E  ,  etc). 

Wind  speeds  and  directions  may  be  entered  in  tabular  form  in  Card  Group  11. 

De velopment  of  Simu lation  Control  Variables 


30.  Simulation  control  variables  determined  by  the  grid  and  depth  char¬ 
acteristics  are  DX  ,  DY  ,  and  TAU  ,  which  are  specified  in  Card  group  4. 

31.  DX  is  the  vertical  spatial  stepsize,  while  DY  is  the  horizontal 
spatial  stepsize.  These  variables  are  determined  from  the  scale  at  which  the 
grid  was  mapped.  In  this  case  the  mapping  was  at  the  1:40000  scale,  where 

1  in.  equals  40,000  in. 


;  >  i  t  hr  rr.il  sp  ue  prototype  distance  (ft)  represented  by  one 
•  e  f  *.*■  h!  in  the  vertical  or  M  direction.  Thus 


nx  =  1!>.  or  3333.33  ft 

I  1  n  .  /  f  t 


't  I  r  tin-  (<>  r  respond  i  ng  value  associated  with  the  Y  or  N  direc- 
t  i  -  >  1 1  "t  the  grid.  In  the  refined  grid  the  same  scale  was  used  for  mapping 
both  directions.  DX  eijua  1  s  DY 

14.  TAl  is  the  time-step  length  and  was  determined  in  the  following 

manne r 


At 


AS  nun  (  ft) 

s/lg(ft/sec2r--dmax(ft) 


where  AS  is  the  minimum  of  DX  times  the  smallest  X  expansion  coefficient 
(Mu)  and  DY  times  the  smallest  Y  expansion  coefficient  (Nu). 

[Mu  is  the  smallest  X  expansion  coefficient;  i.e.,  the  smallest  X  prime 

chosen  for  the  mapping  =  0.10]. 

[Nu  is  the  smallest  y  expansion  coefficient;  i.e.,  the  smallest  y  prime 

chosen  for  the  mapping  =  0.10]. 

g  is  the  acceleration  of  gravity  =  32.2  ft/sec. 

d  is  the  maximum  depth  of  water  assigned  to  a  cell  in  the  depth 

field  =  50  ft. 

At  is  the  explicit  time  step  limit. 

Thus 


At  = 


3333.33  ft*0.1 


V32. 


=  8.3  sec 


2  ft/sec2*50  ft 


This  is  the  maximum  time  step  that  can  be  used  for  explicit  modeling.  Since 
WIFM-SAI.  is  an  implicit  model,  TAU  can  be  up  to  10  times  larger  than  At 
TAU  for  the  refined  grid  was  set  to  60  sec. 

35.  Surface  elevation,  velocity,  and  salinity  stations  within  the  region 
represented  by  the  refined  grid  were  plotted  onto  nautical  chart  11374  (Ed  16). 
This  put  all  the  prototype  data  station  location  and  the  depth  information  for 
the  model  in  one  place.  Station  locations  were  then  plotted  onto  the  refined 


grid  and  global  subgrid.  Table  A5  summarizes  the  locations  and  theLr  N,  M  in¬ 
dices  in  both  grids. 

lb.  This  table  forms  the  basis  for  assignments  to  the  control  variables 
in  Card  Groups  2b,  26,  and  27  which  are  the  plotting  controls. 

NNPOT  -  number  of  surface  elevation  stations 
NVELPN  -  number  of  velocity  stations 
INPOT  -  N  indices  of  surface  elevation  stations 
IMPOT  -  M  indices  of  surfaces  elevation  station.; 


NVCORD  -  N  indices  of  velocity  stations 
MVCORD  -  M  indices  of  velocity  stations 

37.  The  value  of  NFRF.Q  (input  Card  Group  5)  can  be  determined  here. 

_  3600  sec/hr  3600  .  ,  ,  ,  ...  , 

NfRhO  =  -  or  - =  60/hr.  Thus  hydrod  v'nam ;  t  data  will  be 

Iau  sec  bO  sec 

printed  every  60  time  steps  or  once  per  hour  at  each  of  the  data  stations. 

38.  The  N,  M  indices  of  the  hydrodynamic  gages  form  the  two  arrays  NPOT 
and  MPOT  which  are  WIFM-SAL  input  Card  Group  8.  The  total  number  of  hydro- 
dynamic  prototype  gages  is  the  value  of  NGAGE  (input  Card  Group  5). 


Table  A5 

Data  Station  Location  Sununar 


Station 

Location 

N,  M 
Global  Grid 

Indices 

Refined  Grid 

Velocity 

V12 

30° 16 . 14' 

88°40 . 82 ' 

54,27 

3,20 

V13 

30° 15 . 96 ’ 

88°30.61 ' 

62,27 

35,20 

V14 

30° 13.32' 

88°32 . 40 ' 

60,30 

24,25 

V 1 5 . 

A 

30° 12.45' 

88°30 . 75 ' 

62,31 

33,26 

v15b 

30° 12.00' 

88°31 .00' 

62,32 

30,27 

V16a 

A 

30° 14.83' 

88°27 . 07 ' 

67,28 

49,23 

v16b 

30° 14. 95’ 

88°26 . 10' 

68,28 

49,22] 

Surface  Elevation 

T4  Horn  Island 

30° 14 . 1 ’ 

88°39 . 2 ' 

55,29 

4,24 

T5  Pascagoula 

30°20 . 4 ' 

88°32 .0' 

61,21 

25,6 

T6  Petit  Bois 

30° 12 . 2 ' 

88°26 . 5 ' 

67,32 

49,27 

T7  Grand  Batture 

30°20 . 5 ' 

88°24 . 4 ' 

70,21 

49, 52 

Salinity  Transect 

T54 

30° 15 . 3 ' 

88°36 . 25 ' 

8,22 

T56 

30° 12.4' 

88°30 . 7 ' 

33, 26/34, 263 

T58 

30° 14 . 8 ' 

88°30 . 45 ' 

36,23 

T60 

30° 17.05' 

88°30 . 65 ' 

32, 16/33, 16 

T62 

30° 18 . 8 ' 

88°32 . 3 ' 

24,9 

T64 

30°20 . 4 ' 

88°33 . 8 ' 

17,6 

T66 

30° 19 . 85 ' 

88°30 . 8 ' 

31,7 

T68 

30° 16 . 4 ' 

88°26 . 65 ' 

49,19 

1  . 

Station 
grid  . 

V16 

located  beyond 

n 

right  of 

cell  49,22  outside  the  refined 

2 . 

St  at i on 
grid. 

T7 

1  oca  t  ed 

beyond  ri 

ght  of  cell  49.5  outside  the  refined 

3. 

Station 

T56 

located 

at  cell 

hounda ry 

between  33,26  and  34,26. 

4. 

St  a  t ion 

T60 

1  i.cateil 

at  cell 

bounda  ry 

between  32,16  and  33,16. 

APPENDIX  B:  PROGRAM  TIDE 


1.  Program  TIDE  is  a  general  purpose  program  to  predict  the  tidal  char¬ 
acteristics  (elevation  or  current  component)  given  the  harmonic  constants  and 
prediction  period.  The  program  also  will  access  water  surface  elevation  and 
current  component  tapes  containing  the  results  of  the  standard  harmonic  analy¬ 
sis  (filtered  and  unfiltered  values).  Based  upon  the  harmonic  constants  for 
the  analyzed  stations,  the  program  predicts  the  station  tidal  characteristics 
and  outputs  the  value  next  to  the  unfiltered  and  filtered  values  at  each  hour. 
This  serves  as  a  check  on  the  harmonic  analysis.  The  program  also  will  write 
to  a  WIFM  compatible  plot  file  the  predicted  tidal  characteristic  at  user 
specified  locations.  During  the  WIFM  simulation,  this  plot  file  is  accessed 
to  enable  overplotting  of  predicted  tidal  characteristics  on  computed  results 
at  common  locations,  thereby  aiding  the  hydrodynamic  calibration  process. 

2.  The  input  variables  are  defined  followed  by  a  listing  of  the  program. 


Program  TIDE: 

1  YEAR  - 
I MONTH  ■ 
I  DAY 
[HR 

I  CONST 


Start  Time  of  the 
Prediction  Period 


NT  IDE 
NCOMP 


Year  of  the  prediction  (e.g.  1980) 

-  Month  of  the  prediction  (e.g.  9) 

-  Day  of  the  prediction  (e.g.  20) 

-  Hour  of  the  prediction  (e.g.  00) 

-  37  element  array,  specify  3  in  the  appropriate  array  element 
if  the  constituent  in  NCONST  array  is  to  be  considered  in 
developing  node  factors  and  equilibrium  arguments. 

-  Number  of  tide  stations 


Number  of  tidal  components 


1  -  elevationX 

2  -  currents  ) 


NT  -  Tape  logical  unit  number 

NCON  -  Number  of  tidal  constituents  to  include  in  the  predictions 

ICONST(J),  J  -  1 ,  NCON  -  Number  of  the  array  element  in  NCONST 
corresponding  to  each  tidal  constituent  considered. 

1  Access  harmonic  analysis  data  tape 

0  Don't  access  harmonic  analysis  data 
tape 

.ower  hour  of  prediction  period  to  write  to  WIFM  plot  file 
'pper  hour  of  prediction  period  to  write  to  WIFM  plot  file 


1EL  -  Data  Tape  Access  Switch 


NHL 

NIJH 


1WR  - 
NWK  ( I  )  , 


Number  of  stations  to  write  to  WIFM  plot  file 

1=1,  IWR  -  Number  of  the  station  to  write  to  WIFM  plot  file 
based  on  the  input  sequence  of  all  stations 


B 1 


NTT  -  Logical  unit  number  of  WIFM  plot  file 

Subroutine  TAPE:  Harmonic  Analysis  Data  Tape  Format 

TITL(J,  11,  I),  J  =  1,  8  -  Title  on  the  harmonic  analysis  data  tape  of 

the  current  station  information  file 
Note:  (11  =  1,  NST ) ,  where  NST  is  the  number  of  stations 

(1  =  1,  NR),  where  NR  =  2*  NCOMP  (filtered  and  unfiltered  files) 

JDAY  -  Start  Julian  Day 
IHR  -  Start  Hour 
IMIN  -  Start  Minute 
DELT  -  Data  Interval  (Minutes) 


r;uAT  c»  . 


v  r  y 


ITJVVhTiTT 


APPENDIX  C:  CUBIC  POLYNOMIAL  FEATHERING 


!.  I n  employing  tidal  constituent  signals  along  the  grid  boundary, 

'here  is  no  guarantee  that  at  simulation  start  time  these  signals  will  be  zero 

!n  Met,  general  this  wiil  not  bo  the  case*.  Since  water  surface  elevations 

»nd  currents  ar>  ;  o  zero  at  the  start  of  the  simulation,  there  is  a  discs. , 

:  trinity  or  impulse  at  the  boundary;  e.g.,  the  boundary  level  is  not  consistent 

:  !  *  the  :ic  !  tal  •  rid  s  t  ioit .  this  phenomenon  may  lead  :o  oscillations  in  the 

ait.-  r  i  ca  :  salat  .  :•  ,  !i  mu’,  persist  ever  several  cotnputJl  tuna !  cycles. 

,n  or  '  iv’..  I  -ms  problem,  the  t  rut  tidal  signal,  git;  , 

ms  i.r.t  f.v  ,!*,  m  ;  -behaved  signal,  lit)  ,  at  mode  1  start  and  tor  a  :mn.  • 

r  ait  i:  .  ;  •  ;  .  r.gth  At  The  we  1  I  -behaved  signal  i »  at 

’  at.  i '  •;  .  t  a  •  *  it"-  arid  at  me  '  one  1  ,  i  be  v*t  .  i -behaved  signal 

and  one  or  un.io-  d-i  i  vat  ives  equa  1  the  true  signal  and  its  corresponding  deriv- 

3  2 

at  ives  Consider  the  following  cubic  polynomial,  f (T)  =  al  +  bT  +  cf  +  d  , 


f  ( O )  =  o  ->  d  =  o 
f  (T .  )  -  aT"?  +  bT2  +  cT 

i  111 

fO  j)  -  3aT2  +  2bT]  t  c 
f "  (T  )  =  6aT  +  2b 

From  Equations  C.l  through  C.3,  we  obtain  a  ,  h  ,  and  c 

T2 

a  -  f  (T  j )  -  TjfCTj)  +  2--  f'CTjJ/T^ 
i-  =  (f’Tiy  -  6aT j)/2 
c  =  f’fT  )  -  3a T2  -  2bT] 


(C.l) 

(C..2) 

(C.3) 


( C .  4 ) 


IC.5) 


(C.6) 


Thus,  we  have  determined  a  cubic  polynomial  with  value  zero  at  t  -  0  ,  and 
whose  functional  plus  first  two  derivative  values  equal  those  of  the  true 


boundary  signal.  Required  information  for  determining  a  ,  b  ,  and  c  are 
f(Tj)  ,  f ' (tx)  ,  and  f'CTj)  . 

3.  Consider  the  true  boundary  signal  to  be  a  periodic  signal  as  shown 
in  Figure  Cl  below. 

fit) 


4.  Since  the  start  time  may  be  arbitrary,  the  vertical  axis  may  fall 
anywhere  between  1  and  5  and  four  cases  must  be  considered  as  shown  in  Fig¬ 
ure  C-2.  A  level  x  >  0  is  specified  (XLEVEL) ,  this  enables  f(T^)  to  be 
determined  in  each  case.  The  following  approximations  are  employed  for 
f ' ( T j )  and  f"(T  )  . 


f ' (Tj ) 


g(Tj  +  At)  -  g(T1  -  At) 

— 


(C.7) 


fdj) 


g(T1  +  At)  +  g(T1  -  At)  -  2g(T1) 
(At)2 


(C.8) 


Note  f (Tj)  =  gCTj)  • 

5.  The  flow  inputs  were  also  feathered  in  the  following  fashion.  The 
user  specifies  the  number  of  time  steps  (NTID)  over  which  the  feathering  is  to 
take  place.  The  feathered  signal  equals  zero  at  the  start  of  the  simulation. 

At  each  intermediate  time  step  less  than  NTID,  the  feathered  signal  is  ob¬ 
tained  by  linearly  interpolating  between  its  zero  value  at  time  zero  and  the 
true  boundary  value  at  time  step  NTID. 

6.  The  feathering  procedures  are  contained  in  Subroutine  POLY,  which  is 
listed  in  Table  Bl.  This  subroutine  is  called  if  the  user  specifies  IFETR  =  1 
on  input. 


C2 


Table  C-l 

Sul  'I'iHl  L  i  [H-  1’OI.V 

SUHR  C  U  1  t  NF  F'  PL  r 

PARARF7FR  (  NUMP  Tr  6  i  P5  »  NO  I  M  -  II  b  .  YD  I  M-b  S  . \H’  T  1  C  0  )  7  . . 

PARAMETER  < irRC  =  30b,  TO*  1  =  2  C , I PU I a  1  GO  0 ,  I T 0 Sp  1  A , JF U T  =  1 50 0 , JFl S  =  6>  ' 
LOGICAL  1  TAPE, IPT  ,1  DUMP 

cc««i  a  /  w  •'  v  -.1  i  .r  l  t  n,  j  t,  ,  j>- 7  ,v  r  ,'j> >,rw  , 

“if.  M£iT:»,!'ti:  ,  I  rn  n  .  I  r l  "  ,NC  AGi  ,«*  1  ,  •  r  , 

I  '  :  .  •  r  1  ,  j  f  r  ,  A  i  !,,  I  -  t  P  K  ,  >  y  A'  U“  !  .  ■  >  K,  «.  '•  Y  ,  r./p  ,  r  !f  . 

jr  o ,  -  r  !  s  7 ,  ks  ;  v  a  r-:>  •  ► '  ■ .  a  r  a ,  * c  1  \  v  1 .  *  r.c  ,*»pl  v, :  T*‘ r « i*  t  , 

I  C  v !  w  ■  *  I  \  T  i  r,  i  1  •  *■  1  i  *  I  T  v  f  *  !■■!  *■  «  I  ii  J  **  *  r  *  *  u  r-  r  L  'l  '  i  M't'i  Mf  ’  f  K  f»  >i 

1  1 1-  1  1  .  ,  \  (  \  T  ,  T  *  t  !  I  l  V  l  *  K  «•  *  I '  ;  T  9  \  :  »*  ♦  ■  T  ‘  w  *  \  w  A  T  «  11  A  r  %  M  ;  !  M 

*  N  S  <r  •  *  S*  /  •  .  \  l.  l  i  ‘i.  <»*■•'  r.  '  T  «  ’  v  ►  L  <  1  «•  *  1  «  I  * . 1  »  !  '  >_■  W  v.  •  \  T  A  P  u  £  * 

I*  "?.'  .  v.  - fK- •  c, ,  .  i  ^  v  ir-.  -  » r- -v »  .  \  r  a* 

,  Nf  C  •. 

C  u  K  M  '  \  *  :  /  f \;  T,  V  *1  at  %t\  ,wAt  TmP  \  a  ,  P  T  so/'"' 

f 1  s  (  v  * ;  •  '  :  *  r.  f  sf  «»r  >  *  >  ♦  *  pc  r  m  « *  o  t  a  ♦  s  T » i *  r- ■■  i ,  •  >  n  p  c «  t  .0  x  t  r  $  r , 

A  ' '  V  *  V  ;  S  » *  P  *  ♦  U  M  A  ^  0  »  7  I  M  T  •  *f“  iTPR'.  i  •  ,jr«  _  /. 

*  :  A  *.  r  *  >  '  ^  ^  •  •  M  A  <  t  ’ .  1  •  T  •  y  *  ?  a  J  Y  ♦  r  f  *  0  v  f  *  r  Y  f  ♦ .  N  1  »  7  N  2 

•  /  « Yt)*  y 1  •  *!U  ,?:■  •%:  •  T •-  *.  ♦,IL  ?  M  I  T 

•  nM»Msr*?r-T,i  *  r  >  « o '  -■  r  v ,  ^ ;  r  , .  r ,  /  : 

»►'  X  S  «  r  Y  •  4  «■*  A  T  “  S  ^  l  S  A  L  1  !  C  A  L  (  »  C  M  A  X  ,  a  «S  <  L  *  -1  1. 1  T  «  I  f  T  T  R 

C C  M « (-  /  A  A  7 /  M  K  (  I  r  K  C  *  I  f,  x  :  >  •  <  I-  p;  1 1  r  U  T  *  i  I  P  s  J  ♦  h R  G£  (  w?  U 7  f  jr  LS  ) 

t  C  1  (  I  L  7  |  J1  0$  1  t  (.  0  %  L  (  J^U7fJfLSJ  .  _ .  .  ...  _ - _ 

B  F  AD  <  \  l‘T  ^  #  1  P  P  >  X !.  r  V  r  L  i  \  T  C, 

F  0«-  M  A  7  '  f  c  -f  ,  :rj  ) 

7-B  1  NT  VO  i,£-  »M.v:  VII.  «\  TD 


POLY 

'  parms'tt 
P  A  R  H  S  V- 
_ 'C  ALAR__ 

‘  C  t  t  AR 
Pc  A'.  AP 
'  C  A  L  AO 
'PT  A  I.  4  9~ 
r  C  A  L  A  R 
*  ■'  AL  A"  _ 

'CAL  A  R 

f  C  A  |  ‘.o 
rr(|.  Af* 

SC  AL  AR  ' 
SCALAR 

L  cc  a  lap:. 

'■C  a  L  A  » 
a  C  A  L  A  P 

:  c  a  l  a  r_ 
SCALAR 
APR  A  V  •  / 

.  _-‘»rai_L 

SOLV 
cO  L  Y 

•'  C  L  Y 


c  0  P  M  A  T  (  * 

»  i  F  VC  L  MD  *,F9.",I6> 

foly 

.;V'’  s 

DO  1  IT, 

N  T  I  0 

POLY 

•M-  ID 

S  I  G  f,  --  -  1  . 

.-  POLY 

-l.  i : 

IF  (  y  L F  M  ( 

:  . I  )  .L T .SURFS  <2.  I  ) > i ! ON  A  1  , 

f  0  L  v 

2  0 
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L  ’  V  F  L 

IOLY 
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C  a  ‘  A  -  J  .  *  A  -T  T  •  IT  -i  .«  !-'*  T  T 
TT  -  D 

prim  “ju  .  i ,  i  surf  l  « .■* :  > .  j  m  .rc  :nt«i> 
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G  t.  j  ::M 

SU'A'  1  (  .  .1  )  a  TT*  (  S«1  1  .  •  T  4-  .  t.  •  )  _  . . . 

T  T  T  T  ♦  I  A  ii 

F  R  INI  ‘Ci  V.  .  ULRlRM  i.'L  1  “  "  ♦  i  > 

AORYATi*  I  !!•  A  l  Sl(,r..'.(  ♦♦!«.«•  r  I  A  T  M  f.  o  [  :  CALL'S  *  »/  * . 

(]/ 1 1  cl  f  f- .  r  ) 

CCS  T  I  Nl'F 

DC  i  r  h;  »ml0  _  ....... _ 

PR  INT  ‘ 1  P ,  I , COCURGt < Ji  1 ) 1 JAl  ,NT0« 1  ) 

FORWA  T(*  FLCW  I  r.PGT  »»IA**  ORIoIAAL  VALUtS  *t/i 

( 1  > ,  10  (Ffi  .1  ,<  x  )  ,/  n  _  ...  _  - -  - - 

TT  =  C- 

dtpfi  /(Flcat(MD)*tau> 

no  i?  j:i.Nin  _ _ 

OCHRGF«  J.  I  )  =  TT*DTP*DCHRGE  l.i,  I)' 

TT-TT*IAU 

PRINT  <J1  ?  .  I  ♦  (OCHRGl  ( J»  I )  *  dpi  «  Ml"  1  >  _ _ 1^.. _ 

FORYATL*  FLCW  INPUT  «  ♦  I A  *  •  FfATHFRCO  VALUES  *i/t 
<l>.10«Ffi.l.2X),/)> 

COLT  I  Mir  _  _____  .  .  .  _ _ 

RETURN  '  *  . . 
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APPENDIX  D:  SUBROUTINE  ADVBAR 


1.  WIFM-SAL  employs  a  cell  face  flag  convention  to  control  the  hydro- 
dynamic  computations  in  each  sweep  of  the  computational  grid.  The  cell  flag 
codes  are  stored  in  two  arrays:  ICU  for  the  x  sweep  u-face  control  and  ICV 
for  the  y-sweep  v-face  control.  The  ICU  array  consists  of  a  two-digit  pair, 

U  U  ,  while  the  ICV  array  contains  the  two-digit  V.V„  pair  for  each  cell  in 


the  grid.  For  an  open-cell  face  the  f 
second  digit  U^  or  controls  the 

evaluating  the  convective  acceleration 
equations  as  tabulated  below: 

0  -  no  advection 

1  -  x-direction  only 

2  -  y-direction  only 

3  -  both  x-  and  y-directions 

4  -  normal  in  y-direction,  approxi¬ 

mation  in  x-direction 


rst  digit  Uj  or  is  a  six.  The 

advection  approximation  employed  in 
and  eddy-dispersion  terms  in  the  motion 

5  -  normal  in  x-direction,  approxi¬ 

mation  in  y-direction 

6  -  approximation  in  x-direction  only 

7  -  approximation  in  y-direction  only 

8  -  approximation  in  both  x-  and 

y-directions 


At  the  grid  boundaries,  solid  boundaries,  and  cell-face  barriers,  no  advection 
has  been  performed  in  the  computations  for  both  the  global  and  refined  grids. 
Therefore  no  approximations  to  these  terms  have  been  made  (linearized  motion 
equation  has  been  considered)  and  codes  4-8  have  not  been  used  in  this  study. 
Subroutine  ADVBAR,  as  shown  in  Table  D-l,  determines  the  appropriate  codes  for 
U^  and  V  =  (0  ,  1  ,  2  ,  3)  for  linearizing  the  appropriate  motion  equations 
around  cell-face  barriers.  Prior  to  the  development  of  this  routine,  the 
model  user  was  required  to  specify  these  codes  through  input  for  cells  sur¬ 
rounding  cell-face  barriers.  In  the  Mississippi  Sound  global  grid  111  bariers 
(Table  VI-8)  are  employed,  while  on  the  refined  grid  25  barriers  are  employed. 
The  work  required  in  developing  the  codes  by  the  user  would  have  been  indeed 
substantial  in  this  application  and  has  been  eliminated  through  the  use  of 
Subroutine  ADVBAR. 
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